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Zusammenfassung 



Diese Arbeit ist dem Studium von Fluktuationen und Korrelationen zwi- 
schen extensiven Observablen in hadronischen Gleichgewichtssystemen gewidmet. 
Als extensive Observablen bezeichnet man MeBgroBen, die AufschluB iiber die 
„GroBe" des zu beobachtenden Systems geben, wie etwa Teilchenanzahl oder Ge- 
samtenergieinhalt, aber auch (elektrische) Nettoladung. Als Gleichgewichtssyste- 
me bezeichnet man Systeme, die sich in einem Gleichgewichtszustand befinden. 
Also in einem Zustand, der mit (intensiven) Observablen zu beschreiben ist, die 
iiber die Zeit betrachtet, ihren Wert nicht mehr andern. Hiermit konnte die Tem- 
peratur oder die Ladungsdichte des Systems gemeint sein. Unter Hadronen ver- 
steht man Teilchen, denen eine gewisse Substruktur zugeschrieben werden kann, 
die also aus elementaren Bausteinen, den Quarks und Gluonen, aufgebaut sind. 
In der uns im Alltag vertrauten Welt, in der im Wesentlichen die Hadronenvertre- 
ter Neutronen und Protonen vorkommen, spielen diese elementaren Bausteinen 
jedoch keine eigene Rolle. Zusammen mit den Elektronen aus der Familie der 
Leptonen bilden diese beiden Hadronenarten die Atome und Molekiile aus den 
wir und unsere Umgebung bestehen. 

Wenn zwei schwere Ionen (also von ihrer Elektronenhiille befreite Atomker- 
ne) kollidiert werden, entsteht fiir sehr kurze Zeit ein sehr heiBes und dichtes, 
moglicherweise kollektives System. Die Beschreibung dieses Systems ausschlieB- 
lich mit den uns vertrauten hadronischen (und leptonischen) Freiheitsgraden ist 
nun nicht mehr moglich. Die Bausteine dieser Urmaterie sind die Quarks und 
die Feldteilchen der starken Wechselwirkung, die Gluonen. Dieser Zustand kann 
jedoch im Labor nicht direkt beobachtet werden. Nach einer kurzen, aber hefti- 
gen Expansionsphase wird die Teilchendichte zu gering, als dass die Teilchen des 
Systems noch miteinander wechselwirken konnten. Quarks und Gluonen „hadro- 
nisieren" zuriick zu den uns bekannten hadronischen Freiheitsgraden. Aus den 
Fluktuationen und Korrelationen zwischen den Teilchenanzahlen verschiedener 
Arten dieser Hadronen, so wird vermutet, konnen Aussagen iiber die dynami- 
sche Entwicklung, also etwa mogliche Phasenubergange und die effektiven Frei- 
heitsgrade des System gemacht werden. Aufgrund der offensichtlich kraftvollen 
Expansion des Systems kann hier wohl kaum von einem Gleichgewichtssystem 
ausgegangen werden. 

Diese Arbeit beschaftigt sich mit trivialen Beitragen zu Fluktuationen und 
Korrelationen, wie sie durch globale Erhaltungssatze, eingeschrankte Akzeptanz 
im Impulsraum, Resonanzzerfall, oder Quantenstatistik erzeugt werden. Diese Ef- 
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fekte konnen nicht vernachlassigt werden, da deren Auswirkungen von ahnlicher 
Grofienordnung sind, wie die, die aufgrund von Phaseniibergangen oder moglicher 
kritischer Punkte im Phasendiagramm erwartet werden. Die statistischen Eigen- 
schaften von idealen, relativistischen, hadronischen Gleichgewichtsensembles wur- 
den untersucht. Neben den drei Standard-kanonischen Ensembles wurde auch 
eine Klasse von Ensembles mit endlichem thermodynamischen Bad eingefuhrt. 
Die Abhangkeit statistischer Eigenschaften auf (intensive) thermische Parameter 
wurde untersucht. Es wurde argumentiert, dass das entstehende Bild zumindest 
qualitativ auf die Schwerionen-Physik anwendbar ist. 

In Kapitel 2 wurden grofikanonische multivariate Verteilungen von extensiven 
Mefigrofien, also Verteilungen von mehr als einer Zufallsvariable, durch Fourier- 
integration der grofikanonische Zustandssumme erhalten. Eine analytische Ent- 
wicklungsmethode zur Berechnung der Verteilungsfunktionen von kanonischen 
sowie von mikrokanonischen Ensembles des idealen relativistischen Hadronen- 
Resonanz-Gases bei endlichem Vo lumen wurde vorgestellt. Die Einfuhrung der 
Temperatur in die mikrokanonischen Zustandssumme, und die von chemischen 
Potentialen in die kanonische Zustandssumme, haben dazu gefiihrt, dass die grofi- 
kanonische Zustandssumme mit der charakteristischen Funktion der damit ver- 
bundenen grofikanonischen multivariaten Wahrscheinlichkeitsverteilung identifi- 
ziert werden konnte. Mikrokanonische und kanonische Teilchenanzahlverteilun- 
gen konnten somit durch bedingte grofikanonische Wahrscheinlichkeitsverteilun- 
gen definiert werden. Unter bedingter Wahrscheinlichkeit versteht man die Wahr- 
scheinlichkeit in einem grofikanonischen Ensemble etwa eine bestimmte Teilchen- 
anzahl zu beobachten, wahrend andere extensive Mefigrofien, wie etwa globale 
Ladung oder Gesamtenergie, als fest angenommen werden. 

In Kapitel 3 wurden multivariate Verteilungen von extensiven Mefigrofien fiir 
Systeme mit endlichem - an Stelle von einem unendlichem - thermodynamischen 
Bad eingefuhrt. Hierzu wurde ein mikrokanonisches System konzeptionell in zwei 
Subsysteme aufgeteilt. Fiir diese Subsysteme wurde weiter angenommen, dass 
sie im thermodynamischen Gleichgewicht zueinander stehen. Des weiteren sollen 
sie gemeinsame Energie-, Impuls- und Ladungserhaltung respektieren. Teilchen 
konnen nur in einem der beiden Teilsysteme gemessen werden, wahrend das an- 
dere Teilsystem als thermodynamisches Bad fungiert. Wird die Grofie des ersten 
Teilsystems nun als fest angenommen, wahrend die Grofie des zweiten variiert 
werden kann, so kann man die Abhangigkeit der statistischen Eigenschaften eines 
Ensembles von dem beobachtbaren Anteil des System ermitteln. Das heifit, man 
untersucht deren Sensitivitat auf die Anwendung von globalen Erhaltungssatzen. 

Die erzeugten Ensembles sind thermodynamisch aquivalent im dem Sinne, 
dass Mittelwerte extensiver Observablen im beobachteten Teilsystem unverandert 
bleiben, wenn die Grofie des thermodynamischen Bades variiert wird, sofern das 
kombinierte System hinreichend grofi ist. Die drei Standard-kanonischen Ensem- 
bles bleiben dabei spezielle Idealisierungen von physikalischen Systemen. Die 
allgemeineren Ensembles mit endlichem thermodynamischem Bad sollten daher 
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ebenfalls von phanomenologischem oder konzeptionellem Interesse sein. Diese er- 
sten beiden Kapitel bilden somit die mathematische Grundlage fiir die analyti- 
schen Berechnungen und Monte Carlo Simulationen, die fiir diese Arbeit durch- 
gefiihrt wurden. 

Die Analyse von Hadronen-Resonanz-Gas-Ereignissen beginnt mit dem Stu- 
dium des grofikanonischen Ensembles in Kapitel 4. Das grofikanonische Ensemble 
gilt als das am leichtesten zugangliche unter den Standard-kanonischen Ensem- 
bles. Aufgrund der Annahme eines unendlichen thermodynamischen Bades, sind 
die Besetzungszahlen in den einzelnen Impulszustanden der Teilchen miteinan- 
der unkorreliert. Somit erscheinen auch die Teilchenanzahlen von je zwei ver- 
schiedenen Gruppen von Teilchensorten miteinander unkorreliert. Aufgrund der 
Annahme von unkorrelierten Besetzungszahlen, ergibt sich, dafi alle extensiven 
Mefigrofien, mit Ausnahme des Volumens von Ereignis zu Ereignis (oder von Mi- 
krozustand zu Mikrozustand) variieren. Der Energieinhalt und die Teilchenanzahl 
des Systems sind somit stark miteinander korreliert, wahrend die durchschnittli- 
che Energie pro Teilchen mit der Teilchenanzahl unkorreliert ist. Die elektrische 
Ladung, Baryonenzahl und Seltsamkeit des Systems sind korreliert mit der Hadro- 
nenanzahl, weil einige Teilchenarten mehrere dieser Ladungen tragen. Verschie- 
dene Teilchenarten haben verschiedene Quantenzahlkonfigurationen und folgen 
aufgrund ihrer Masse unterschiedlichen Impulsspektren. Die Korrelation zwischen 
Baryonenzahl und Seltsamkeit im beobachteten Subsystem hangt dann, wie die 
Korrelation von Energie und Impuls oder von Energie und Teilchenanzahl, davon 
ab, welcher Teil des Impulsspektrums der Messung zuganglich ist. 

In Kapitel 5 wurden multivariate Wahrscheinlichkeitsverteilungen von exten- 
siven Mefigrofien zu ihren mikrokanonischen Grenzwert extrapoliert. Zu diesem 
Zweck wurden iterativ Stichproben von Ereignissen erzeugt, und fiir abnehmen- 
de Grofie des thermodynamischen Bades analysiert. Die Verteilungsfunktion der 
extensiven MeBgrofien, die gewichtet wurden, konvergiert zu einer 5-Funktion, 
wahrend die Positionen der Mittelwerte konstant geblieben sind. Wahrend der 
Transversalimpuls pro Teilchen und die Teilchenanzahl im grofikanonischen En- 
semble noch unkorreliert waren, so gilt diese Annahme nicht mehr fiir Systeme mit 
endlichem Warmebad. Durch die sukzessive Konzentration auf Ereignisse in der 
unmittelbaren Nahe eines gewahlten Gleichgewichtswertes zeigen sich die Aus- 
wirkungen globaler Erhaltungssatze. Dies verlauft in einer systematischen Weise, 
so dass die Extrapolation von Observablen von ihren grofikanonischen zu ihren 
mikrokanonischen Grenzwerten moglich wurde. Ein Nachteil ist, dass die statisti- 
sche Unsicherheit, verbunden mit endlichen Stichproben, wachst, wenn man sich 
dem (mikro-)kanonischen Grenzwerten nahert. 

In Kapitel 6 wurden Teilchenanzahlffuktuationen und Korrelationen fiir ein 
neutrales Hadronen-Resonanz-Gas mit begrenzter Teilchenakzeptanz im Impuls- 
raum unter Berucksichtigung der Auswirkungen von Resonanzzerfallen unter- 
sucht. Die Extrapolationsmethode wurde angewandt, um kanonische und mikro- 
kanonische Grenzwerte von Beobachtunggrofien zu ermittlen. Ein Vergleich mit 
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analytischen, asymptotischen Losungen fiir primordiale Verteilungen in begrenz- 
ter Teilchenakzeptanz zeigen eine sehr gute Ubereinstimmung der beiden Metho- 
den. Je groBer die Anzahl der zu erhaltenen extensiven MeBgroBen, desto groBer 
ist auch die statistische Unsicherheit, die mit endlichen Stichproben von Ereig- 
nissen verbunden ist. Mikrokanonische Effekte jedoch werden durch den Monte 
Carlo Ansatz korrekt wiedergegeben. Werden Erhaltungssatze eingeschaltet, so 
werden Teilchenanzahlfluktuationen und Korrelationen modifiziert. Impulsraum- 
effekte bei Teilchenanzahlfluktuationen und Korrelationen ergeben sich aufgrund 
von Erhaltungssatzen. Fiir ein ideales, primordiales, groBkanonisches Ensemble 
in der Boltzmann-Naherung (der Ausgangspunkt), sind die Teilchenanzahlvertei- 
lungen unkorrelierte Poissonverteilungen, unabhangig von der gewahlten Akzep- 
tanz im Impulsraum, da davon ausgegangen wurde, dass Teilchen unabhangig 
voneinander produziert werden. Das Erfordernis der Energie-, Impuls- und La- 
dungserhaltungssatze fiihrt zu unterdriickten Fluktuationen und verstarkten Kor- 
relationen zwischen den Teilchenanzahlen von zwei verschiedenen Gruppen von 
Teilchen am „oberen" Ende des Impulspektrums, im Vergleich zum „unteren" En- 
de des Impulsspektrum, vorausgesetzt, einen nicht vernachlassigbarer Teil eines 
isolierten Systems wird beobachtet. Resonanzzerfalle anderen diese Trends nicht. 

Kapitel 7 ist dem mikrokanonischen Ensemble gewidmet. Ein vereinfach- 
tes physikalisches System wurde gewahlt, um eine einfachere Diskussion zu 
ermoglichen. Aufgrund der zur Verfugung stehenden analytischen Losungen, 
konnten Fermi-Dirac und Bose-Einstein-Effekte in die Analyse miteinbezogen 
werden. Bose-Einstein-Verstarkung und Fermi-Dirac-Unterdriickung der Teil- 
chenanzahlfluktuationen sind besonders stark in Impulsraumsegmenten, in denen 
die Besetzungszahlen in den einzelnen Impulszustanden groB sind. Dieser Effekt 
ist deutlich starker als der, der in friiheren Berechnungen von Fermi-Dirac und 
Bose-Einstein Effekten auf Teilchenanzahlfluktuationen im vollen Phasenraum er- 
mittelt wurde. Fiir Systeme in kollektiver Bewegung wurde festgestellt, dass die 
Rolle, die kinematische Erhaltungssatze spielen, besonders wichtig ist. Fluktua- 
tion und Korrelation von MessgroBen sind Lorentz-invariant, sofern Impulserhal- 
tung entlang der kollektiven Bewegungsrichtung beriicksichtigt wird. SchlieBlich 
wurde festgestellt, daB auch im thermodynamischen Grenzwert Langstreckenkor- 
relationen zwischen getrennten Regionen im Impulsraum verbleiben. Teilchenan- 
zahlen in unterschiedlichen Intervallen in der Rapiditat, im transversalen Impuls- 
raum, oder im Azimut, haben einen nicht verschwindenden Korrelationskoeflizi- 
enten. 

In Kapitel 8 wurde das Phasendiagramm fiir das Hadronen-Resonanz-Gas- 
Modell in seiner Abhangigkeit von Temperatur und baryonchemischem Potential 
untersucht. GroBkanonische Ladungskorrelationen und Fluktuationen sind unter- 
schiedlich in den vier verschiedenen Ecken des Phasendiagramms. Wie in der Ak- 
zeptanzanalyse ist die Korrelation zwischen zwei Ladungen stark, wenn Teilchen, 
die beide Ladungen tragen, reichlich vorhanden sind. Bei niedriger Temperatur 
und niedrigem baryonchemischem Potential dominieren Mesonen iiber Baryonen, 
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bei hoher Temperatur und hohem baryonchemischem Potential hingegen tragen 
Baryonen einen Grofiteil der Gesamtentropie des Systems. Diesem Zusammen- 
hang entsprechend verhalten sich Fluktuationen und Korrelationen der Ladun- 
gen systematisch. Auch Teilchenanzahlfluktuationen und Korrelationen im kano- 
nischen und mikrokanonischen Ensembles tragen dem Rechnung. Der Einfluss ei- 
nes bestimmten Erhaltungsatzes auf Teilchenanzahlfluktuationen ist stark, wenn 
die Teilchen der analysierten Arten reichlich vorhanden sind, und somit einen 
wesentlichen Teil der Gesamtladung und -energie tragen. Ein Vergleich zwischen 
Ensembles mit und ohne Energie-, Ladungs- oder Impulserhaltung zeigt subtile 
Unterschiede, die unter anderem am Beispiel von Resonanzzerf alien in kanoni- 
schen und mikrokanonischen Ensembles herausgearbeitet wurden. Resultierende 
primordialen Teilchenanzahlkorrelationen ergeben sich nicht aufgrund von loka- 
len Interaktionen zwischen den Bestandteilen, sondern aufgrund von global im- 
plementierten Erhaltungssatzen fur Energie und Ladungen. 

In Kapitel 9 wurden Teilchenanzahlfluktuationen und Korrelationen fiir ther- 
mische Parameter analysiert, die der chemischen Freeze-out Linie zentraler 
Schwerionenkollisionen folgen. Modellparameter wurden von friiheren Hadronen- 
Resonanz-Gas-Modell-Vergleichen mit experimentellen Messungen von mittleren 
Hadronproduktionsraten ubernommen. Diese Messungen stellen eine Verbindung 
zwischen den experimentellen Kontrollparametern Kollisionsenergie und GroBe 
der zu kollidierenden Ionen, und der Region im Phasendiagramm, die durch 
das Experiment sondiert wird, her. Einem ersten Vergleich mit experimentellen 
Daten zufolge ergibt sich eine gute Ubereinstimmung mit Hadronen-Resonanz- 
Gas-Berechnungen. Insbesondere die mikrokanonischen Formulierung des Modells 
scheint qualitative und quantitative Eigenschaften der Daten gut zu reproduzie- 
ren. 

Die hier vorgestellten Berechnungen beschreiben qualitative Auswirkungen 
auf Teilchenanzahlfluktuationen und Korrelationen. Diese Auswirkungen erge- 
ben sich allein aus den Prinzipien der statistischen Mechanik und aufgrund der 
Anwendung globaler Erhaltungssatze. Es wird interessant sein, zu sehen, ob diese 
qualitativen Auswirkungen in weiteren experimentellen Messungen von Fluktua- 
tionen und Korrelationen sichtbar werden. Wenn ja, konnten diese Effekte wohl 
von ahnlicher Grofienordnung sein, wie die Signale „neuer" Physik. Das Tren- 
nen von dynamischer Evolution des Systems und statischen Erhaltungssatzen 
ist dann eine wichtige, wenn auch nicht-triviale Aufgabe. Gleichgewichtskorre- 
lationen sind Restkorrelationen, die zuriick bleiben, nachdem das System seine 
(Entwicklungs) Geschichte „vergessen" hat. Wenn das System „beobachtet" oder 
„gemessen" wird, bevor es einen Gleichgewichtszustand erreicht hat, dann werden 
Korrelationen aufgrund der Anfangskonfiguration verbleiben. Aber selbst wenn 
sich das System wahrend seiner Lebensdauer fern von jedem Gleichgewichtspunkt 
bewegt, so spielen Korrelationen aufgrund der globalen Erhaltungssatze eine wich- 
tige Rolle. 
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Chapter 1: Introduction 



Chapter 1 
Introduction 

This thesis is dedictaed to the study of fluctuation and correlation observ- 
ables of hadronic equilibrium systems. The statistical hadronization model of 
high energy physics, in its ideal, i.e. non-interacting, gas approximation will be 
investigated in different ensemble formulations. The hypothesis of thermal and 
chemical equilibrium in high energy interaction will be tested against qualitative 
and quantitative predictions. 

The statistical hadronization model, first introduced by Fermi [1] and Hage- 
dorn [2], has been surprisingly successful during the last couple of decades in 
describing fundamental properties of systems created in heavy ion collisions, cos- 
mic rays, and elementary particle reactions. In the context of heavy ion collisions 
it has been applied to an extensive set of data on hadron production, ranging 
form the center of mass energies of the experiments at the SIS [3, 4, 5], AGS [6], 
SPS [7, 8, 9], and most recently, RHIC [10] facilities. Model predictions for the 
upcoming LHC and future FAIR [11, 12, 13, 14, 15, 16, 17] experiments largely 
follow these trends. A systematic evolution of thermodynamic parameters, as 
collision energy (and size of colliding ions) is changed [18, 19, 20, 21, 22], has 
allowed to establish the 'chemical freeze-out line', which is now a commonly ac- 
cepted ingredient in the phase diagram of strongly interacting matter. More 
controversially this model has also been applied to a range of elementary collision 
systems [23, 24, 25], where only few particles are produced, and the picture of a 
gas of hadrons can hardly be suitable. The remarkable ability of the statistical 
model to explain these data has lead to the suggestion [23, 24, 25, 26, 27, 28], 
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that thermal (or phase space dominated) particle production is a general prop- 
erty of the hadronization process itself, rather than the result of a long sequence 
of microscopic interactions. This thesis will not argue about possible physical 
interpretations [29] of the partition function of statistical mechanics. It is, how- 
ever, noted that in order to apply a semi-classical approximation a volume of 
0(10 fm 3 ) seems to be sufficient [30]. 

Somewhere above this chemical freeze-out line in the phase diagram a phase 
transition from hadronic degrees of freedom to a phase of deconfmed quarks and 
gluons, generally termed the quark-gluon plasma is conjectured; more specifically, 
a first order phase transition at low temperature and high baryon chemical po- 
tential, and a cross-over at high temperature and low baryon chemical potential. 
In between, a second order endpoint or a critical point might emerge. For recent 
reviews see [31, 32, 33, 34, 35, 36, 37, 38]. One of the answers still outstanding 
in high energy physics is then the one of a possible formation of a deconfmed 
state of matter, and the nature of the transition between phases. The growing 
interest in the study of event-by-event fluctuations in strong interactions is, thus, 
motivated by expectations of anomalies in the vicinity of the onset of deconfine- 
ment [39, 40, 41, 42, 43] and in the case when the expanding system goes through 
the transition line between quark-gluon plasma and hadron gas [44, 45, 46]. In 
particular, a critical point of strongly interacting matter may be accompanied by 
a characteristic power-law pattern in fluctuations [47, 48, 49]. Recently, it has 
been suggested that correlations across a large interval of rapidity could also arise 
from color glass condensate initial conditions [50, 51]. 

In recent years a wide range of experimental measurements of fluctua- 
tions of particle multiplicities [52, 53, 54, 55, 56, 57, 58, 59], transverse mo- 
menta [60, 61, 62, 63, 64] and multiplicity correlations in rapidity [65, 66, 67] 
have been reported, leading to a lively discussion regarding their physical inter- 
pretation. The most promising region in the phase diagram for observation of 
critical phenomena seems to be accessible to the SPS accelerator [68, 69, 70]. A 
new SPS scan program [71] for different ion sizes as well as center-of-mass energies 
has been proposed to study strongly interacting systems at different energy and 
net-baryon densities, and life times. Also a second scan program [72], lowering 
the RHIC colliding beam energy to probe the same domain of the phase diagram, 
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is under discussion. This should be as well the main motivation for further in- 
vestigation of properties of statistical ensembles. 

Fluctuations of, and correlations between, various experimental observables 
are believed to have the potential to reveal new physics. They are amongst the 
most promising candidates suggested to be suitable for signaling the formation 
of new states of matter, and transitions between them. For recent reviews here 
see [41, 42, 43, 44, 46, 47, 48, 49, 73, 74]. In particular, multiplicity and charge 
fluctuations have been proposed to be a good discriminating tool between quark- 
gluon plasma and hadron gas [39, 40]. Provided the signal survives the phase 
transition [75], and subsequent evolution of the system. Hence, in order to prop- 
erly assess the discriminating power of such observables, one might firstly want 
to asses the magnitude of 'trivial' physical effects, such as the ones induced by 
global conservation laws, quantum statistics, resonance decays, kinematical cuts, 
finite spatial extension, etc. The statistical properties of a sample of events are 
certainly not solely determined by critical phenomena. To get a reliable indica- 
tion of new physics, it is therefore important to note that most fluctuation and 
correlation observables are also sensitive to some 'baseline' contributions that, 
nevertheless, can have non-trivial behavior. 

A rather general observation regarding statistical samples can be made: The 
statistical properties of a sample of events depend on the rules chosen to select 
this sample from an even larger sample, and on the degree of completeness of 
the information available about the sample. In the context of heavy ion collision 
physics these two aspects roughly translate into centrality class construction and 
particle acceptance in momentum space. 

Centrality selection is discussed first. Two heavy ions collide with relativis- 
tic momenta. Being extended objects, they can do so in many different ways. 
Roughly the following rule should apply (in the average sense): The larger the 
interaction region, the more particles of each particle species are produced. The 
problem is now that the initial state of a collision cannot be observed directly. 
All that can be observed is its final state. From this, one can then infer the 
likelihood of a certain initial state. Yet, each single possible final state observ- 
able will generally suggest a slightly different initial state. Hence, the need to 
average over centrality classes. The problem is then, within any such centrality 
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class will be events with rather different initial states, altering the true correlation 
between two observables. Results will depend, in general, on which trigger was 
chosen to construct the centrality class or sample of events [76, 77]. Centrality 
selection, although being a crucial experimental issue, is foreign to the statistical 
hadronization model, and will not be considered in this thesis. 

The second problem is particle acceptance. This term should denote particle 
identification and momentum measurement. Centrality selection is now ignored 
and a perfectly prepared initial state is assumed. Two limiting cases can be 
explored. The first one being the ideal detector. All final state particles are ob- 
served. Any correlation is measured to any degree. The opposite limit would be 
a very bad detector. Capable of only detecting a particle every once in a while. 
This detector could surely measure the ratios of the occurrence of particles of dif- 
ferent species. But it would be completely unable to inform one on how particle 
species within one event are correlated. Any realistic detector is in between these 
limits. And, hence, will disagree with either limit, in basically any observable. 
The statistical hadronization model provides a natural framework for such a dis- 
cussion. 

Strictly speaking the model, at least in the from presented here, is not di- 
rectly comparable to data of heavy ion collisions. It does certainly no justice to 
the complexity and the dynamical evolution of the system it seeks to describe. 
The model exhibits no collective flow, or expansion. The size (or volume) of the 
system is assumed to be the same for all events. Furthermore, no distinction is 
made between chemical and kinetic freeze-out. On the other hand, the medium 
created during the collisions of two heavy ions is rapidly expanding, while the 
initial state can only be accessed indirectly. Yet, essentially the model seeks to 
describe properties of many-particle systems, and, here in this thesis, their sta- 
tistical properties. The statistical hadronization model, might not be a bad place 
to start such a task. 

The purpose of this thesis is the calculation of 'baseline' contributions, 
on top of which one hopes to find unambiguous signals of a phase transi- 
tion [44, 45, 46], a critical point [47, 48, 49], or thermal/chemical (local or global) 
non-equilibrium [78, 79, 80]. I.e. to study these baseline correlations in a lim- 
iting case: that of a thermalized relativistic ideal (no inter-particle interactions) 
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quantum gas, for which we want to assess the importance of globally applied 
conservation laws, quantum statistics, resonance decays, and kinematical cuts for 
fluctuation and correlation observables. In this case, all observables are calcula- 
ble simply using statistical mechanics techniques. Such an approach has a long 
and distinguished history of calculating particle multiplicities in hadronic colli- 
sions [1, 9, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90]. Given its success in describing 
experimentally measured average hadron yields, and its ability to reproduce low 
temperature lattice susceptibilities [91], the question arises as to whether fluctu- 
ation and correlation observables also follow its main line. Critical phenomena 
(and many more), however, remain beyond the present study. 

Conventionally in statistical mechanics three standard ensembles are dis- 
cussed; the micro canonical ensemble (MCE), the canonical ensemble (CE), and 
the grand canonical ensemble (GCE). In the MCE 1 one considers an ensemble of 
micro states with exactly fixed values of extensive conserved quantities (energy, 
momentum, electric charge, etc.), with 'a priori equal probabilities' of all micro 
states (see e.g. [92]). The CE introduces the concept of temperature by intro- 
duction of an infinite thermal bath, which can exchange energy (and momentum) 
with the system. The GCE introduces further chemical potentials by attaching 
the system under consideration to an infinite charge bath 2 . Only if the experi- 
mentally accessible system is just a small fraction of the total, and all parts have 
had the opportunity to mutually equilibrate, can the appropriate ensemble be 
the grand canonical ensemble. 

The main focus of the past study of the statistical hadronization model has 
been on the mean multiplicities of produced hadrons. However, there is a quali- 
tative difference in the properties of mean values and event-by-event fluctuations 
about these mean values in statistical mechanics. In the case of the ensemble 
averages, results obtained in the GCE, CE, and MCE approach each other in the 
large volume limit. One refers here to as the thermodynamical equivalence of 
statistical ensembles. However, even in this limit, these ensembles have different 

lr The term MCE is also often applied to ensembles with energy but not momentum conser- 
vation. 

2 Note that a system with many charges can have some charges described via the CE and 
others via the GCE. 
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properties with respect to fluctuations and correlations [93, 94, 95]. In the MCE, 
energy and charge are exactly fixed. In the CE, charge remains fixed, while energy 
is allowed to fluctuate about some average value. Finally, in the GCE the require- 
ment of exact charge conservation is dropped, too. One may also consider isobaric 
ensembles [96], or even more general 'extended Gaussian ensembles' [97, 98]. In 
previous articles [93, 94, 95, 96, 98, 99, 100, 101, 102, 103, 104, 105, 106] it was 
shown that these differences mean that, in particular, multiplicity fluctuations 
and correlations are ultimately ensemble specific. I.e. depend on how the system 
under investigation is prepared. 

The observation that fluctuations of certain quantities are constraint by con- 
servation laws is not new, nor restricted to heavy ion physics. The Italian physicist 
Ugo Fano wrote back in 1947 in the abstract of his article [107] 'Ionization Yield 
of Radiations. II. The Fluctuations of the Number of Ions': 

The ionization produced by individual fast charged particles is frequently 
used as a measure of their initial energy; fluctuation effects set a theo- 
retical limit to the accuracy of this method. Formulas are derived here 
to estimate the statistical fluctuations of the number of ions produced by 
constant amounts of radiation energy. The variance of the number of 
ionizations is found to be two or three times smaller than if this number 
were governed by a Poisson distribution. An improved understanding is 
gained of the statistical treatment of fluctuation phenomena. 

So, it has been well understood, already some 60 years ago, that global con- 
servation laws affect the statistical properties of a system. Similar observation 
was made in the fluctuation of the number of atoms forming a Bose-Einstein 
condensate. The 'grand canonical fluctuation catastrophe' at the Bose-Einstein 
condensation point [108] is avoided by a micro canonical formulation [109]. 

This thesis is organized as follows: In Chapter 2 grand canonical joint distri- 
butions of extensive quantities are obtained by Fourier integration of the grand 
canonical partition function. An analytical expansion method for calculation of 
distributions at finite volume for the canonical as well as the micro canonical 
ensembles of the ideal relativistic hadron resonance gas will be presented. In 
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Chapter 3 joint distributions of extensive quantities are then considered for sta- 
tistical systems with finite, rather than infinite, thermodynamic bath. These 
two chapters form the mathematical basis for analytical calculations and Monte 
Carlo simulations performed in this thesis. The analysis of hadron resonance gas 
events will start with the study of the grand canonical ensemble in Chapter 4. 
The grand canonical ensemble is considered to be the most accessible amongst 
the standard canonical ensembles. In Chapter 5 joint distributions of extensive 
quantities will be extrapolated to the micro canonical limit. In Chapter 6 mul- 
tiplicity fluctuations and correlations are studied for a neutral and static hadron 
resonance gas with limited acceptance in momentum space. The effects of res- 
onance decay will be considered, and the extrapolation scheme will be applied 
to obtain canonical and micro canonical ensemble limits. Chapter 7 is dedicated 
to the micro canonical ensemble itself. A simplified physical system is chosen 
to allow for smoother discussion. Owing to available analytical solutions, Fermi- 
Dirac and Bose-Einstein effects, as well as collective motion, are included into 
the analysis. In Chapter 8 the temperature and baryon chemical potential phase 
diagram of the hadron resonance gas model is explored. Lastly, in Chapter 9, 
multiplicity fluctuations and correlations will be analyzed for thermal parameter 
sets following the chemical freeze-out line. A first comparison to available exper- 
imental data suggests good agreement with hadrons resonance gas calculations. 
A summary, Chapter 10, will close the thesis. Further technical details of the 
calculations and Monte Carlo simulations, are presented in the Appendix. 
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Chapter 2 

Grand Canonical Partition 
Function 

Aim of this chapter is to introduce a technique for calculation of grand canon- 
ical probability distributions 1 of various extensive quantities at finite volume. 
The method is based on Fourier analysis of the grand canonical partition func- 
tion. Taylor expansion of the generating function is used to separate contributions 
to the partition function in their power in volume. Laplace's asymptotic expan- 
sion is employed to show that any equilibrium distribution of multiplicity, charge, 
energy, etc. tends to a multivariate normal distribution in the thermodynamic 
limit. Gram-Charlier expansion allows additionally for calculation of finite vol- 
ume corrections. Analytical formulas facilitate inclusion of resonance decay (in 
full acceptance), or inclusion of finite acceptance effects (without resonance decay) 
directly into the system partition function. Multiplicity distributions in (micro) 
canonical ensembles can then be defined through conditional grand canonical 
distributions. 

The discussion of equilibrium systems with finite, rather than infinite, ther- 
modynamic bath, in Chapter 3, will provide a description for a Monte Carlo 
approach, capable of assessing the effects resonance decay in finite acceptance. 
These two chapters will prepare the mathematical framework for the following 
discussions of statistical fluctuations and correlations in hadronic equilibrium sys- 
tems. 

■"■The term 'distribution' is used for simplicity, rather than the more correct 'probability 
density function'. 
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2.1 Probability Distributions 

In textbooks on statistical mechanics (see e.g., Ref. [110, 111, 112]) often first 
the MCE is introduced, where exact conservation laws for energy-momentum 
and particle number are imposed on a collection of micro states. Relaxing the 
constraints for energy and momentum constitutes the CE, while allowing addi- 
tionally particle number to fluctuate about some mean value introduces the GCE. 
In a relativistic gas of hadrons quantum numbers (charges), rather than particle 
numbers, will be the conserved quantities. 

Here it will prove to be of considerable advantage to start off with the GCE for- 
mulation and imposing exact conservation laws thereafter. Generally, the (micro) 
canonical partition function is obtained from the grand canonical one by multipli- 
cation with Kronecker (or Dirac) 5-functions which pick out a set of micro states 
consistent with a particular conservation law. It is often more economical to use 
Fourier representations of ^-functions, rather than the 5-function themselves. 

The basic idea is to define the probability of finding the system in a state 
with a given number of particles Na of some species A at some fixed value of 
conserved (electric) charge, Q, i.e. the CE distribution P ce (NA), in terms of 
the GCE distributions, P gC e(NA,Q) and P gce (Q)- In general one may write for 
the multiplicity distribution P C£ (Na) of a CE with conserved electric charge Q: 

number of all micro states with Q and Na 

Pce(N A ) = — — — 1 . 2.1 

number of all micro states with Q 

Likewise, one can write for the CE joint multiplicity distribution P ce (N a, N b) of 
particle species A and B: 

P (N N ) number of all micro states with Q, Na and Nb ^ 2) 

number of all micro states with Q 

The number of all micro states with electric charge Q, and multiplicities Na 
and Nb of a system with temperature T = and volume V is given by the CE 
partition function Z(V, (3, Q, Na, Nb). Similarly, Z(V,(3,Q) denotes the num- 
ber of micro states with fixed electric charge Q, but arbitrary multiplicities Na 
and Nb, for the same physical system. 

Throughout this thesis the following conventions are used: Ensembles are to 
be identified by the arguments of their partition functions. Extensive quantities, 
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denoted by capital letters in the argument, are to be considered as exactly con- 
served. Intensive quantities, denoted by small letters in the argument, indicate 
that the corresponding extensive quantity is conserved only in the average sense. 
The partition function Z(V, E, Q) denotes a MCE with conserved energy E and 
electric charge Q. The partition function Z(V, (3, Q) belongs the a CE with an 
infinite heat bath at temperature The energy content of the volume V 

fluctuates then about some mean value (E). Lastly, Z(V,j3,fi), is the number 
of micro states available to a system with an infinite heat and charge bath at 
temperature and chemical potential \i. All constraints on the micro states 
are dropped, and a weighted average over all charge and energy configurations is 
taken. 

Particle number is here included in the argument of partition functions, de- 
spite the fact that no (non- vanishing) particle number specific chemical potential 
is introduced. Fourier integrals associated with particle multiplicity are going 
to be solved along with those associated with conserved quantities. In the next 
chapter, particle number is deligated to an index of the partition function, as its 
integrals do not need to be solved for the Monte Carlo approach. 

The strategy to calculate joint multiplicity distributions could thus be the 
following (in principle also valid at finite volume): 

p (N N) Z(y,p,Q,N A ,N B ) 

PUN A ,N B ) - z{ViM) > (2-3) 

Z(V,P,Q,N A ,N B ) Z(V,(3,fi) 
Z(V,P,u.) e Q»P Z(V,(3,Q)' [ - } 

= P gce (Q, Na, Nb) Pgcl(Q) = P a ce(N A ,N B \Q) . (2.5) 

In order to get from Eq.(2.3) to Eq.(2.5) both canonical partition func- 
tions Z(V, f3,Q, N A , N B ) and Z(V,/3,Q) are divided by their GCE counter- 
part Z(V,/3,/j,) and multiplied by the Boltzmann weig ht e Q ^. The first 
term on the right hand side of Eq.(2.4) then equals the GCE joint distribu- 
tion Pgce{Q, Na, N b ), while the second term is just the inverse of the GCE charge 
distribution P gce (Q). Their ratio is the (normalized) GCE conditional distribution 
of particle multiplicities Na and N B at fixed electric charge Q, P gce (N A , N B \Q), 
and equals the CE distribution P ce (N A , N B ) at the same value of Q. This result 
is independent of the choice of chemical potential /i. 

The problem of finding a solution, or a (large volume) approximation, to 
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the CE distribution P C£ (Na,N b ) is now turned into the problem of finding a 
solution or approximation to the GCE distribution of multiplicities N A and N B , 
and charge Q, P gce (Q, N A , N B ). 

It is worth noting that Eqs. (2.3-2.5) £1X6 cLS well the basis for any Monte 
Carlo approach [113, 114]. A sampling distribution, usually taken from a Boltz- 
mann GCE system, is used to generate a {iVj}-tuple of particle multiplicities of 
all species % considered. All 'events' consistent with certain constraints, like a 
set of conserved charges, are accepted, while the rest is rejected. On the basis of 
this set of all accepted 'events' one constructs an ensemble by using a suitable 
re-weighting scheme to account for quantum statistics and proper normalization. 
For the Monte Carlo approach presented in the next chapter an unconstrained, 
i.e. grand canonical, sample of events, P gce (Q, {A^}), is generated. Constrained 
distributions could then be defined through: 

P(Q,{Ni}) = W(Q) P gce (Q,{N t }) , (2.6) 

where a calculated weight factor W(Q) is employed to project out a set of events 
in the vicinity of an equilibrium value Q eq . 

An immediate consequence of Eqs. (2.3-2.5) is that temperature and chemical 
potentials appear in this formulation of (micro) canonical distributions (as well 
as in the Monte Carlo [113, 114], microscopic correlator [99] and saddle point ex- 
pansion [115] approaches). At first sight this seems to be a serious problem and 
an unnecessary complication of the initial task of finding a reasonable approxima- 
tion to CE and MCE partition functions. However, the main technical challenge 
when numerically integrating the original version of the (micro) canonical parti- 
tion function arises from a heavily oscillating integrand. Auxiliary parameters /3 
and n will produce a very smooth function, for which approximation schemes can 
be used. In taking the ratio, Eq.(2.4), artificially introduced temperature and 
chemical potential drop out. The quality of the approximation on the other hand 
will crucially depend on their choice. 

In Section 2.2 the generating function of the GCE distribution of extensive 
quantities is introduced. In Section 2.3 it is shown that the requirement of max- 
imizing the generating function of the charge distribution at some given equilib- 
rium point leads to a unique determination of thermal parameters, and moreover 
constitutes the optimal choice for the approximation scheme. Emerging thermo- 
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dynamic relations are discussed in Section 2.4. 



2.2 Generating Function of the Charge Distri- 
bution 

The GCE partition function of an ideal relativistic gas with volume V, local 
temperature T = chemical potentials fij and collective four velocity reads 
(the system four-temperature [113] is /3 M = (iu^): 



Z(V,(3,u^,fij) = exp 



VV (fru^fij) 



(2.7) 



where ^ ((3, u^, fij) is a sum over the single particle partition functions ipi (/3, u^, fij) 
of all particle species i considered in the model: 

^ {13, u^, fij) = ^ ^i(/3>«A»>A*j) • ( 2 - 8 ) 

i 

The single particle partition function ipi ((3, u^, fij) of particle species i is given 
by a Jiittner distribution: 

^ m m , fij) = ^3 |d 3 p In (l ± e ^ + * * , (2.9) 

where are the components of the four momentum, q? are the components of 
the charge vector, and ^ is the degeneracy factor. The upper sign refers to Fermi- 
Dirac statistics, while the lower sign refers to Bose-Einstein statistics. The case 
of Maxwell-Boltzmann statistics is analogous. 

For a hadron resonance gas, including hadrons and their resonances up to 
a mass of about 2.5 GeV, i.e. excluding the charm quark sector, the vector 
of chemical potentials fij and the 'charge' vector q\ of particle species % are 
introduced 2 : 

fij = (fi B , fi s , fi Q , fi NA , fi NB ) q{ = (b h Si, q u n A (Q), n B {Q)) , (2.10) 



2 Finitc acceptance effects are discussed here. Modifications to Eqs. (2.9,2. 10) which allow for 
inclusion of resonance decay in full acceptance are presented in Appendix C. 
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where /is, and /iq are the baryon, strangeness, and electric charge chemi- 
cal potentials, respectively. The parameters u. Na and /ijv s are particle-specific 
chemical potentials, and could denote out of chemical equilibrium multiplicities 
of species 'A' and 'B', similar to phase space occupancy factors 7 S [116] and 7 g 
[117, 118]. Throughout this thesis out of equilibrium effects are neglected, and 
thus flN A = fl n b = 0. 

In addition, 6 i; Sj, and g« are the baryonic charge, the strangeness, and the 
electric charge of a particle of species i. Q is the momentum space bin in which 
one is set to measure particle multiplicity. n^(fi) = 1 if the momentum vec- 
tor of the particle is within the acceptance, ua{^) = if not. The charge 
vector q\ also contains, to maintain a common notation for all particle species 
considered in Eq.(2.8), the 'quantum' number hb(O). If one were set to deter- 
mine the joint distribution of positively versus negatively charged hadrons, for 
instance, the 7r + particle would have q n + = (0,0,1,1,0), while the K~ particle 
would have qx- = (0, —1, —1, 0, 1), see Chapters 6 and 7. 

One may also be interested in correlations of, for instance, the systems net- 
baryon number B and net-strangeness S, as e.g. in Refs.[74, 91]. In this case, 
the A particle, with q\ = (1? —1) 0, 1, —1), would be counted in groups A and B, 
provided the momentum vector is within the acceptance Q. The S~ particle on the 
other hand carries two strange quarks and would have = (1, —2, —1,1, —2), 
see Chapters 4 and 5. 

The generating function of the charge distribution in the GCE is introduced 
by the substitutions in Eq.(2.9): 



/3 fij ->■ (3 fij + i(f)j , 



(2.11) 
(2.12) 



The yet un- normalized joint probability distribution of extensive quantities Q\ P M 
in the GCE is then given by the Fourier transform of Eq.(2.7) after substitutions 
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Eqs.(2. 11,2.12): 



oo 



z 



d J <f) 



j 




(27T) 



4 



— 7T 



— OO 



x exp W (P,^, % ap, <j>j) 



(2.13) 



Depending on the system under consideration, the vector of extensive quanti- 
ties Q J and corresponding Wick rotated fugacities 4>j could read: 



Q j = (B, S, Q,N A ,N B ) fa = (cf> B , <f> s , cf> Q , <p NA , <p NB ) . (2.14) 



Here B is the net-baryon number, S is the net-strangeness, and Q is the elec- 
tric net-charge of the system. Together with particle numbers Na and Nb this 
would be a 5-dimensional distribution in the case of a CE hadron resonance gas. 
Additionally for four-momentum conservation, yielding a 9-dimensional Fourier 
transform Eq.(2.13) for a MCE hadron resonance gas, the vectors P M and are 
introduced: 



where E is the energy and P x , P y , and P z are the components of the collective 
momentum of the system, while a M are the corresponding fugacities. 

The distinction between discrete (Kronecker 5) and continuous quantities 
(Dirac 5) is not relevant for the large volume approximation, where particle num- 
ber is a continuous variable to be integrated over. Proceeding by Taylor expansion 
of Eq.(2.8), it is convenient to include discrete and continuous quantities into a 
common vector notation: 




a 




(2.15) 





(2.16) 



The dimensionality of the vector Q l is denoted as I = 2 + 3 + 4 = 9 for 
a MCE hadron resonance gas. Now, expanding the cumulant generating func- 
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tion, \l> (j3, u^fj.j] 9[), in a Taylor series yields: 

oo . n 

* (Av^O - E ^ K « ,b, *" ,In ^ V- A . ( 2 - 17 ) 



72! 
n=0 



where the elements of the cumulant tensor, defined by: 



de h de h ...de l , 



(2.18) 

!=0( 



Generally cumulants are tensors of dimension L and order n. The first cumulant 
is then a vector, while the second cumulant is a symmetric L x L matrix. A good 
approximation to Eq.(2.13) around the point Q l eq = (Q 3 eq , Peq), can be found in 
terms of a Taylor expansion of Eq.(2.8) in Q\ = (<f)j, a^), if (see Section 2.3): 



dQ} 



0, . (2.19) 



e <y 



Implicitly, Eq.(2.19) does not define chemical potentials fij and four-temperature 
Pn = fiiin, but corresponding Lagrange multipliers, which maximize the ampli- 
tude of the Fourier spectrum Z Ql (V, fij) of the generating function for a 
desired value of QjL = (Q 3 eq ,Pe„)- Their values generally differ from the GCE 
set {fi jU^^/jLj), however they coincide in the thermodynamic limit. Lagrange 
multipliers are relevant for finite volume corrections in Section 2.3. The tempera- 
ture (3 and the four-velocity are not five independent Lagrange multipliers, as 
the four- velocity vector has unit length u^u^ = 1. The integrand of Eq.(2.13) is 
sharply peaked at the origin <pj = ot^ = in the thermodynamic limit. The main 
contribution therefore comes from a very small region [115]. To see this, a second 
derivative test can be done on the integrand of Eq.(2.13) taking into account the 
first two terms of Eq.(2.18), see Appendix B. 

The cumulant of th order, kq, is just the logarithm of the GCE partition func- 
tion divided by the volume, Z (V, f3,u^, fij) = exp(V^fi;o)- Hence, after extending 
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the limits of integration to ±00, which will introduce a negligible error, one finds: 



d L e 



exp 



iQ l 0i 



7 



n=l 



(2.20) 



It is worth noting that the parts of the integrand of Eq.(2.13) related to discrete 
quantities were 27r-periodic, while the integrand of Eq.(2.20) is a superposition of 
oscillating and decaying modes. Spelling out the first two terms of the summation 
yields: 



Z Q (V,(3,u„^) ~ Z(V,f3,u„ H ) 



—00 
h,h 



d L e 



exp 



-iQ l e l 



(2.21; 



n=3 



I 



Performing now a change of variables will simplify this integral: 



= Wa t k e k , 



(2.22) 



where a t k is the square root of the second rank tensor K2: 



(J, 



1/2 



(2.23) 



The new integral measure d Lq 9 then equals to: 



dH = detlvWl d L 9 = V L/2 detlcrl d L 6 



(2.24) 



Lastly in terms of this transformation normalized cumulant tensors A n with com- 
ponents: 



A 



ll,l<2,...,ln 



K 



fel,fc2,.,fen f„-l\ h f-l\ I 



(»-v (--v-.-k 1 ) 



]\ In 



(2.25) 



are introduced. The new variable £ will be a measure for the distance of the 
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actual charge vector Q k to the peak Vk\ of the distribution: 

i l = {&-VK*)(<r-\ l V- 1 0. 
Including above steps at once yields: 



(2.26) 



V L / 2 det\a\ J (2tt) l 



exp 



A 



oo 

ll,la,...,ln 



^— ' n! 



n=3 



(2.27) 



Eq.(2.27) is the starting point for obtaining an asymptotic solution in this sec- 
tion, as well as for finite volume corrections in Section 2.3. Through coordinate 
transformation Eq.(2.22) terms were separated in their power in volume. Thus, 
as system size is increased, influence of higher order normalized cumulants A n 
decreases, allowing for truncation of the summation for sufficiently large volume. 

A few words on physical units are in order. Discussion of the canonical en- 
semble shall suffice. The single particle partition function Eq.(2.9) ^[fm -3 ], and 



therefore all cumulant elements Eq.(2.1{ 



" [fm 3 ] in the CE. Conse- 



quently entries in Eq.(2.23) are er- fe [fm -3 / 2 ]. The normalization in Eq.(2.27), 
V J I 2 det |er|, for J-dimensional a, as well as the new variable of integration 
Eq.(2.22) $i, are hence dimensionless. The elements of the inverse sigma ten- 
sor are (cr -1 )^ [fm 3 / 2 ] and, thus, the elements £ J of the vector Eq.(2.26) will 
be dimensionless. Finally, the elements of the normalized cumulants, Eq.(2.25), 
are A^ lj2v " Jn [fm~ 3+3n / 2 ], which is canceled by the factor y- n / 2 + l i n the summa- 
tion in Eq.(2.27). Thus all terms involved in Eq.(2.27) are dimensionless. 

As discussed, for V — > oo one can discard terms of V~ x l 2 and higher in 
Eq.(2.27), and consider only the first two cumulants for the asymptotic solution: 



Z(V,(3, Ull , N ) j dH 
V L ' 2 det\a\ J (27r y 



cxp 



(2.28) 



This is the characteristic function of a multivariate normal distribution [119]. 
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Completing the square the integral (2.28) can be solved: 



Z Q (V, f3, Up, fij) 



(2nV) L/2 det |cr| 



exp 



(2.29) 



The asymptotic 
is then given by: 

P 9 ce(Q l ) 



solution of the GCE joint distribution of extensive quantities Q l 



Z^\V,^,u„ H ) 
Z(V,p,Up,u.j) 



(2%V) L/2 det 1^1 



exp 



2 « 



(2.30) 



Mean values in the thermodynamic limit are given by the first Taylor expansion 
terms, (N A ) = Vk* a , (Q) = Vnf, (E) = Vnf, etc. and converge to GCE values. 
To obtain a joint (two-dimensional) particle multiplicity distribution one has to 
take a two-dimensional slice of the (L dimensional) GCE distribution, Eq.(2.30), 
around the peak of the extensive quantities which one is considering as exactly 
fixed. Please note that Eq.(2.29) (albeit in different notation) was used clS cLIl clS - 
sumption in the microscopic correlator approach [99, 100, 101, 102]. More details 
of the calculation, in particular on the connection between the partition func- 
tions Z pii ^ (V, (3, Up, Hj) and the conventional version Z(V, P^, Q j ) [113, 114], 
can be found in Appendix A. Joint particle multiplicity distributions are con- 
structed from Eq.(2.30) in Appendix D. Details of the cumulant tensor are dis- 
cussed in Appendix C. A closed formula for the scaled variance of CE or MCE 
particle multiplicity fluctuations can be found in [105]. 



2.3 Finite Volume Corrections 

In considering finite system size effects on distributions, the region where 
the thermodynamic limit approximation is valid is left. Chemical potentials fij 
and four-temperature /3 M do not correspond anymore to the physical ones, which 
would be found in the GCE, but have to be thought of as Lagrange multipliers, 
used to maximize the partition function for a given (micro) canonical state. First 
some volume dependent correction terms are derived. A condition is obtained 
that defines the values of Uj and The correct choice allows to write down 
the thermodynamical potentials, the Helmholtz free energy F for the CE, and 
the entropy S for the MCE, in terms of the generalized partition function. Some 
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general criterion for the validity of the expansion is given. Approximations will 
be compared to analytical CE and MCE solutions of multiplicity distributions. 



2.3.1 Gram-Charlier Expansion 

In Section 2.2 it was shown that in the thermodynamic limit any equilib- 
rium distribution can be approximated by a multivariate normal distribution, 
Eq.(2.29). Further parameters, describing the shape of the distribution, like skew- 
ness (%), or excess/kurtosis (/C4), tend to zero as volume is increased. Returning 
to Eq.(2.27) with a number L of extensive quantities: 



Z Q (V,f3,u^fij] 



V L / 2 det\a\ J (2vr) J 



exp 



°° \h,h,...,ln 

z — ' n\ 

n=3 



. (2.31) 



Finite volume corrections will be obtained by Gram-Charlier expansion [120, 
121, 122, 123]. The charge vector Q 1 denotes a vector of L Abelian charges, and 
could read for a hadron resonance gas Q l = (B, S,Q, E, P x , P y , P z ), but could, 
in principle, also include particles multiplicities. Expanding the exponential in 
terms of powers in volume, one finds: 



z Q o/,/?, Wj ; 



Z(V,^,u^^) f dH 



V L ' 2 detkl 



(2 



exp 



7T 



3! 



VV2 



4! 



V 



+ 



1 \li,h,h \h, hM -6„q „q 



2! 3! 



r 



(2.32) 



Correction terms in Eq. (2.32) can be obtained by differentiation of exp [— 
with respect to the elements One can thus reverse the order by first integrating 
and then again differentiating. Using generalized Hermite polynomials: 



1 \ n 

-1) exp 



d" 



2 J dt h dZ h ...dZ ln 



exp 



(2.33) 
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with the adjusted shorthand notation for the contractions: 



MO 
K (0 



^3 

3! 
4! 



# 3 (0 



h,h,h 



(2.34) 
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(2.36) 



the partition function for finite volume can be approximated by: 



Z Q (V, (3, Up, fij) 
1 + 



Z(V,f3,Up,fij] 



x 



{2nV) L/2 det \a 
h 3 (0 h 4 (0 ^ 5 (0 
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^1/2 
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V 
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^3/2 
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o (v- 2 ) 



(2.37) 



Considering the simplest case of only one conserved charge, it is evident from 
Eq.(2.33), that the first order correction term in Eq.(2.37) is a polynomial of or- 
der 3 in £, while the second order correction term is a polynomial of order 4, etc. 
Hence for large values of e.g. a charge, energy, momentum, and multiplicity 
state far from the peak of the distribution will lead to a bad approximation, and 
even to negative values for P(Q l ). The validity of this approximation is thus 
restricted to the central region of the distribution. CE and MCE results will be 
compared to scenarios which are accessible to analytical methods in Section 2.3.3. 
In order to distinguish approximations which include corrections up to different 
orders in volume in Eq.(2.37), the asymptotic solution is denoted as CLT (cen- 
tral limit theorem), including terms up to O (V" 1 / 2 ) as GC3 (Gram-Charlier 3), 
including terms up to O (V -1 ) as GC4, and including terms up to O (V -3 / 2 ) 
as GC5. 
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2.3.2 Chemical and Thermal Equilibrium 



Here the question is addressed of how to choose the optimal values for (3^ 
and fj,j. The postulate is that a (micro) canonical equilibrium state should be 
as well the most likely state in the GCE. On the other hand, it is apparent that 
the expansion works best around the peak of the distribution. Hence /3 M and /Xj 
are chosen such that the partition function is maximized at some equilibrium 
point Q l eq . Taking terms up to 0{V~ 1 ^ 2 ) into account, the first derivative of the 
partition function Eq.(2.20) reads: 



dQ} 



1 



(2tt) l/2 T/(£+i)/2 det |cr| 



exp 



"5 ?<■ 



(2.38) 



3!^ V ' 1 L ' \kxMMM 



+ o (v- 1 ) 



Lagrange multipliers (or chemical potentials) should be chosen such that the first 
derivative Eq.(2.38) of Z Ql with respect to the conserved quantities Q l vanishes, 
hence Eq.(2.37) is maximized at the point Q l eq : 



dQ} 



(2.39) 



Using only the asymptotic solution, valid in the thermodynamic limit, this con- 
dition leads to: 

-l\k 



Zi = (Q k -VK l!k ) (tr" 1 )* 



0/ 



(2.40) 



Hence the partition function is maximal at the point Qk, e q — V K i,k- Charge 
and energy density correspond thus to the GCE values, and each component 
of ^ — > , and (3^ — > (3£ ce . While, when taking the first finite volume correction 
term in Eq.(2.38) into account, one obtains: 



A 



fel,fe2,fe3 



3\VV 







(2.41) 



rather than Eq.(2.40), and jJ ^ /i^ ce , and /3 M ^ /3 



gee 



ki,k 2 ,k 3 ,k4 

For calculation of distribu- 



tions for systems of finite volume one should therefore find chemical potentials 
that satisfy condition (2.39). A technical comment is in order. From Eq.(2.38) 
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it is evident that the first order correction term to the derivative of the partition 
function is a polynomial of order 4 in £, while the second one is of order 5, etc. 
It is therefore crucial to find in numerical calculations the correct maximum. 

2.3.3 Quality of Approximation 

To test the quality of the approximation for multiplicity distributions at finite 
volume for (very) small systems, the analytical solution for a classical CE particle- 
anti-particle gas [94], and a classical MCE (without momentum conservation) 
ultra-relativistic gas [1, 95] are compared to Eq.(2.37). 

Fig. (2.1) shows the multiplicity distribution of positively charged particles 
for two system sizes in the exact form and in different orders of approxima- 
tion Eq.(2.37) {top), and the ratio of approximation to exact solution (bottom). 
In Fig. (2.2) the same physical system is shown for a (relatively large) positive 
electric net-charge. Due to a one-to-one correspondence between the distribu- 
tions of negatively (suppressed) and positively (enhanced) particles one finds 
the multiplicity distribution P(N + ) generally more narrow than in the case of a 
neutral system. In particular towards the edge of the body of the distribution 
the approximation becomes worse. For the MCE massless gas approximations to 
multiplicity distributions (top) and ratios to the exact solution (bottom) are again 
compared in Fig. (2.3) for two system sizes. 

A few general comments attempt to summarize. The first observation is that 
indeed as system size is increased a better description of the central region is found 
in terms of the asymptotic solution CLT. The second observation is that even for 
systems with a very small number (in the order of 5) of produced particles one 
finds a good approximation in terms of Gram-Charlier expansion. In particular 
GC5 provides a very accurate description of the central region with deviations in 
the order of a few percent. This is quite remarkable given the fact that multiplicity 
distributions for such small systems are not smooth and continuous functions of 
multiplicity, while the approximation Eq.(2.37) is. Furthermore, implicitly the 
concepts of chemical potential and temperature are introduced for systems with 
small particle number, which may be in contradiction to the common believe that 
these parameters can only be meaningful when the number of involved particles 
becomes very large, i.e. in the thermodynamic limit. The last observation is 
that, indeed, see bottom rows of Figs. (2. 1-2.3), finite volume corrections (given 
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Figure 2.1: (Top): The CE ir + multiplicity distribution in Boltzmann approx- 
imation for temperature T = 160 MeV, electric net-charge Q — 0, and vol- 
ume V = 100 fm 3 (left), or volume V = 500 fm 3 (right). Exact solutions (solid), 
and various approximations, CLT (dash), GC3 (dash-dot), GC4 (dash- dot- dot), 
and GC5 (dot) are shown. (Bottom): same as (top), but the ratios of exact 
solution to approximations are shown. 

in terms of polynomials) lead only to good results for the central region of the 
distribution. 

To give an estimate for a region in which the approximation is reliable, one 
notes that the finite volume approximation scheme begins to break down when 
the first expansion term in Eq. (2.37) becomes unity. In the one-dimensional case 
this would be: 

h (U ax )_ _ Q (1) _ (2 42) 
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Figure 2.2: Same as Fig. (2.1), but for an electric net-charge of Q = 5 (left), or 
an electric net-charge of Q = 25 (right). 



Approximating the Hermite polynomial H 3 ~ £ , one can get an estimate for ^ max : 

/ o| \ 1/3 

Ua X — ( T~ J ^ 1/6 • ( 2 - 43 ) 



A 3 / 



While, when switching back to the definition of £ = w ^j^ 1 , Eq.(2.26), the width 
of the central region can be estimated by: 

iQ-WciU^ _ f3\ \ 1/3 y 2 /3 _ / 2 44 x 



Hence the width of the central region scales as V 2 ^ 3 , while the width of the dis- 
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Figure 2.3: (Top): MCE multiplicity distribution for an ultra-relativistic gas in 
Boltzmann approximation. For energy E — 1 GeV, and volume V = 50 fm 3 (left), 
or for E = 3 GeV, and V = 150 fm 3 (right). Exact solutions (solid), and various 
approximations, CLT (dash), GC3 (dash-dot), GC4 (dash- dot- dot), and GC5 
(dot) are shown. (Bottom): same as (top), but the ratios of exact solution to 
approximations are shown. 

tribution should scale as V 1 ^ 2 and the approximation should be quite good. Even 
though larger volumes work better, they will still be sufficiently small enough 
to allow for calculation of distributions relevant for heavy ion collisions. How- 
ever, it is stressed that there is no simple criterion for what is a 'small' or 
a 'large' volume for a particular physical system. Formally the existence and 
fmiteness of (at least) the first three cumulants k is sufficient for application of 
the asymptotic solution [120, 121, 122, 123]. Considering the simple case of a 
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multiplicity distribution of Bose-Einstein particles in the GCE one finds from 
Eq.(2.9) < k 2 ' N < Kz' N,N < ■ ■ ■ ■ Hence, in particular when finite chemical 
potentials are involved, cumulants are growing with order, implying that apart 
from mean value and variance further parameters like skewness and excess (or 
kurtosis) [119, 124] of the distribution will remain important quantities. For 
chemical potentials approaching the Bose-Einstein condensation point, one finds 
for instance the cumulant K^ ,Ar diverging. The expansion discussed here is hence 
not valid in this regime. 



2.4 Temperature and Chemical Potential 

The introduction of chemical potentials in the CE and temperature in 
the MCE was first and foremost a mathematical trick which allowed to con- 
veniently integrate partition functions for which otherwise no analytical solution 
could be obtained. However the generalized partition function is self-consistent 
and not in contradiction to the common definition of temperature and chemical 
potential. It is shown in the following that the definition of and \x through 
Eq.(2.39) coincides with expressions well known from textbooks [110, 111, 112]. 

Canonical Ensemble 

The canonical partition function known from textbooks and the generalized 
version discussed here are connected as follows: 

Z(V,Q\(3) = Z Qj (V : /3,u.j) e~ Qj ^ . (2.45) 

The Helmholtz free energy F is the thermodynamic potential relevant for CE, 

F ee -T\nZ CE . (2.46) 

Employing the first law of thermodynamics dE = TdS — PdV + [ijdQ\ and 
using F = E — TS for the free energy, where P is the pressure, and E and S are 
total energy and entropy, respectively, one can write for the differential of the 
free energy dF = —SdT — PdV + fijdQi . The Lagrange multiplier fj,j associated 
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with conserved charge Qi is thus defined by: 

JSlI 6 = „. (2 47) 



where advantage of condition (2.39) was taken. Thus, the correct choice for the 
effective chemical potential is given by the extremum, = 0, of the charge dis- 
tribution, which coincides with fj,j = (j^frj > Eq.( 2.47). In the thermodynamic 
limit this is equivalent to, Eq.(2.40), (Aj —> fij, gC e- The subscripts indicate that the 
derivative with respect to the conserved charge has to be taken at fixed values 
of V and (3. In the non-relativistic case, where particle number N, rather than 
quantum numbers, is conserved, the corresponding relation to Eq.(2.47) would 

be im) Y ji = ^ t 110 ' m > 112 1- 

Micro Canonical Ensemble 

The standard MCE partition function can be obtain by Fourier integration of 
the generating function of the GCE four-momentum distribution and multiplica- 
tion with the inverse Boltzmann factor, 

Z (V, P^) = Z p " (V, p, u„) e p - u ^ . (2.48) 

The relevant thermodynamic potential in the MCE is the entropy S, 

S = hxZ^P^) . (2.49) 

The Lagrange multiplier for the temperature is [110, 111, 112, 113, 125, 126]: 

wX = — = v " p ' (2 - 50) 

where again condition (2.39) was employed, ^§p— = 0. The subscript in Eq.(2.50) 
is used to indicate that the derivative with respect to the components of P^ has 
to be taken at fixed volume V. Thus Eq.(2.50) resembles the optimal choice of 
an effective temperature for the approximation scheme. In the thermodynamic 
limit, V -> oo, one finds f}» -» /3£ ce , due to Eq.(2.40). 
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Grand Canonical Ensemble 

Conventionally, e.g. in textbooks, first the MCE is introduced. Summation 
over energy - with temperature being a Lagrange multiplier - which is used to 
maximize the entropy, introduces the CE. Additionally dropping the constraint 
of exact charge conservation leads to the GCE. Here the chemical potentials fij 
are the Lagrange multipliers. The GCE partition function can be written as: 

Z(V,/3,u M ,^) = Yl ^ mP^ Qi) (2.51) 
{<2 3 } 

= eQ ^ jP e ~ P " U ^ z (V,P",Q j ) ■ (2.52) 

{Qi} {PM} 

While, in the notation used here, this line would read: 

Z{V,P,u^m) = ^Z&Wfru^iMj) (2.53) 
{Q j } 

= EE Z QJ ^(V,f3,u„ H ) . (2.54) 

{Qi} {P"} 

The thermodynamic potential for the GCE is the grand potential Q: 

Q = -T\nZ(V,f3,u^fij) = -T\nZ(V, /?, 9{) . (2.55) 

The relations (2.47) and (2.50) are essentially Legendre transformations from F 
and S to Z (V,Qi , (3) and Z(V,P^), respectively. The inverse Legendre trans- 
formations were obtained in the saddle-point expansion method (see Refs. [113, 
114, 115]). Lastly it is emphasized that both methods are complementary from 
a thermodynamic point of view. 

Thus, if exact solutions of the canonical or micro canonical partition func- 
tions were available this reversal would not have been necessary. However, this 
redefinition of the GCE partition function is entirely consistent and simplifies cal- 
culations considerably. Whenever an exact solution to the generalized partition 
function is possible, all relations above would hold exactly. This interpretation 
of the GCE partition function as the generating (or characteristic) function of a 
statistical system should be quite useful, even in more general cases than the one 
presented here. 
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2.5 Discussion 

An analytical expansion method for calculation of distributions at finite vol- 
ume for the canonical as well as the micro canonical ensembles of the ideal rela- 
tivists hadron resonance gas has been presented. The introduction of tempera- 
ture into the micro canonical partition function and chemical potentials into the 
canonical partition function have lead to the identification of the grand canonical 
partition function with the characteristic function of associated joint probabil- 
ity distributions. Micro canonical and canonical multiplicity distributions could, 
thus, be defined through conditional probability distributions. Hence the prob- 
ability of finding a certain multiplicity, while other parameters (global charge or 
energy) were taken to be fixed. 

In considering finite volume corrections to the system partition function, thus 
relaxing the assumption of thermodynamic equivalence of different statistical en- 
sembles, one is lead to demand that the partition function should be maximized 
for a particular set of conserved charges. It turned out that this requirement is 
entirely equivalent to the well known textbook definitions of chemical potential 
in the canonical ensemble as the derivative of Helmholtz free energy with respect 
to conserved charge and temperature in the micro canonical ensemble through 
differentiation of entropy with respect to conserved energy. 

This method is based on Fourier analysis of the grand canonical partition 
function. Conventionally one would not introduce chemical potentials and tem- 
perature into these calculations. However, one then faces the problem of a heavily 
oscillating (or even irregular) integrand, making numerical integration unpracti- 
cal. Artificially introduced temperature and chemical potentials, correctly chosen, 
produce a very smooth integrand allowing for expansion of the integrand in pow- 
ers of volume. Analytical solutions to asymptotic distributions could thus be 
found in terms of Laplace's expansion, while finite volume corrections could be 
obtained from Gram-Charlier expansion. A comparison with available analytical 
solutions to simple statistical systems suggests that good results can be expected 
even for rather small system size. One drawback is that the results can only be 
applied to the central region of the distribution, owning to the fact that finite 
volume correction terms appear in the form of Hermite polynomial of low order. 

Another, rather practical, drawback is that the effects of resonance de- 
cay [101, 103, 105] and limited acceptance [106, 127] are hard to consider simul- 
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taneously with the approach presented here. This will be solved by the Monte 
Carlo approach, considered in the next chapter. 
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Chapter 3 

Monte Carlo Approach 

A statistical hadronization model Monte Carlo event generator provides the 
means for studying fluctuation and correlation observables in equilibrium systems. 
Data analysis can be done in close relation to experimental analysis techniques. 
Imposing global constraints on a sample is always technically more challenging. 
Direct sampling of MCE events (or micro states) has only been done in the non- 
relativistic limit [128, 129]. Sample and reject procedures, suitable for relativistic 
systems, become rapidly inefficient with increasing system size. However, they 
have the advantage of being very successful for small system sizes [113, 114, 115]. 

In this chapter a different approach is taken: the GCE is sampled, events 
are then re-weighted according to their values of extensive quantities, and the 
sample-reject limit (MCE) is approached in a controlled manner. In this way 
one can study the statistical properties of a global equilibrium system in their 
dependence on the size of its thermodynamic bath. As any of the three stan- 
dard ensembles remain idealizations of physical systems, one might find these 
intermediate ensembles to be of phenomenological or conceptual interest too. 

In Section 3.1 statistical ensembles with finite thermodynamic bath are dis- 
cussed. The GCE Monte Carlo sampling procedure is described in Section 3.2. 

3.1 Statistical Ensembles with Finite Bath 

The starting point taken is similar to the one chosen by Patriha [92], and 
Challa and Hetherington [97]. However quickly a different route is taken. Two 
systems, described by their micro canonical partition functions, i.e. the number 
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of micro states for two separate systems, are assumed. The first system is further 
assumed to be enclosed in a volume V\ and to have fixed values of extensive quan- 
tities Pf = (Ex, P Xt i, P y> i, P z ,i), and Q{ = (E>i, Si, Qi), while the second system 
is assumed to be enclosed in a volume V 2 and to have fixed values of extensive 
quantities Pf = (E 2 , P x , 2 , P y ,2, Pzfi), and Qi = (B 2 ,S 2 ,Q 2 ), where E is the en- 
ergy of the system, P x , y , z are the components of its three-momentum, and B, S, 
and Q, are baryon number, strangeness and electric charge, respectively. Thus: 

Z(Vi,P?,Q{) = Z Ni (Vi,P?,Q{) , and Z(V 2 , P£, Qi) , (3.1) 

{N{} 

where Z N i(Vi, Pf , Q{) denotes the number of micro states of system 1 with addi- 
tionally fixed multiplicities N[ of particles of all species i considered in the model. 
Supposing that system 1 and system 2 are subject to the following constraints: 

V g = V x + V 2 , (3.2) 
P? = Pi + P? , (3-3) 

Qi = Q{ + Qi • (3.4) 

Then, the partition function Z(V g , P£, Q J g ) of the joint system is constructed as 
the sums over all possible charge and energy-momentum split-ups: 

z(v g ,p»,Qi) =J2Y, z ( v 9 - v ^ P 9- p ^Ql-QD Z(Vi,P?,Q{)- (3-5) 

{Q{} i p D 

The distribution of extensive quantities in the subsystem Vi is given by the ratio 
of the number of all micro states consistent with a given charge and energy- 
momentum split-up and a given set of particle multiplicities to the number of all 
possible configurations: 

P(P?,Q{,N{) = Z{V ° - V ^f " Q{) Z N{ (Vi,P?,Q{) ■ (3-6) 

Lastly, a weight factor W(Vi, Pf, Q{; V g , P^, Q J „) is defined such that: 

P(P?,Q{,N[) = W(Vi,P?,Q{;V g ,P»,Qi) Z N i(Vi,P?,Q{) . (3.7) 
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By construction, the first moment of the weight factor is equal to unity: 

w = EE w{ y^ p ^ ^ Q j 9 ) ZNi(y u pf, q{) = i , (3.8) 

{Q{}{ p ?} 

as the distribution is properly normalized. 

The weight factor W(Vi, Pf , Q{; V g , P£, Qi) generates an ensemble with sta- 
tistical properties different from the limiting cases of vanishing, V g — > V\ (MCE), 
and of an infinite thermodynamic bath, V g — > oo (GCE). This effectively allows 
for extrapolation of GCE results to the MCE limit. In principle any other (ar- 
bitrary) choice of W(Vi, Pf , Q{; V g , P{f, Q 3 g ) could be taken. This thesis confines 
itself to the situation discussed above. It is worth noting that all micro states 
consistent with the same set of extensive quantities (Pf, Q±) still have 'a priori 
equal probabilities'. 

In the large volume limit, ensembles are equivalent in the sense that densities 
are the same. The ensembles defined by Eq.(3.7) and later on by Eq.(3.11) are no 
exception. If both V\ and V g are sufficiently large, then the average densities in 
both systems will be the same, Q g /V g and P g /V g respectively. The system in V\ 
will hence carry on average a certain fraction: 

A = Vx/Vg , (3.9) 

of the total charge Q J g and four-momentum Pjf, i.e.: 

(Q{) = A Q j g , and (Pf) = A P» . (3.10) 

By varying the ratio A = Vi/V g , while keeping (Q{) and (Pf) constant, one can 
thus study a class of systems with the same average charge content and four- 
momentum, but different statistical properties. In the thermodynamic limit (it 
is enough to demand that V\ is sufficiently large) a family of thermodynamically 
equivalent (same densities) ensembles is generated. 
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3.1.1 Introducing the Monte Carlo Weight W 

Since Eq.(3.7) poses a formidable challenge, both mathematically and numer- 
ically, one may write instead: 

P{P^Q{,N{) = W r:x ^^(\]:V :l W/',) P gce (P?,Q{,Ni\P,u„ H ) , 

(3.11) 

where the distribution of extensive quantities Pf, Q{ and particle multiplicities N{ 
of a GCE system with temperature T = volume Vi, chemical potentials fij 
and collective four-velocity w M is given by: 

where /ij = (/ie, /is, fig) summarizes the chemical potentials associated with 
baryon number, strangeness and electric charge in a vector. The normalization in 
Eq.(3.12) is given by the GCE partition function Z(Vi, f3, u M , /ij), i.e. the number 
of all micro states averaged over the Boltzmann weights e~ PlU ^ and e^ 1 ^ 13 : 

ZiV^fi^tij) = E E E e ~ Pl ^ eQ{NP Z Ni {V u P?,Q{) . (3.13) 
The new weight factor }/\? p i >Qi'> p a '^(Vi; V g \/3, m m , /ij) now reads: 



z(v g ,Pt;,Qi) 

In the case of an ideal (non-interacting) gas, Eq.(3.14) can be written, see Chap- 
ter 2, and Appendix A, as: 

x Z^-Pi^^l (Vg -V U 0, M^, 

Z p ^(V g ,^,^) 

The advantage of Eq.(3.11), compared to Eq.(3.7), is that the distribution 
P gce (Pi, Q{, N[\(3, u^, /ij) can easily be sampled for Boltzmann particles, while 
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a suitable approximation for the weight )/\? p i'Qi> p s '^(Vx; Vg\/3, u^, [ij) is avail- 
able. Again, by construction, the first moment of the new weight factor is equal 
to unity: 

< w > = E E E vv p -^ p *^ = i. 

(3.16) 

In principle, Eq.(3.7) and Eq.(3.11) are equivalent. In fact, Eq.(3.7) can be 
obtained by taking the limit (//#, /is, ^q) = (0,0,0), = (1,0,0,0), and /? — > 
of Eq.(3.11). However, as one can easily see, (W n ) ^ (W n ). Higher, and in 
particular the second, moments of the weight factors W and W are a measure 
of the statistical error to be expected for a finite sample of events. The larger 
the higher moments of the weight factor, the larger the statistical error, and the 
slower the convergence with sample size. Please see also Appendices F and G. 

As GCE and MCE densities are the same in the system V g , these values are 
effectively regulated by intensive parameters (3, fij and w M . In order to study a 
system with average (Q{), it is most economical to sample the GCE with (Q{) 
and calculate weights according to Eq.(3.15). This will result in a low statistical 
error for finite samples (as shown in later sections), and allow for extrapolation 
to the MCE limit. 

Firstly the weight factor Eq.(3.15) will be calculated, and then the relevant 
limits are taken. With the appropriate choice of 0, fij and u M the calculation of 
Eq.(3.15) can be done with the method presented in Chapter 2. 

3.1.2 Calculating the Monte Carlo Weight W 

For the Monte Carlo approach in this chapter, the total number of (poten- 
tially) conserved extensive quantities in a hadron resonance gas is L = J + 4 = 
3 + 4 = 7, where J = 3 is the number of charges (B, S, Q) and there are four 
components of the four-momentum. Particle number integrals are not going to 
be solved here. Including all extensive quantities into a single vector: 



Q l = (Q j , P») = (B, S, Q, E, P x , P y , P z ) , 



(3.17) 
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the weight Eq.(3.15) can be expressed as: 



W^HVr^u,,^ = Z(V^u^ 3 ) x Z& '~^1 Yhh^il .(3.18) 

The general expression for the partition function Z Ql (V, (3,u^, fij) in the large 
volume limit is given by Eq.(2.29): 



Z Q (V^u^/ij) 



exp 



(2ttV) l / 2 det \a 
where, unlike in Eq.(2.26), the volume is split off: 

while the definition of the sigma tensor, Eq.(2.23), remains: 



1 1 f l f 



(3.19) 



(3.20) 



4 



-1/2 



(3.21) 



Here K\ and are the GCE vector of mean values and the GCE covariance 
matrix respectively. The values of (3, \i 3 - and are chosen such (Eq.(2.39)) that: 



dz Ql 



dQ} 



2 



0, 



(3.22) 



The approximation (3.19) gives then a reliable description of Z Ql 3 around the 
equilibrium value Q l g = V 9 k,[, provided V g is sufficiently large. The charge vec- 
tor, Eq.(2.26), is then equal to the null- vector £j — 0/ (Q l g — V g n\). For the 
denominator in Eq.(3.18) one then finds: 



while for the numerator one obtains: 



4=Q' g , eq (2nV g ) L / 2 det\a\ 



exp [ 



(3.23) 



_ Z(Vg - Vl.^U^jlj) 

~ (27^-^))^ d et | a | 
1 1 



x exp 



2 {Vg ~ V 



(3.24) 
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where the charge vector Eq.(3.20), in Eq.(3.24) is: 

e = (AQ 2 ) k . 
Then, using the condition Qg = Qg t£q = V 3 k\ yields: 



(3.25) 



(AQ 2 



[Qg-Ql) -{Vg-V^Kl = -{Ql-VxKt] 



(3.26) 



Substituting Eq.(3.23) and Eq.(3.24) into Eq.(3.18) results in: 



=Ql 



X 



Z(Vg,f3,U^fij) 

(2TrV g ) L / 2 det |a| 
{2n{V g -V l )) L l 2 det|a| 

1 1 



x exp 



2 (V g - V,) 



(3.27) 



The GCE partition functions are multiplicative in the sense that their product 
is Z(Vi, (3, Up, fij) Z(V g — Vx,/3, Up, jij) = Z(V g , (3, u^, fij), and thus the first term 
in Eq.(3.27) is equal to unity. Now using Eq.(3.9), A = V\/V g , one can re-write 
Eq.(3.27) as: 



W^HV^VglP,^,^) 



o l =o l 

^3 ~^9,e<l 



(1 - A) L / 2 



x exp 



(3.28) 



\ (i=a) vk ^ 



Model parameters are hence the intensive variables inverse temperature /3, four- 
velocity vF 1 and chemical potentials p? , which regulate energy and charge den- 
sities, and collective motion. Provided V\ is sufficiently large, a family of ther- 
modynamically equivalent ensembles is defined, which can now be studied in 
their dependence of fluctuation and correlation observables on the size of the 
bath Vi = V g — V\. Hence, one can test the sensitivity of such observables, for 
example, to globally applied conservation laws. The expectation values are then 
identical to GCE expectation values, while higher moments will depend crucially 
on the choice of A. 
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3.1.3 The Limits of W 

The largest weight is given to states for which = 0, i.e. with extensive 
quantities Q[ = Qi eq ,- Hence, the maximal weight a micro state (or event) at a 
given value of A = V\/V g can assume is Wmax 9 {Vi, V g \/3, fij) = (1 — \)~ L ' 2 . 
Taking the limits of Eq.(3.28), it is easy to see that: 

hm W Ql ^ Q HVi;V g \/3, UlMN ) = 1. (3.29) 

A — ^0 

I.e. for A = the GCE is sampled, and all events have a weight equal to unity. 
Hence, one also finds (W 2 ) = 1 and therefore ((AW) 2 ) = 0, implying a low 
statistical error. For A — » 1, effectively the sample-reject limit [113, 114, 115] is 
approached. Accordingly: 

hm W^iVuVgl^u^fij) oc 5(Q[-V 1 k 1 1 ) . (3.30) 

A — y 1 

However, as now not all events have equal weight, ((AW) 2 ) grows and so too 
the statistical error of finite samples. Also, the larger the number L of extensive 
quantities considered for re-weighting, the larger will be the statistical uncer- 
tainty. 

3.2 The GCE Sampling Procedure 

The Monte Carlo sampling procedure for a GCE system in the Boltzmann 
approximation is now explained. The system to be sampled is assumed to be 
in an equilibrium state enclosed in a volume V\ with temperature T = and 
chemical potentials fij = (fJ>B, A*s> A*q)- Additionally, the system is assumed to 
be at rest. The four-velocity is then w M = (1, 0, 0, 0) and the four-temperature 
is /3 M = (/?, 0, 0, 0). In this case, multiplicity distributions are Poissonian, while 
momentum spectra are of Boltzmann type. The GCE sampling process is com- 
posed of four steps, each discussed below. 

Multiplicity Generation 

In the first step, multiplicities iV{ are randomly sampled of all particle species % 
considered in the model. The expectation value of the multiplicity of thermal 
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Boltzmann particles in the GCE is given by: 

(Nl) = |4 m * Z 3 " 1 K * . (3.31) 



Multiplicities {Nl} n are randomly generated for each event n according to Pois- 
sonians with mean values (Nl): 

P(N{) = ^ e-M> . (3.32) 

In the above, and gi are the mass and spin-degeneracy factor of a particle of 
species i respectively. The chemical potential /ij = fijqj = ush + UsSi + p>Q<li, 
where q\ = (fr^Sj,^) represents the quantum number content of a particle of 
species i. 

Momentum Spectra 

In the second step, momenta are generated for each particle according to 
a Boltzmann spectrum. For a static thermal source spherical coordinates are 
convenient: 

g = ZT- W ^ ,3,3) 
These momenta are then isotropically distributed in momentum space. Hence: 

Px = bl sin^ cos</> , (3.34) 
p y = \p\ sin 9 sin (ft , (3.35) 
Pz = bl COS6 1 , (3.36) 



e = VH 2 + m ?> ( 3 - 37 ) 

where p x , p y , and p z are the components of the three-momentum, e is the energy, 
and \p\ = y/p* + Py+ p 2 z is the total momentum. The polar and azimuthal angles 
are sampled according to: 

6 = cos' 1 [2 (x - 0.5)] , (3.38) 
= 2 7r(x-0.5) , (3.39) 
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where x is uniformly distributed between and 1. Additionally, the transverse 
momentum p T and rapidity y are calculated for each particle: 

Pt = yfti+P* , (3-40) 

Finally, particles are distributed homogeneously in a sphere of radius r\ and decay 
times are calculated based on the Breit-Wigner width of the resonances. 



Resonance Decay 

The third step (if applicable) is resonance decay. Following the prescription 
used by the authors of the THERMINATOR package [130], only 2 and 3 body 
decays are performed, while unstable daughter particles are allowed to decay in a 
successive manner. Only strong and electromagnetic decays are considered, while 
weak decays are omitted. Particle decay is first calculated in the parent's rest 
frame, with daughter momenta then boosted into the lab frame. Finally, decay 
positions are generated based on the parent's production point, momentum and 
life time. 

Throughout this thesis, always only the lightest states of the following 
baryons: 

p n A £+ TT E~ E° Q~ (3.42) 
and their respective anti-baryons, as well as following mesons: 

7T+ tT tt° K + K- K° K° (3.43) 

are considered as stable. The system could now be given collective velocity u M . 

Re- weighting 

In the fourth step, the values of extensive quantities are calculated for the 
events generated by iteration over the particle list of each event. For the val- 
ues of extensive quantities Q[ n = {B 1>n , S 1)Tl , Q 1;U , E 1>n , P x>1 , n , P y ,i, n , P z ,i,n) in 
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subsystem V\ of event n one may write: 

= E <4 > ( 3 - 44 ) 

particles i n 

where q- n = (b in , Si n ,qi n ,ei n ,p x ,i n ,Py tin ,p z , in ) is the 'charge vector' of particle % in 
event n. Based on Q\ n the weight w n is calculated for the event: 

Wn = W^Hv^Vglfru^u.j) , (3.45) 

according to Eq.(3.28). All micro states (or events) with the same set of extensive 
quantities Q l l n are still counted equally. 

Monte Carlo Distributions 

The Monte Carlo output is essentially a distribution Pmc{X\, X2, X3, ...) of a 
set of observables X±, X 2 , X 3 , etc. For all practical purposes this distribution is 
obtained by histograming all events n according to their values of X\ >n , X2, n , Xs tn , 
etc. and their weight w n . One can then define moments of two observables Xi 
and Xj through: 

(X?X?) = X ? X T PMc{X u Xi) . (3.46) 

Xi,Xj 

Additionally, the variance {(AXiY) and the covariance (AXiAXj) respectively 
are, defined as: 

((AX,) 2 ) = (X 2 ) - (X t ) 2 , and (3.47) 
(AXiAXj) = (XX,) - (X)(X J ) . (3.48) 

In the following, the correlation coefficient pij, defined as: 



\/((AA- i ) 2 ){fAA- J ) 2 > 



3.3 Discussion 
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and later on the scaled variance Ui, given by: 

are used for quantification of fluctuations and correlations. The scaled variance Ui, 
is often also called the 'Fano Factor' [107], and is a useful measure for the width 
of a distribution, provided Xj > 0, such as in the case of particle multiplicity. 



3.3 Discussion 

In this chapter distributions of extensive quantities in equilibrium systems 
with finite thermodynamic bath have been discussed. A recipe for a thermal 
model Monte Carlo event generator has been presented. The event generator is 
capable of extrapolating fluctuation and correlation observables for Boltzmann 
systems of large volume from their GCE values to the MCE limit. This approach 
has a strong advantage compared to analytical approaches or standard micro 
canonical sample-and-reject Monte Carlo techniques, in that resonance decays as 
well as (very) large system sizes can handled at the same time. 

To introduce the scheme, a micro canonical system has been conceptually 
divided into two subsystems. These subsystems are assumed to be in equilibrium 
with each other, and subject to the constraints of joint energy-momentum and 
charge conservation. Particles are only measured in one subsystem, while the 
second subsystem provides a thermodynamic bath. By keeping the size of the 
first subsystem fixed, while varying the size of the second, one can thus study the 
dependence of statistical properties of an ensemble on the fraction of the system 
observed (i.e. assess their sensitivity to globally applied conservation laws). The 
ensembles generated are thermodynamically equivalent in the sense that mean 
values in the observed subsystem remain unchanged when the size of the bath is 
varied, provided the combined system is sufficiently large. 

The Monte Carlo process can be divided into four steps. In the first two steps 
primordial particle multiplicities for each species, and momenta for each particle, 
are generated for each event by sampling the grand canonical partition function. 
In the third step resonance decay of unstable particles is performed. Lastly, the 
values of extensive quantities are calculated for each event and a corresponding 
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weight factor is assigned. All events with the same set of extensive quantities 
hence still have 'a priori equal probabilities'. In the limit of an infinite bath, 
all events have a weight equal to unity. In the opposite limit of a vanishing 
bath, only events with an exactly specified set of extensive quantities have non- 
vanishing weight. In between, extrapolation is done in a controlled manner. The 
method is particularly efficient for systems of large volume, inaccessible to sample- 
and-reject procedures, and agrees well, where available, with analytic asymptotic 
micro canonical solutions. 

Given the success of the hadron resonance gas model in describing experi- 
mentally measured average hadron yields, and its ability to reproduce low tem- 
perature lattice susceptibilities, the question arises as to whether fluctuation and 
correlation observables also follow its main line. The first and, in particular, 
second moments of joint distributions of extensive quantities will be studied in 
the following. The focus will be mainly on particle multiplicity distributions and 
distributions of 'conserved' charges. In particular, the effects of resonance decay, 
conservation laws, and limited acceptance in momentum space are discussed. Due 
to the Monte Carlo nature, data can be analyzed in close relation to experimental 
analysis techniques. The hadron resonance gas is an ideal testbed for this type 
of study, in that it is simple and intuitive. 
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Chapter 4 

Grand Canonical Ensemble 

The grand canonical ensemble of the hadron resonance gas model will be the 
starting point of the discussion of statistical properties of equilibrium systems. 
The GCE is considered to be the most accessible amongst the standard ensembles. 
Its discussion will provide the basis for the study of canonical and micro canonical 
ensembles in the following chapters. The GCE is the correct ensemble to choose, 
if the observed subsystem is embedded into a much larger heat and charge bath. 
Primordial momenta of particles within a small subvolume are then uncorrelated 
with each other in the Boltzmann approximation. 

Extensive quantities are fluctuating, while corresponding intensive variables 
are constant. Hence, in a subvolume one finds the energy content (extensive) 
to fluctuate, while temperature (intensive) is constant. Chemical potentials are 
constants, but the corresponding net-charge content is fluctuating. And naturally 
in a relativistic gas, the systems total particle number and integrated transverse 
momentum (no Lagrange multiplier associated with them) also fluctuate. How- 
ever, extensive quantities do not merely fluctuate, but are also correlated with 
each other. 

On the other hand, correlations observables do not depend on where they are 
measured in coordinate space, due to the infinite bath assumption. Observable 
correlation and fluctuation signals are, however, sensitive on the location and size 
of the available acceptance window in momentum space within which particles 
can be measured. 

In Section 4.1 joint distributions of fully phase space integrated extensive 
quantities are considered. In Section 4.2 momentum spectra of primordial and 
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final state hadrons are analyzed. Correlations between extensive quantities are 
studied in their dependence on the acceptance window in momentum space in 
Sections 4.3 and 4.4. 

4.1 Joint Distributions of Extensive Quantities 

A neutral and static GCE hadron resonance gas with local temperature 
T = 0.160 GeV and volume V = 2000 fm 3 has been sampled according to the 
procedure described in Section 3.2. The values of net-charges and energy and 
momentum Q\ n are calculated by iteration for the particle list of each event, 
according to Eq.(3.44). In the GCE all events have weight w n — 1, Eq.(3.45). 
In Figs. (4. 1-4.3) the event output is summarized. GCE joint distributions of 
extensive quantities are shown for samples of 5 ■ 10 5 events. 

Charge fluctuations directly probe the degrees of freedom of a system, i.e. 
they are sensitive to its particle mass spectrum (and its quantum number config- 
urations). The correlation coefficient and variances carry information about ori- 
entation and elongation of the probability distributions in Figs. (4.1-4.3). Firstly, 
charge correlations and the contributions of different particle species to the co- 
variance (AXiAXj) , Eq.(3.48), and therefore to the correlation coefficient p,jj, 
Eq.(3.49), are considered. 

All baryons have baryon number b = +1. Baryons can only carry strange 
(valence) quarks, i.e. their strangeness is always s < 0. Anti-baryons have b = 
— 1, and s > 0. I.e., both groups contribute negatively to the baryon-strangeness 
covariance, and so (ABAS) < 0, and therefore pss < 0, as seen in Fig. (4.1) (left) 
and indicated by the solid lines in Fig. (4.5). 

Positively charged baryons and their anti-particles contribute positively to 
the baryon-electric charge covariance (ABAQ), while negatively charged baryons 
(and their anti-particles) contribute negatively. Two observations can be made 
on the hadron resonance gas mass spectrum: there are more positively charged 
baryons than negatively charged ones, and their average mass is lower. I.e., in a 
neutral gas (//g = Pq = — 0) the contribution of positively charged baryons 
(and negatively charged anti-baryons) dominates and therefore (ABAQ) > 
and pbq > 0, as seen in Fig.(4.1) (right) and indicated by the solid lines in 
Fig.(4.6). 
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B B 

Figure 4.1: Grand canonical joint distribution of baryon number B and 
strangeness S (left), and grand canonical joint distribution of baryon number B 
and electric charge Q (right), for a neutral hadron resonance gas at T = 0.160 GeV 
and V = 2000 fm 3 . Here 5 • 10 5 events have been sampled. 



x10' 3 xlO" 3 




S N 

Figure 4.2: Grand canonical joint distribution of strangeness S and electric 
charge Q (left), and grand canonical joint distribution of energy E and total 
particle number N (right). The rest as in Fig. (4.1). 

Mesons and their anti-particles always contribute positively to the strangeness- 
electric charge correlation coefficient psq. Electrically charged strange mesons are 
either composed of an w-quark and a s-quark, or of an w-quark and a s-quark 
(and superpositions thereof). Their contribution to (ASAQ) is in either case 
positive. On the baryonic side, only the S + (as well as its degenerate states and 
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N N 

Figure 4.3: Grand canonical joint distribution P(N, Pt) of total particle num- 
ber N and total transverse momentum Pt {left), and grand canonical joint dis- 
tribution P(N,pf) of total particle number N and mean transverse momentum 
pr = Pt/N (right). The rest as in Fig. (4.1). 

their respective anti-particles) has a negative contribution to (ASAQ) , while all 
other strangeness carrying baryons have either electric charge q = —1, or q = 0. 
Therefore, one finds psq > 0, as seen in Fig. (4.2) (left) and indicated by the solid 
lines in Fig. (4. 7). 

Since there are more particles (at \xb = 0, and T = 0.160 GeV) which carry 
electric charge than particles which carry strangeness, one finds the electric charge 
variance to be larger than the strangeness variance ((AQY) > ((AS) Z ). Similarly, 
as there are more strangeness or anti-strangeness carrying particles than thermal 
baryon and anti-baryons, strangeness fluctuations are found to be stronger than 
baryon number fluctuation, ((AS) 2 ) > ((AB) 2 ). 

The correlation between particle number N and energy E in Fig. (4.2) (right) 
is more obvious. The more particles are in a subvolume, the more energy it will 
contain at a certain temperature P^ 1 . The total transverse momentum Pt can 
be seen as a measure for the kinetic energy in the system, Fig. (4.3) (left). And 
naturally, the more particles are inside the subvolume, the more kinetic energy is 
stored in them. These extensive quantities are strongly correlated in the GCE. 
Due to the independent sampling of primordial particle momenta in the GCE, 
however, the mean energy per particle e = E/N and the mean transverse momen- 
tum pr = Pt/N, are event-by-event un-correlated with the particle number N of 
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the events. Here, the mean transverse momentum versus particle number distri- 
bution is shown in Fig. (4.3) (right). 

4.2 Momentum Spectra 

To prepare a more detailed analysis of charge fluctuations and correlations, 
inclusive momentum spectra are analyzed. In Fig. (4.4) transverse momentum and 
rapidity spectra of positively charged hadrons, both primordial and final state, 
are shown for a static thermal system. 

Based on these momentum spectra acceptance bins ApT,i and Ay, are con- 
structed following [98, 106, 127] and [59, 131]. Momentum bins are chosen such 
that each of the five bins contains on average one fifth of the total yield of posi- 
tively charged particles. The values defining the bounds of the momentum space 
bins Apx t i and Ay^ are summarized in Table 4.1. 




Pt [GeV] ■ y 

Figure 4.4: (Left:) Transverse momentum spectrum of positively charged 
hadrons, both primordial and final state. (Right:) Rapidity spectrum of posi- 
tively charged hadrons, both primordial and final state. Here 2 ■ 10 6 events have 
been sampled for a static neutral hadron resonance gas in Boltzmann approxi- 
mation with temperature T = 0.160 GeV. 

Resonance decay shifts the transverse momentum distribution to lower aver- 
age transverse momentum (pt) and widens the rapidity distribution of thermal 
'fireballs' [132, 133]. Transverse momentum bins of final state spectra are hence 
slightly 'contracted', while rapidity bins get slightly 'wider', when compared to 
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their respective primordial counterparts. 





Pt.i [GeV] 


Pt,2 [GeV] 


PT,3 [GeV] 


Vta [GeV] 


PT,5 [GeV] 


PT.e [GeV] 


Primordial 


0.0 


0.22795 


0.36475 


0.51825 


0.73995 


5.0 


Final state 


0.0 


0.17105 


0.27215 


0.38785 


0.56245 


5.0 




2/i 


2/2 


2/3 


2/4 


2/5 


2/6 


Primordial 


-5.0 


-0.4275 


-0.1241 


0.1241 


0.4273 


5.0 


Final state 


-5.0 


-0.5289 


-0.1553 


0.1551 


0.5289 


5.0 



Table 4.1: Transverse momentum and rapidity bins Ap<r,i = [pt,i,Pt,i+i] and 
A|/j = [yi, jji+i], both primordial and final state, for a static neutral hadron reso- 
nance gas in Boltzmann approximation with temperature T = 0.160 GeV. 



Resonance decay combined with transverse as well as longitudinal flow pro- 
vides a rather good description of experimentally observed momentum spectra 
in relativistic heavy ion collisions at SPS and RHIC energies [130, 134, 135]. 
The momentum spectra discussed here, on the other hand, contain no flow, and 
results, thus, cannot be directly compared to experimental data or transport sim- 
ulations. However, qualitatively one might observe the effects discussed in the 
following sections. 

4.3 Correlations between Charges 

An interesting example of quantities - for which the measured value depends 
on the observed part of the momentum spectrum - are the correlation coefficients 
between the charges baryon number B, strangeness S and electric charge Q. Also 
variances and covariances of the baryon number, strangeness, and electric charge 
distribution are strongly sensitive to the acceptance cuts applied. Their values 
are additionally rather sensitive to the effects of globally enforced conservation 
laws. If the size of the 'bath' is reduced, a change in one interval of phase space 
will have to be balanced (preferably) by a change in another interval, and not by 
the (finite) 'bath'. 

The correlation coefficients pbs, Pbq, and Psq in limited acceptance bins Ap T ti 
and At/i, as defined in Table 4.1, are considered in the grand canonical ensem- 
ble. Primordial particles in one momentum bin are then sampled independently 
from particles in any other momentum space segment, due to the 'infinite bath' 
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assumption. Nevertheless, the way in which quantum numbers are correlated is 
different in different momentum bins 1 , as different particle species have, due to 
their different masses, different momentum spectra. Resulting correlation func- 
tions are hence spectra, rather than simply values. 

Figs. (4.5-4.7) show the correlation coefficients between baryon number and 
strangeness, pss, between baryon number and electric charge, Pbq, and between 
strangeness and electric charge, psq, as measured in the acceptance bins Ap T)i 
and Ai/i defined in Table 4.1, both primordial and final state. The average baryon 
number, strangeness, and electric charge in each bin is equal to zero, since the 
system was assumed to be neutral. The primordial values (15 bins) shown in 
Figs. (4.5-4.7) are calculated using analytical solutions. The error bars on many 
data points are smaller than the symbol used. 
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Figure 4.5: Baryon-strangeness correlation coefficient pbs in the GCE in limited 
acceptance windows, both primordial and final state in transverse momentum 
bins Apx,i {left), and in rapidity bins Ayi [right). Horizontal error bars indicate 
the width and position of the momentum bins (And not an uncertainty!). Vertical 
error bars indicate the statistical uncertainty of 20 Monte Carlo runs of 10 5 events 
each. The marker indicates the center of gravity of the corresponding bin. The 
solid lines show the fully phase space integrated GCE result. 



Tables 4.2 to 4.4 summarize the transverse momentum and rapidity depen- 
dence of the correlation coefficients pss, Pbq, and psq. The statistical error 
quoted corresponds to 20 Monte Carlo runs of 10 5 events each. The analytical 



-"Tor the extensive quantities energy E and longitudinal momentum P z the correlation in 
limited acceptance bins ApT,i and Ayi is more apparent than for quantum numbers. 
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Figure 4.6: Baryon-electric charge correlation coefficient pbq in the GCE in lim- 
ited acceptance windows, both primordial and final state in transverse momentum 
bins ApT,i {left), and in rapidity bins Ayi (right). The rest as in Fig. (4.5). 
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Figure 4.7: Strangeness-electric charge correlation coefficient psq in the GCE in 
limited acceptance windows, both primordial and final state in transverse momen- 
tum bins ApTj (left), and in rapidity bins Ay^ (right). The rest as in Fig. (4.5). 



values (5 bins) listed in the tables are calculated using the momentum bins de- 
fined in Table 4.1. Mild differences between Monte Carlo and analytical results 
are unavoidable. The calculated values are also not exactly symmetric in Ayi, 
as the exact size of the acceptance bins constructed is sensitive to the number of 
bins used for the calculation of the momentum spectra. Correlation coefficients p 
are also rather sensitive to exact bin size, and the fourth digit becomes somewhat 
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PBS 


A 

Apr.l 


A 

Ap T . 2 


Apr. 3 


Ap T a 


A 

Apr,5 


n calc 
Pprim 

Pprim 

P final 


-0.2479 
-0.248 ±0.003 
-0.216 ±0.002 


-0.2641 
-0.264 ±0.003 
-0.220 ±0.003 


-0.2864 
-0.286 ±0.003 
-0.241 ± 0.004 


-0.3188 
-0.319 ±0.002 
-0.269 ±0.003 


-0.3839 
-0.385 ±0.002 
-0.335 ±0.003 


PBS 






Ay 3 


Ay 4 


Ay 5 


Pprim 
Pprim 
P final 


-0.2407 
-0.241 ±0.003 
-0.191 ±0.002 


-0.3345 
-0.334 ±0.003 
-0.300 ±0.002 


-0.3536 
-0.353 ±0.003 
-0.328 ±0.002 


-0.3345 
-0.335 ±0.003 
-0.299 ±0.002 


-0.2408 
-0.240 ±0.003 
-0.190 ±0.002 



Table 4.2: Baryon-strangeness correlation coefficient pss i n the GCE in transverse 
momentum bins Apx,i and rapidity bins Ayi, both primordial and final state. For 
comparison, analytical values p^im f° r primordial correlations are included. The 
statistical uncertainty corresponds to 20 Monte Carlo runs of 10 5 events each. 



PBQ ApT,l Apr,2 Apr,3 ApTA ApT,5 



„calc 
Pprim 

Pprim 

P final 


0.1120 
0.113 ±0.002 
0.112 ±0.003 


0.1271 
0.126 ±0.002 
0.120 ±0.003 


0.1420 
0.143 ±0.003 
0.138 ±0.003 


0.1579 
0.158 ±0.002 
0.164 ±0.003 


0.1781 
0.178 ±0.003 
0.221 ±0.003 


PBQ 


Ayi 


Aj/2 


Ay 3 


Ay 4 


Ay 5 


n calc 
Pprim 

Pprim 

P final 


0.1160 
0.116 ±0.002 
0.118 ±0.003 


0.1601 
0.160 ±0.003 
0.192 ±0.003 


0.1658 
0.166 ±0.003 
0.202 ±0.003 


0.1601 
0.159 ±0.003 
0.192 ±0.003 


0.1160 
0.117 ±0.002 
0.119 ±0.003 



Table 4.3: Baryon-electric charge correlation coefficient pbq in the GCE in trans- 
verse momentum bins Ap^ j and rapidity bins A?/j, both primordial and final 
state. 



Psq Apta Ap T:2 Ap T ,3 Ap TA Ap T £ 



n calc 
Pprim 

Pprim 

P final 


0.2831 
0.284 ±0.003 
0.243 ±0.003 


0.3033 
0.304 ±0.003 
0.254 ±0.003 


0.3150 
0.314 ±0.003 
0.276 ±0.003 


0.3185 
0.319 ±0.002 
0.292 ± 0.003 


0.3055 
0.305 ±0.002 
0.303 ±0.002 


PSQ 


Ayi 


Ay 2 


Ay 3 


Ay 4 


Ay 5 


n calc 
Pprim 

Pprim 

P final 


0.2934 
0.294 ±0.003 
0.255 ±0.002 


0.3137 
0.314 ±0.003 
0.299 ±0.003 


0.3104 
0.310 ±0.002 
0.297 ±0.003 


0.3137 
0.312 ±0.003 
0.298 ±0.003 


0.2934 
0.292 ±0.002 
0.255 ±0.003 



Table 4.4: Strangeness-electric charge correlation coefficient psq in the GCE in 
transverse momentum bins ApT,i and rapidity bins Ayi, both primordial and final 
state. 
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unreliable. 

The rapidity dependence of pss, Pbq, and psq is now considered. Strange 
baryons are, on average, heavier than non-strange baryons, so their thermal rapid- 
ity distributions are narrower. The kaon rapidity distribution is then, compared 
to baryons, again wider. A change in baryon number (strangeness) at high \y\ is 
less likely to be accompanied by a change in strangeness (baryon number) than 
at low \y\. The value of pbs, therefore, drops toward higher rapidity, as shown in 
Fig. (4.5) (right). 

By the same argument, one finds a weakening of the baryon-electric charge 
correlation pbq at higher rapidity, Fig. (4.6) (right)), as the rapidity distribution 
of electrically charged particles is wider than the one of baryons. The strangeness- 
electric charge correlation coefficient psq first rises mildly and then drops some- 
what stronger towards higher rapidity. As one shifts ones acceptance window to 
higher values of \y\, first the contribution of baryons (in particular S + ) decreases 
and, as the meson contribution grows, psq rises slightly. Towards the highest \y\, 
pions again dominate and de-correlate the quantum numbers. 

The transverse momentum dependence can be understood as follows: heav- 
ier particles have higher average transverse momentum (pt) and, hence, their 
influence increases towards higher p?. Heavy particles, in particular baryonic 
resonances, often carry several quantum numbers, causing the correlation coeffi- 
cients to grow. 

The contribution of strange baryons compared to non-strange baryons grows 
for higher transverse momentum, since strange baryons have on average larger 
mass than non-strange baryons. The correlation coefficient pbs, thus, becomes 
strongly negative at high px- As the contribution of baryons compared to mesons 
grows stronger towards larger p?, a change in baryon number (electric charge) is 
now more likely to be accompanied by a change in electric charge (baryon num- 
ber) than at low pr, and pbq increases with pr (the A resonances 2 ensure it keeps 
rising). For the Ap^ j dependence of psq one notes that one of the strongest con- 
tributors at higher p? is the Q~ , with a relatively low mass of mn- = 1.672 GeV. 
So after a rise, psq drops again towards highest pt, due to an increasing S + 
contribution 3 . 

Since resonance decay has the habit of dropping the lighter particles (mesons) 

included in the THERMITS particle table up to the A(2420) . 
included in the THERMUS particle table up to the £(2030). 
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at low pt and higher \y\, while keeping heavier particles (baryons) at higher pj< 
and at mid-rapidity, none of the above arguments about the transverse momen- 
tum and rapidity dependence are essentially changed by resonance decay. The 
correlation coefficient pss becomes more negative towards higher p?, while be- 
coming weaker towards higher \y\. Similarly, pbq grows larger at high px and 
drops with increasing y. The larger contributions of baryons to the high p? tail 
of the transverse momentum spectrum, and their decreased contribution to the 
tails of the rapidity distribution, compared to mesons, are to blame. The bump 
in the pt dependence of psq, caused by the S + , has vanished, as the S + is only 
considered as stable in its lightest version with mass mg+ = 1.189 GeV. The 
small bump in the y dependence of psq, however, stays. The correlation is first 
increased by a growing kaon contribution and then again decreased by a growing 
pion contribution at larger rapidities. 

4.4 Correlation Functions 

Above results are now applied to measures of correlations used in heavy ion 
phenomenology. The measurement of charge correlations appears to be a good 
discriminating tool between different physics scenarios, such as the ones provided 
by lattice QCD [136], effective field theories [137], or non-equilibrium transport 
models [138]. Acceptance effects are studied in this section. A discussion of 
the dependence of fully phase space integrated quantities on temperature and 
chemical potential can be found in Chapters 8 and 9. In this context variances 
of marginal distributions are studied. 

Following the definitions of Ref. [139] one can relate the correlation coefficients 
Pbs an d Psq to the baryon number strangeness correlation measure Cbs'- 



In the ideally weakly coupled QGP the fully phase space integrated values are for 





and the electric charge strangeness correlation measure Cqs- 




(4.2) 
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the baryon number strangeness correlation measure C^ P = 1, and for the electric 
charge strangeness correlation measure C® GP = 1 [139, 140]. On the other hand, 
if correlations were determined by hadronic degrees of freedom, the fully phase 
space integrated correlation measures would be C^g G ~ 0.66, and Cq^ g ~ 1.2 re- 
spectively, for a hadron resonance gas at T = 0.170 GeV and \ib = 0.0 GeV [139]. 
The 'quark Molecular Dynamics' (qMD) model [138] is in agreement with hadron 
resonance gas results below the critical temperature T c , while following lattice 
QCD calculations above T c . Additionally, the baryon number electric charge 
correlation measure Cbq is defined to be: 



Since the motivation in [139] for having the baryon number variance {(AB) 2 ) 
not in the denominator, was that the electrically neutral neutron is hard to mea- 
sure, the covariance in Eq.(4.3) is divided by the electric charge variance ((AQ) 2 ), 
which is experimentally more accessible. The fully phase space integrated hadron 
resonance gas estimate, Fig. (4.9), is C^q G ~ 0.23. A QGP estimate is not re- 
ported. 

In Figs. (4.8-4. 10) the correlation measures Cbs (baryon number - strangeness), 
Cbq (baryon number - electric charge), and Csq (strangeness - electric charge), 
are shown as measured in the acceptance bins Ap^j and Ayi defined in Table 4.1, 
both primordial and final state. The average baryon number, strangeness, and 
electric charge in each bin is equal to zero, since the system the assumed to be 
neutral. The analytical primordial values (15 bins) shown in Figs. (4.8-4.10) are 
calculated using analytical methods. It is mentioned that, again, the error bars 
on many data points are smaller than the symbol used. Fully phase space inte- 
grated values, indicated by solid lines, are compatible of hadron resonance gas 
results of Ref. [139]. Nevertheless the results here show a strong dependence on 
the acceptance cuts applied. The fully phase space integrated values give a good 
estimate of the average. The data in limited acceptance are consistently above the 
weakly coupled QGP estimates [139, 140] in the case of Cbs, or below as for Cqs- 
The momentum space dependence of charge correlations in the QGP has not yet 
been calculated. In the hadron resonance gas the momentum space dependence 
is caused by different momentum spectra due to different hadron masses. 
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Figure 4.8: Baryon-strangeness correlation measure Cbs m the GCE in limited 
acceptance windows, both primordial and final state, in transverse momentum 
bins ApT,i {left), and in rapidity bins Ayi (right). Horizontal error bars indicate 
the width and position of the momentum bins (And not an uncertainty!). Vertical 
error bars indicate the statistical uncertainty of 20 Monte Carlo runs of 10 5 events 
each. The marker indicates the center of gravity of the corresponding bin. The 
solid lines show the fully phase space integrated GCE result. 
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Figure 4.9: Baryon-electric charge correlation measure Cbq in the GCE in limited 
acceptance windows, both primordial and final state, in transverse momentum 
bins ApT,i (left), and in rapidity bins Ayi (right). The rest as in Fig.(4.8). 



One notes a qualitatively similar behavior of the correlation coefficients p 
shown in Figs. (4.5-4.7) to the correlation measure C. To illustrate their dif- 
ferences, joint distributions are considered in limited acceptance. The correla- 
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Figure 4.10: Strangeness-electric charge correlation measure Csq in the GCE 
in limited acceptance windows, both primordial and final state, in transverse 
momentum bins Apx,i (left), and in rapidity bins Ayi (right). The rest as in 
Fig. (4.8). 



tion measures p and C are similar in their Ap^ j and Ayi dependence, and are 
essentially explained by arguments above. Additionally, the momentum space 
dependence of variances in Eqs.(4.1-4.2) needs to be taken into account. 

Fig. (4.11) shows the GCE joint baryon number strangeness distribution as cal- 
culated in the acceptance windows Apx^ (top) and Ayi (bottom), while Fig. (4. 12) 
summarizes the primordial baryon number variance ((AB) ), stran eeness vari- 
ance ((AS) 2 ), and the baryon number strangeness covariance (ABAS) in the 
GCE in transverse momentum bins Apx,i (left) and rapidity bins Ayi (right). 

The marginal variances ((AB) 2 ) and ((AS*) 2 ) get wider in the transverse 
momentum direction Ap Ti , Fig. (4. 11) (top), while the correlation pss gets 
stronger. The contribution of baryons (compared to all hadrons) increases to- 
wards higher transverse momentum, which causes the baryon distribution to 
widen. Strangeness carrying particles are heavier than non-strange particles, so 
their variance increases, too. Together they drive up their covariance ( AB AS) . In 
rapidity Ay^ direction, Fig. (4. 11) (bottom), the arguments are somewhat reversed. 
Particles carrying baryon number and strangeness are predominantly produced 
(thermally) around mid-rapidity. At larger rapidity, the pions take over, and the 
baryon number and strangeness variances decrease. And as the kaon rapidity dis- 
tribution is wider than the one of the A, one finds the correlation pss to decrease 




Figure 4.11: (Top:) Joint baryon number strangeness distribution in transverse momentum bins, from left to right Apr 
to ApT,5- (Bottom:) Joint baryon number strangeness distribution in rapidity bins, from left to right Ay± to Ay§. 
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towards higher \y\. 
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Figure 4.12: Primordial baryon number variance ((AB) 2 ), stran eeness variance 
((AS) 2 ), and the baryon number strangeness covariance (ABAS) in the GCE in 
limited acceptance windows: (Left:) transverse momentum bins Apx,i- (Right:) 
rapidity bins Ayi. Horizontal error bars indicate the width and position of the 
momentum bins (And not an uncertainty!). Vertical error bars are smaller than 
the statistical uncertainty of 20 Monte Carlo runs of 10 5 events each. The marker 
indicates the center of gravity of the corresponding bin. 

Similarly, in full acceptance one finds for the variances of the marginal charge 
distributions in limited acceptance ((AB) 2 ) < ((AS) 2 ) < ((AQ) 2 ). Yet, the de- 
pendence of their ratios on transverse momentum and rapidity cuts is stronger 
for some and weaker for others. Baryons dominate at high pt; strangeness car- 
rying particles are heavier; accordingly both variances increase at large pt com- 
pared to the electric charge variance (not shown). Heavier particles live around 
mid-rapidity. The ratios ((AB) 2 ) / ((AQ) 2 ) and ((AB) 2 ) / ((AS) 2 ) get larger at 
larger and \y\, while the ratio ((AQ) 2 ) / ((AS) 2 ) gets smaller at larger momen- 
tum. Therefore, the momentum dependence is stronger for Cbs and Cbq than 
for p BS and pbq, while being weaker for Cqs than for p S g. 



4.5 Discussion 



In this chapter samples of grand canonical Monte Carlo events have been an- 
alyzed. Joint distributions of extensive quantities have been considered with full 
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acceptance, or with limited acceptance, in momentum space assumed. Attention 
was given to the effects of resonance decay. 

The correlations and fluctuations of extensive quantities of a statistical system 
are determined by the degrees of freedom, i.e. particles and their masses and 
quantum numbers, available to the system. So is, for instance, the correlation 
between energy and particle multiplicity dependent on the type of particle species 
measured and the experimental acceptance window assumed. 

Due to the assumption of an infinite thermodynamic bath, occupation num- 
bers of individual momentum levels are un-correlated with each other. Likewise, 
particle multiplicities of two distinct groups of particles appear un-correlated. To- 
gether with occupation number fluctuations unconstrained by global constraints, 
all extensive quantities, bar the volume, fluctuate on an event-by-event or micro 
state-by-micro state basis. Therefore, the energy content and particle multiplic- 
ity are strongly correlated, while the average energy per particle is un-correlated 
with multiplicity. The net-charge content of baryon number and strangeness in 
the system are correlated because some hadron species carry both charges. Dif- 
ferent hadron species have different quantum numbers and different spectra. The 
correlation of baryon number and strangeness depends, like the one of energy and 
momentum, or energy and particle multiplicity, on which part of the spectrum is 
accessible for measurement. 

The correlation coefficients between quantum numbers increase towards 
higher pr, since heavier particles have higher average transverse momentum (pr)- 
Heavy particles, in particular baryonic resonances, often carry several quantum 
numbers, causing the correlation coefficients to grow. For a thermal fireball this 
effect is strongest in a mid-rapidity window. 

The values of p or C after resonance decay are directly sensitive to how the 
data is analyzed. In the above study only final state particles (stable against 
strong decays) are analyzed. One could, however, also reconstruct decay positions 
and momenta of parent resonances and could then count them as belonging to 
the acceptance bin the parent momentum would fall into. In the situation above, 
however, this would again yield the primordial scenario. If reconstruction of 
resonances is not done, one is sensitive to charge correlations carried by final 
state particles. As in the primordial larger acceptance bin effectively 

averages over smaller bins. However, the smaller the acceptance bin, the more 
information is lost due to resonance decay. In full acceptance, final state and 
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primordial correlation coefficients ought to be the same, since quantum numbers 
(and energy-momentum) are conserved in the decays of resonances (weak decays 
omitted). 
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Chapter 5 

Extrapolating to the MCE 

In the previous chapter the grand canonical ensemble and limited acceptance 
effects have been studied. In this chapter, fully phase space integrated exten- 
sive quantities are considered and studied in their dependence on the size of the 
thermodynamic bath for a neutral and static system. Firstly, mean values, co- 
variances, and correlations coefficients of joint distributions of extensive quantities 
are analyzed in Section 5.1. Then, in Section 5.2, the joint distributions them- 
selves, and, in Section 5.3, the weight factor are considered more closely. This 
discussion will prepare the analysis of multiplicity fluctuations and correlations, 
Chapter 6, and the effect of conservation laws on them. 

5.1 Fully Phase Space Integrated Quantities 

Fully phase space integrated grand canonical results are extrapolated to the 
micro canonical limit. For this purpose hadron resonance gas events for various 
values of A = V\/V g are iteratively generated, re- weighted, and analyzed. By 
construction of the weight factor W, Eq.(3.28), the extrapolation proceeds in a 
systematic fashion such that, for instance, particle momentum spectra as well as 
mean values of extensive quantities remain unchanged. On the other hand, all 
variances and covariances of extensive quantities subject to re-weighting converge 
linearly to their micro canonical values. 

This can be seen from the form of the analytical large volume approximation 
to the grand canonical distribution of (fully phase space integrated) extensive 
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quantities P gce (Q[) (from Eq.(2.30)): 

where the variable £ z is given by Eq.(3.20). Now taking the weight factor W\, 
Eq.(3.28), (<t and ^ are the same in both equations) one obtains for the distri- 
bution P\(Q[) of extensive quantities Q± in subsystem 1: 

Px(Q[) * wf i;S * P gce (Q[) 

l 

~ (2vr(l - A)Vx) L / 2 detH 6XP 

This is essentially the same multivariate normal distribution as the grand canon- 
ical version P gC e{Q[), however linearly contracted. Monte Carlo results will then 
be compared to Eq.(5.2). 

Again a static and neutral hadron resonance gas with collective four- 
velocity = (1,0,0,0), chemical potential vector /ij = (0,0,0), local temper- 
ature T = = 0.160 GeV, and volume V\ = 2000 fm 3 is considered. This is 
a system large enough 1 for using the large volume approximation worked out in 
Section 2.2. 

In Figs. (5.1) and (5.2) the results of Monte Carlo runs of 2.5 • 10 4 events each 
are shown. Each value of A has been sampled 20 times to allow for calculation 
of a statistical uncertainty estimate. 19 different values of A have been studied. 
In this case study, the extensive quantities baryon number B, strangeness S, 
electric charge Q, energy E, and longitudinal momentum P z are considered for 
re- weighting. Conservation of transverse momenta P x and P y can be shown not to 
affect the Apx,i and Ai/i dependence of multiplicity fluctuations and correlations 
studied in Chapter 6. Their Ayi dependence is, however, rather sensitive to P z 
conservation. Angular correlations, considered in Chapter 7, on the other hand, 
are strongly sensitive to joint P x and P y conservation [106, 127]. 

Fig. (5.1) (left) summarizes the results for mean values of baryon number (B), 
strangeness (S), electric charge (Q), energy (E), and the momenta (P x ) and (P z ). 

1 Generally it is not easy to define when a system is 'large enough' for the large volume 
approximation to be valid. Here, good agreement is found with asymptotic analytic solutions. 
Charged systems, or Bose-Einstein/Fermi-Dirac systems, usually converge more slowly to their 
asymptotic solution. 
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Figure 5.1: Mean values (left) and variances (right) of various extensive quanti- 
ties, as listed in the legends, as a function of A. Each marker and its error bar 
represents the result of 20 Monte Carlo runs of 2.5 ■ 10 4 events each. 19 different 
equally spaced values of A have been investigated. Solid lines indicate GCE val- 
ues (left), or linear extrapolations from the GCE value to the MCE limit (right). 
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Figure 5.2: Covariances (left) and correlation coefficients (right) between various 
extensive quantities, as listed in the legends, as a function of A. Solid lines 
indicate linear extrapolations from the GCE value to the MCE limit (left), or 
GCE values (right). The rest as in Fig.(5.1). 



The solid lines represent GCE values. Only the expectation value of energy is 
not equal to 0, since the system sampled is assumed to be static and neutral 
with T ^ 0. The evolution of the respective variances is shown in Fig. (5.1) (right). 
Variances of extensive quantities subject to re-weighting converge linearly to 
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as A goes to 1. One notes that ((AP X ) 2 ) remains constant (within error bars), 
as this quantity is not re-weighted in this case study. Error bars on many data 
points are smaller than the symbol used. 

In Fig. (5.2) (left) the evolution of covariances (ABAS), (ABAQ), (ASAQ), 
and (AEAQ) is shown as a function of the 'size of the bath'. As seen, the 
covariances between quantities considered for re-weighting converge linearly to 0. 
In a neutral system, covariances between energy and charges are equal to 0. 
As an example, (AEAQ) is shown. In a static system, also the covariances 
between momenta and any other extensive quantity are equal to 0. As an example 
here (AEAP Z ) is shown. The correlation coefficients, Eq.(3.49), on the other 
hand, remain constant as a function of A, as shown in Fig. (5.2) (right). The 
values of fully phase space integrated correlation coefficients pss, Pbq, and psq 
can be compared to the GCE results denoted by the solid lines shown in Figs. (4.5- 
4.7) in Section 4.3. 

5.2 Probability Distributions 

It is also interesting to study probability distributions as a function of the 
size of the bath A. Firstly joint distributions of charges, P\(B,S), P\(B,Q), 
and P\(S,Q) from (left) to (right), are attended to in Fig. (5. 3). Mean values of 
net-charges stay constant throughout the extrapolation, see also Fig. (5.1) (left). 
Their variances and covariances, compare Fig. (5.1) (right) and Fig. (5.2) (left), 
converge to as the size of the bath is reduced. All three quantities, B, S, and Q 
are re- weighted, and thus correlation coefficients amongst them, Fig. (5.2) (right), 
stay constant. Hence events in a small region around (Bi), (Si), and (Qi), are 
highlighted, while correlations between them are left unchanged 2 . 

In Fig. (5.4) (left) the joint energy and particle number distribution P\(E, N) 
is shown. Here only one of the two quantities was re-weighted. In the GCE one 
finds particle number to be strongly correlated with energy. The more particles 
are in a box of volume V\, the more energy is contained inside this box. This 
distribution does obviously not converge to a 5-function, as the size of the bath is 

2 Sometimes it may not be clear which ensemble one should choose to apply to a given 
statistical system. The three standard ensembles remain particular idealizations. Also the 
scenario discussed here is a particular limit. The system under investigation could have a bath 
which might impose its own correlations. Or in practical terms, have a more general form of 
the weight W^ 1 then the one discussed here. 
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BBS 

Figure 5.3: Evolution of the joint distribution of baryon number and 
strangeness, P\(B,S) (left), baryon number and electric charge, P\(B,Q) (cen- 
ter), and strangeness and electric charge, P\(S,Q) (right), with the size of the 
bath A = Vi/Vg. Here 5 • 10 5 events have been sampled for each value of A. 
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Figure 5.4: Evolution of the joint distributions of energy and particle number, 
P\(E,N) {left), particle number and transverse momentum, P\(N,Pt) (center), 
and particle number and mean transverse momentum, P\(N,p~r) (right), with the 
size of the bath A = V\/V g . Here 5 • 10 5 events have been sampled for each value 
of A. 
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reduced. In direction of energy the distribution gets successively narrower, while 
retaining some width along the particle number direction. Along the way, the 
correlation coefficient for finite size of the bath changes. Therefore, the correlation 
coefficient pen — > as A — > 1. 

Finally, distributions with neither quantity being re-weighted are considered. 
The distributions of particle number and total transverse momentum P\(N, Pr), 
Fig. (5.4) (center), and of particle number versus mean transverse momentum 
P\(N,pr~ = Pt/N), Fig.(5.4) (right), will serve as examples. Although the total 
transverse momentum Pt is strongly correlated with particle number N in the 
GCE, the mean transverse momentum pr = Pt/N is not. The more particles 
are inside a box, the more kinetic energy is contained within its boundaries. Due 
to the infinite heat bath assumption, individual particle momenta are however 
un-correlated with each other. Measuring a particle with a certain momentum 
vector, does not constrain the remaining system in any way or shape. As the size 
of the bath is now reduced, this does not hold true anymore. The total trans- 
verse momentum is still positively correlated with particle number, yet weaker 
as A — >■ 1. The mean transverse momentum distribution makes this rather 
plain. At fixed energy, the larger the particle multiplicity, the more the energy 
goes into mass, and the less energy into thermal motion, or vice versa. The 
correlation coefficient is hence negative pjvpr < 0, and the more particles are 
measured in an event, the less kinetic energy the individual particle will carry on 
average. Please note, that although the distribution P\(N,pr) changes, neither 
of the mean values does; (N) and (pr) stay constant. The expectation value of 
mean transverse momentum at fixed particle number (pt)n, obtained from the 
conditional distribution P(pr\N), however, decreases as N increases. 

5.3 Monte Carlo Weight Factor 

The Monte Carlo sample was therefore successively transformed. With de- 
creasing size of the bath, A — > 1, larger and larger weight is given to events in the 
immediate vicinity of the equilibrium expectation value, and smaller and smaller 
weight to events away from it. The distribution of extensive quantities consid- 
ered for re-weighting (a multivariate normal distribution in the GCE in the large 
volume limit) hence gets contracted to a ^-function with vanishing variances and 
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Figure 5.5: Evolution of the weight factor distribution with the size of the 
bath A = Vi/Vg. The occurrence of certain weight factors is shown with the 
value of the extensive quantities energy E (left), baryon number B (center), and 
particle number N (right). The extensive quantities B, S, Q, E, P z are consid- 
ered for re-weighting. Here 5 • 10 5 events have been sampled for each value of A. 
The solid vertical lines indicate the maximal weight w™ ax = (1 — A) L ^ 2 . 
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Figure 5.6: Evolution of the weight factor distribution with the size of the 
bath A = Vi/Vg. The occurrence of certain weight factors is shown with the 
value of the extensive quantities energy E (left), baryon number B (center), and 
particle number N (right). Only the extensive quantities B, S, Q are considered 
for re- weighting. Here 5 • 10 5 events have been sampled for each value of A. The 
solid vertical lines indicate the maximal weight w™ ax = (1 — A)~ L ^ 2 . 
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covariances. I.e., successively the properties of events which have very similar 
values of extensive quantities are highlighted. This will have a bearing on charge 
correlations and, in particular, multiplicity fluctuations and correlations which 
will be discussed in the following sections. 

Considering the evolution of the weight factor distribution with the size of 
the bath A = V\/V g , Fig. (5.5) shows the occurrence of certain weight factors w n , 
Eq.(3.45), with the value of the extensive quantities energy E\^ n (left), baryon 
number B l n (center), and particle number Ni >n (right). 

In the GCE, A = all events have 
weight equal to unity as discussed in 
Chapter 3. In the limit A — > 1, i.e. 
approaching the MCE, one finds for 
quantities which are re-weighted, that 
basically all events have weight equal 
to 0, except for those in the immedi- 
ate vicinity of the equilibrium value. 
Not all events with the same value of 
energy E^ n receive the same weight 
Wr? 1 '^ s . The energy value might be 
close to the peak. If the baryon num- 
ber (and strangeness, longitudinal mo- 
mentum, etc.) value are now also close 
to their equilibrium values, the events 
obtain a higher weight, if the baryon 
value is far away, a lower weight. The 
largest weight an event can assume is, 
Eq.(3.28), w^ ax = (1 — A) . For intermediate A a small hill, where most events 
are located, emerges. As A is increased most events receive smaller and smaller 
weight and this small hill gets pushed to the left. The evolution of the baryon 
number distribution, Fig. (5.5) (center) is similar. 

The extensive quantity particle number is not re-weighted. Hence its distri- 
bution does not converge to a 5-function. The connection between the extensive 
quantity N and other (re-weighted) quantities manifests itself indirectly. Events 
have certain energy, baryon number or strangeness. These quantities are corre- 
lated with particle number. Events with unusually large or small multiplicities 
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Figure 5.7: First and second moment of 
the weight factor Eq. (3.28) function 
of A. The rest as in Fig. (5.1). 
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now receive mostly small weight factors, because chance has it, that they also 
have values of energy, net-baryon number or net-strangeness which are far from 
the equilibrium expectation value. But even in the limit A — > 1 a certain width 
remains, and no sharp tip forms, as particle number continues to fluctuate in the 
MCE. In Fig. (5.7) the second moment of the weight factor, Eq.(3.28) is shown as 
a function of A for both the extrapolation to the MCE and to the CE. A large 
second moment (W 2 ) implies a large statistical uncertainty and, hence, usually 
requires a larger sample. 

For comparison, the evolution of weight factors for an extrapolation to the 
CE limit is presented in Fig. (5.6). Events are re- weighted only according to the 
values of their extensive quantities baryon number B\^ n , strangeness Si >n , and 
electric charge Qi, n . As in Fig. (5.5), the distributions for the extensive quantities 
energy E, baryon number B, and particle multiplicity iV are shown for differ- 
ent values of A. The energy distribution, Fig. (5.6) (left), is much wider than 
the one in Fig. (5.5) (left). The systems energy content continues to fluctuate, 
((AE)) 7^ 0, in the CE. The baryon number distribution Fig. (5.6) (center) on the 
other hand converges, as in the MCE, to a 5-function. The particle number distri- 
bution Fig. (5.6) (right) is then also wider than for an extrapolation to the MCE, 
Fig. (5.5) (right). The last observation is that, indeed, the maximal weight w™' 1 * 
an event can receive (at some fixed value of A) is smaller for the extrapolation to 
the CE, Fig. (5.6), than for an extrapolation to the MCE, Fig. (5.5). 

5.4 Discussion 

In this chapter fully phase space integrated grand canonical results have been 
extrapolated to the micro canonical limit. For this purpose samples of events 
for various values of A = V\/V g have been iteratively generated, re-weighted, and 
analyzed. By construction of the weight factor W the extrapolation proceeds in 
a systematic fashion such that, for instance, particle momentum spectra, as well 
as mean values of extensive quantities, remain unchanged. 

The GCE distribution of extensive quantities considered for re- weighting con- 
verges to a 5-function. Mean values as well as correlation coefficients stay con- 
stant, while variances and covariances converge to 0. Events in the vicinity of the 
equilibrium expectation value are projected out, the rest is suppressed. 
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Although the extensive quantities particle multiplicity N and total transverse 
momentum are not re-weighted, their fluctuations and correlations are never- 
theless strongly affected. The effects of decreasing size of the heat bath emerge 
because N and Pt are correlated with E, Q, etc., which are re-weighted. 

As the system under investigation gets larger, a smaller fraction of events 
will have values of Ei jn ,B l n , Qi >n , etc. in the vicinity of a desired equilibrium 
value, and the straight sample and reject method becomes ever more inefficient. 
This method can handle large system sizes, due to additional information ob- 
tained from the extrapolation. However, as A grows, so too does the statistical 
uncertainty. In the limit A — > 1, one approaches a sample-reject type of for- 
malism. One cannot, therefore, directly obtain the micro canonical limit for the 
large system size studied here, as this is prohibited by available computing power. 
However, extrapolation to this limit is possible. It is pointed out in this context 
that the intermediate ensembles, between the limits of GCE and MCE, may be 
of phenomenological interest, too. 
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Chapter 6 

Multiplicity Fluctuations and 
Correlations 

In this chapter joint multiplicity distributions of charged hadrons are analyzed 
in different ensembles with limited acceptance in momentum space assumed. Mul- 
tiplicity fluctuations and correlations are qualitatively affected by the choice of 
ensemble and are directly sensitive to the fraction of the system observed. For 
vanishing size of ones acceptance window, one would lose all information on how 
the multiplicities of any two distinct groups Ni and Nj of particles are correlated, 
and measure pij = 0. This information, on the other hand, is to some extent 
preserved in p BS , Pbq, and psq, i.e. the way in which quantum numbers are 
correlated, if at least occasionally a particle is detected during an experiment. 

The Monte Carlo scheme is employed further. In Section 6.1, the joint distri- 
butions of positively and negatively charged particles in momentum bins ApT t i 
and Ayi are constructed. Then, in turn, primordial and final state GCE results 
on the scaled variance u, Eq.(3.50), and the correlation coefficient p, Eq.(3.49), 
are extrapolated to the MCE limit in Section 6.2. To complement the above 
study of micro canonical effects, canonical (charge conservation) effects on final 
state multiplicity fluctuations are studied in Section 6.3. Further discussion of 
MCE effects is presented in Chapter 7. 
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6.1 Grand Canonical Ensemble 

In Fig. (6.1) the Apx,i (left) and Ayi (right) dependence of the GCE scaled 
variance u + of positively charged hadrons are shown, both primordial and fi- 
nal state. In the primordial Boltzmann case one finds no dependence of mul- 
tiplicity fluctuations on the position and size of the acceptance window. The 
observed multiplicity distribution is, within error bars, a Poissonian with scaled 
variance u + = 1. In fact, in the primordial GCE Boltzmann case any selection 
of particles has u — 1, as particle multiplicity, as well as momenta, are sampled 
independently. 
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Figure 6.1: GCE scaled variance u + of multiplicity fluctuations of positively 
charged hadrons, both primordial and final state, in transverse momentum 
bins ApT,i (^ e ft) and rapidity bins Ay^ (right). Horizontal error bars indicate 
the width and position of the momentum bins (And not an uncertainty!). Verti- 
cal error bars indicate the statistical uncertainty of 20 Monte Carlo runs of 2 ■ 10 5 
events each. The markers indicate the center of gravity of the corresponding bin. 
The solid line indicates the final state acceptance scaling estimate. 

Fig. (6.2) shows the Ap?,i (left) and Ay^ (right) dependence of the GCE cor- 
relation coefficient p + _ between positively and negatively charged hadrons, both 
primordial and final state. In the primordial Boltzmann case one finds also no 
dependence of multiplicity correlations on the position and size of the accep- 
tance window. The observed joint multiplicity distribution is a product of two 
Poissonians with correlation coefficient p + _ = 0. 

Fig. (6.3) presents the Ap T)i (left) and Ai/i (right) dependence of the GCE 
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Figure 6.2: GCE multiplicity correlations p + _ between positively and nega- 
tively charged hadrons, both primordial and final state, in transverse momentum 
bins ApT,i {left) and rapidity bins Ayi (right). The rest as in Fig. (6.1). 
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Figure 6.3: GCE scaled variance u± of multiplicity fluctuations of all charged 
hadrons, both primordial and final state, in transverse momentum bins Apx^ (left) 
and rapidity bins Ayi (right). The rest as in Fig. (6.1). 



scaled variance u± of all charged hadrons, both primordial and final state. As 
the correlation coefficient between positively and negatively charged hadrons 
is p + _ = in the primordial GCE one also finds u± = u + = 1. 

Resonance decay is the only source of correlation in an ideal GCE Boltzmann 
gas. Neutral hadrons decaying into two hadrons of opposite electric charge are 
the strongest contributors to the correlation coefficient p+-. The chance that 
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both (oppositely charged) decay products are dropped into the same momentum 
space bin is obviously highest at low transverse momentum (i.e. the correlation 
coefficient is strongest in Ap^i). The rapidity dependence is somewhat milder 
again, because heavier particles (parents) are dominantly produced around mid- 
rapidity and spread their daughter particles over a range in rapidity. One notes 
that the scaled variances and correlation coefficients in the respective acceptance 
bins in Figs. (6. 1,6. 2) are generally larger than the acceptance scaling procedure 1 
suggests, with the notable exception of p + _(Ap T 5 ). 

There is a simple relation connecting the scaled variance of the fluctuations 
of all charged hadrons u± to the fluctuations of only positively charged par- 
ticles u + via the correlation coefficient p + _ between positively and negatively 
charged hadrons in a neutral system: 



The contribution of neutral parent particles to oo± is stronger than to uj + , as both 
oppositely charged decay products go into the analysis rather than just one of 
them. Therefore, the effects of resonance decay on the Apx,i dependence of u± are 
considerably enhanced compared to the one of u + , and generally u± > u + , as the 
correlation coefficient p + _ remains positive in the final state GCE, in Fig. (6.3). 
Compared to this, the final state values of u±, u + and p + _ remain rather flat 
with Ay { in the GCE. 

If one would construct now a larger and larger number of momentum space 
bins of equal average particle multiplicities, one would successively lose more 
and more information about how multiplicities of distinct groups of particles are 
correlated, and approach the Poissonian limit also for final state particles. 

6.2 Micro Canonical Ensemble 

In the very same way in which fully phase space integrated extensive quantities 
were extrapolated to the MCE limit in Chapter 5, now multiplicity fluctuations u + 
and correlations p + _ in transverse momentum bins Apx,i an d rapidity bins Ay^ 
for a hadron resonance gas are extrapolated from the GCE (A = 0) to the MCE 

1 For the acceptance scaling approximation it is assumed that particles are randomly detected 
with a certain probability q — 0.2, independent of their momentum. 
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(A — > 1). Analytical primordial MCE results are obtained in the infinite volume 
approximation [106, 127], providing some guidance as to asses the accuracy of 
the extrapolation scheme. For final state fluctuations and correlations in limited 
acceptance, on the other hand, no analytical results are available. 

Mean values of particle numbers of positively charged hadrons (N + ) and neg- 
atively charged hadrons (A r „) in the respective acceptance bins, defined in Ta- 
ble 4.1, remain constant as A goes from to 1, while the variances ((AN + ) 2 ) 
and ((AN_) 2 ), and covariance (AN + AN^) converge linearly to their respective 
MCE limits. The correlation coefficient p + „ between positively and negatively 
charged hadrons, on the other hand, will not approach its MCE value linearly, as 
discussed in Chapter 5. 

Primordial 

In Fig. (6.4) the primordial scaled variance u + of positively charged hadrons in 
transverse momentum bins Apx,i {left) and rapidity bins Ayi (right) is shown as a 
function of the size of the bath A = Vi/V g , while in Fig. (6.5) the dependence of the 
primordial correlation coefficient p + „ between positively and negatively charged 
hadrons in transverse momentum bins Apx,i {left) and rapidity bins Ay, [right) 
is presented function of A. 

The results of 8 ■ 20 Monte Carlo runs of 2 • 10 5 events each are summarized 
in Table 6.1. The system sampled was assumed to be neutral rij = (0,0,0) and 
static = (1,0,0,0) with local temperature = 0.160 GeV and a system 
volume of V\ = 2000 fm 3 . 8 different values of A have been studied. The last 
marker (A = 1) denotes the result of the extrapolation. Only primordial hadrons 
are analyzed. Values for both Ap?,i and Ayi bins are listed. Analytical num- 
bers are calculated according to the method developed in [106, 127], using the 
acceptance bins defined in Table 4.1, and are shown for comparison. The effects 
of energy-momentum and charge conservation on primordial multiplicity fluctu- 
ations and correlations in finite acceptance will also be discussed in Chapter 7. 

Firstly fully phase space integrated results, also later discussed in Chapters 8 
and 9, are attended to. The scaled variance of multiplicity fluctuations is lowest 
in the MCE due to the requirement of exact energy and charge conservation, 
somewhat larger in the CE (see also Section 6.2), and largest in the GCE, as now 
all constraints on the micro states of the system have been dropped [101, 103, 105]. 
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Figure 6.4: Evolution of the primordial scaled variance u + of positively charged 
hadrons with the Monte Carlo parameter A = V\/V g for transverse momentum 
bins ApT,i {left) and for rapidity bins At/i (right). The solid lines show an analytic 
extrapolation from GCE results (A = 0) to the MCE limit (A — > 1). Each marker 
and its error bar, except the last, represents the result of 20 Monte Carlo runs 
of 2 • 10 5 events. 8 different equally spaced values of A have been investigated. 
The last marker denotes the result of the extrapolation. 




Figure 6.5: Evolution of the primordial correlation coefficient p + _ between posi- 
tively and negatively charged hadrons with the Monte Carlo parameter A = V\jV g 
for transverse momentum bins Apx^i (left) and for rapidity bins Ayi (right). The 
rest as in Fig. (6.4). 

The fully phase space integrated MCE and CE correlation coefficients between 
oppositely charged particles are rather close to 1. Doubly charged particles allow 
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for mild deviation, as also the A ++ resonance is counted as only one particle. 



Primordial 




Ap T ,2 


Apr,3 


Ap T ,4 


ApT,5 


, ,£fce 
w mce 

mce.c 
UJ + 


1.000 ±0.002 
0.889 ± 0.007 
0.8886 


1.000 ±0.002 
0.880 ±0.007 
0.8802 


1.000 ±0.002 
0.869 ±0.007 
0.8682 


1.000 ±0.002 
0.850 ±0.006 
0.8489 


1.000 ±0.002 
0.798 ±0.007 
0.7980 


qce 

p+- 

n mce 
mce.c 

p+- 


0.000 ±0.002 
0.094 ±0.005 
0.0935 


-0.000 ±0.002 
0.085 ± 0.006 
0.0844 


-0.000 ±0.002 
0.072 ± 0.006 
0.0730 


0.000 ±0.002 
0.056 ±0.006 
0.0554 


0.000 ±0.001 
0.003 ±0.005 
0.0040 


Primordial 


Aj/i 


Aj/2 


A2/3 


A2/4 


Ay 5 


U} + 
mce.c 


1.000 ±0.002 
0.795 ±0.006 
0.7950 


1.000 ±0.002 
0.835 ± 0.007 
0.8350 


1.000 ±0.003 
0.853 ±0.008 
0.8521 


1.000 ±0.002 
0.834 ± 0.006 
0.8351 


1.000 ±0.002 
0.794 ±0.007 
0.7949 


qce 

p+- 

n mce 

P+- 

mce.c 

p+- 


-0.000 ±0.001 
-0.013 ±0.005 
-0.0135 


0.000 ±0.002 
0.040 ±0.006 
0.0406 


0.001 ±0.002 
0.061 ±0.006 
0.0616 


0.000 ±0.002 
0.041 ±0.006 
0.0406 


-0.000 ±0.002 
-0.012 ±0.006 
-0.0135 



Table 6.1: Summary of the primordial scaled variance u + of positively charged 
hadrons and the correlation coefficient p + _ between positively and negatively 
charged hadrons in transverse momentum bins Apx,i and rapidity bins Arji. Both 
the GCE result (A = 0) and the extrapolation to MCE (A = 1) are shown. The 
uncertainty quoted corresponds to 20 Monte Carlo runs of 2 • 10 5 events (GCE) or 
is the result of the extrapolation (MCE). Analytic MCE results oj™ ce ' c and p™ c ^ ,c 
are listed too. 

The transverse momentum dependence can be understood as follows: a change 
in particle number at high transverse momentum involves a large amount of en- 
ergy, i.e., in order to balance the energy record, one needs to create (or annihilate) 
either a lighter particle with more kinetic energy, or two particles at lower px- 
This leads to suppressed multiplicity fluctuations in high Apx,i bins compared to 
low Apx,i bins. By the same argument, it seems favorable, due to the constraint 
of energy and charge conservation, to balance electric charge, by creating (or 
annihilating) pairs of oppositely charged particles, predominantly in lower Apj^i 
bins, while allowing for a more un-correlated multiplicity distribution, i.e. also 
larger net-charge (SQ = N + — N_) fluctuations, in higher Apx t i bins. 

For the rapidity dependence similar arguments hold. Here, however, the 
strongest role is played by longitudinal momentum conservation. A change in 
particle number at high y involves now, in addition to a large amount of energy, 
a large momentum p z to be balanced. The constraints of global P z conservation 
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are, hence, felt least severely around \y\ ~ 0, and it becomes favorable to bal- 
ance charge predominantly at mid-rapidity (p+_ larger) and to allow for stronger 
multiplicity fluctuations (cu + larger) compared to forward and backward rapidity 
bins. 

In a somewhat casual way one could say: events of a neutral hadron resonance 
gas with values of extensive quantities B, S, Q, E and P z in the vicinity of (Q\) 
have a tendency to have similar numbers of positively and negatively charged 
particles at low transverse momentum px and rapidity y and less strongly so at 
high pt and \y\. 

Final State 

Now the extrapolation of final state multiplicity fluctuations and correlations 
to the MCE limit is attended to. An independent Monte Carlo run for the same 
physical system was done, but now with only stable final state particles 'detected'. 

In Fig. (6.6) the final state scaled variance u + of positively charged hadrons 
in transverse momentum bins Apx^i {left) and rapidity bins Ayi (right) is shown 
as a function of A, while in Fig. (6.7) the dependence of the final state correla- 
tion coefficient p + _ between positively and negatively charged hadrons in trans- 
verse momentum bins Apx^i (left) and rapidity bins Ayi (right) on the size of the 
bath A = Vi/Vg is shown. 

The Apx ,i and Ayi dependence on A of the final state MCE scaled variance u + 
is qualitatively similar to that of the primordial versions, Fig. (6.4), and is essen- 
tially also explained by the arguments of the previous section. The effects of 
charge and energy-momentum conservation work very much the same way as 
before, and it still seems favorable to have events with wider multiplicity distri- 
butions at low pt and low y, and more narrow distributions at larger p? and 
larger \y\. The dependence of the final state correlation coefficients p + _ on A, 
Fig. (6. 7), is a slightly different to the primordial case, Fig. (6. 5). However, in the 
MCE limit, events still tend to have more similar numbers of oppositely charged 
particles at low p? and low y, than at large pt and large \y\. 

The effects of resonance decay are qualitatively different in the MCE, CE, 
and GCE (see also Section 6.3). Again, firstly attending to fully phase space in- 
tegrated multiplicity fluctuations discussed in [103, 101] and Chapter 9. The final 
state scaled variance increases in the GCE and CE compared to the primordial 
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Figure 6.6: Evolution of the final state scaled variance oo + of positively charged 
hadrons with the Monte Carlo parameter A = V\/V g for transverse momentum 
bins ApT,i {left) and for rapidity bins Ayi (right). The solid lines show an analytic 
extrapolation from GCE results (A = 0) to the MCE limit (A — > 1). Each marker, 
except the last, represents the result of 20 Monte Carlo runs of 2 • 10 5 events. 
8 different equally spaced values of A have been investigated. The last marker 
denotes the result of the extrapolation. 




Figure 6.7: Evolution of the final state correlation coefficient p + _ between posi- 
tively and negatively charged hadrons with the Monte Carlo parameter A = V\/V g 
for transverse momentum bins ApT,i (left) and for rapidity bins Ayi (right). The 
rest as in Fig. (6.6). 

scaled variance. Multiplicity fluctuations of neutral mesons remain unconstrained 
by (charge) conservation laws. However, they often decay into oppositely charged 
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particles, which increases multiplicity fluctuations of pions, for instance. In the 
MCE, due to the constraint of energy conservation, the event-by-event fluctua- 
tions of primordial pions are correlated to the event-by-event fluctuations of, in 
general, primordial parent particles, and u^ mal < oj prim is possible in the MCE. 
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Figure 6.8: MCE scaled variance u + of multiplicity fluctuations of positively 
charged hadrons, both primordial and final state, in transverse momentum 
bins ApTj {left) and rapidity bins Ay^ (right). Horizontal error bars indicate the 
width and position of the momentum bins (And not an uncertainty!). Vertical 
error bars indicate the statistical uncertainty quoted in Table 6.2. The markers 
indicate the center of gravity of the corresponding bin. The solid and the dashed 
lines show final state and primordial acceptance scaling estimates respectively. 



Figs. (6.8-6.10) compare the final state Apxj (left) and Ar/i (right) dependence 
of the MCE scaled variance of positively charged hadrons u + , the MCE correla- 
tion coefficient between positively and negatively charged hadrons p+~, and the 
MCE scaled variance of all charged hadrons u±, respectively to their primordial 
counterparts. The results of 8 • 20 Monte Carlo runs of 2 • 10 5 events each for a 
static and neutral hadron resonance gas with T = 0.160 GeV are summarized in 
Table 6.2. 

A few words to summarize Figs. (6.8-6. 10): resonance decay and (energy) con- 
servation laws work in the same direction, as far as the transverse momentum 
dependence of the scaled variances u>+, u± and the correlation coefficient p + _ 
is concerned. Both effects lead to increased multiplicity fluctuations and an in- 
creased correlation between the multiplicities of oppositely charged particles in 
the low pr region, compared to the high px domain. Compared to this, the 
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Figure 6.9: MCE multiplicity correlation coefficient p H between positively and 

negatively charged hadrons, both primordial and final state, in transverse mo- 
mentum bins ApT,i {left) and rapidity bins Ayi (right). The rest as in Fig. (6. 8). 
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Figure 6.10: MCE scaled variance u± of multiplicity fluctuations of all charged 
hadrons, both primordial and final state, in transverse momentum bins Apx,i (left) 
and rapidity bins Ayi (right). The rest as in Fig. (6.8). 



MCE Ayi dependence of uj + , p + _, and u>± is mainly dominated by global conser- 
vation of P z . Resonance decay effects, see also the GCE versions, Figs. (6. 1-6. 3), 
act more equal across rapidity, than in transverse momentum. 

Again, the scaled variance of all charged particles is found to be larger than 
the scaled variance of only positively charged hadrons oj± > u + , except for 
when p + _ < 0, i.e when the multiplicities of oppositely charged particles are 
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Final state 


A 11 rr 1 

1 — i y i , i 


'— *Pi .A 


A'f)T Q 


A?1t /I 
1 ,4 


A 11 rr e 
i — x t'l ,0 


u + 

. ,mce 

U/_l_ 
+ 


1.031 ±0.002 
0.904 ±0.007 


1.026 ±0.002 
0.884 ±0.007 


1.020 ±0.002 
0.872 ±0.007 


1.015 ±0.002 
0.847 ±0.007 


1.010 ± 0.002 
0.778 ± 0.006 


qce 

p+- 

„mce 
P+- 


0.163 ±0.001 
0.143 ±0.005 


0.107 ±0.001 
0.088 ±0.005 


0.109 ±0.001 
0.090 ± 0.005 


0.075 ±0.002 
0.049 ± 0.006 


0.052 ±0.002 
-0.010 ±0.006 


Final state 


Aj/i 




Ay 3 


Ay 4 


Ay 5 


, .gee 

UJ + 


1.017 ±0.002 
0.771 ±0.007 


1.023 ±0.002 
0.840 ±0.006 


1.024 ±0.002 
0.859 ±0.007 


1.023 ±0.003 
0.839 ± 0.007 


1.017 ±0.002 
0.770 ±0.006 


qce 

p+- 


0.100 ±0.001 
-0.027 ±0.005 


0.116 ± 0.001 
0.069 ±0.005 


0.115 ± 0.002 
0.092 ±0.006 


0.115 ±0.002 
0.069 ±0.006 


0.100 ±0.001 
-0.027 ±0.005 



Table 6.2: Summary of the final state scaled variance u + of positively charged 
hadrons and the correlation coefficient p + „ between positively and negatively 
charged hadrons in transverse momentum bins Ap^ j and rapidity bins A?/j. Both 
the GCE result (A = 0) and the extrapolation to the MCE (A = 1) are shown. 
The uncertainty quoted corresponds to 20 Monte Carlo runs of 2 • 10 5 events 
(GCE) or is the result of the extrapolation (MCE). 

anti-correlated, as for instance in Ap<r j5 , Ay 1; and Ay 5 . In contrast to that, how- 
ever narrowly, u± > 1 in the lowest transverse momentum bin Ap^i- Acceptance 
scaling provides a good estimate of the average in transverse momentum. In ra- 
pidity, however the approximation somewhat overshoots the limited acceptance 
results, pointing to correlations amongst momentum bins, Chapter 7. 

6.3 Canonical Ensemble 

To complement the discussion of multiplicity fluctuations and correlations 
for a hadron resonance gas in transverse momentum and rapidity segments, the 
canonical ensemble limit is studied next. Events will be re- weighted according to 
their extensive quantities baryon number B, strangeness S, and electric charge Q. 
Conservation laws for energy E and longitudinal momentum P z (as well as for 
transverse momenta P x and P y ) are hence not applied. Again some guidance for 
assessing the quality of the extrapolation is provided by analytical methods. Two 
independent (final state and primordial) Monte Carlo runs of 8 • 20 • 2 • 10 5 events 
each have been generated. The same system as before is considered and the same 
acceptance cuts are applied. 

Resonance decay is a process respecting the conservation laws for charge, en- 
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ergy and momentum. So despite the fact that charge conservation laws do not 
correlate primordial particle momenta, these additional correlations are intro- 
duced into the final state momentum space dependence, as in the GCE. 

In Figs. (6. 11-6. 13) the final state scaled variance u + of positively charged 
hadrons, the final state correlation coefficient p + _ between positively and neg- 
atively charged hadrons, and the final state scaled variance u± of all charged 
hadrons, in transverse momentum bins Ap^i {left) and rapidity bins Ay^ (right) 
are shown as a function of A. Primordial results are not included. 

Figs. (6.14-6.16) compare the final state Apx,i (left) and Ayi (right) dependence 
of the CE scaled variances u + , and u±, and the CE correlation coefficient p + _ 
respectively to their primordial counterparts. The results of 8 • 20 Monte Carlo 
runs of 2 • 10 5 events each for a static and neutral hadron resonance gas with 
temperature T = 0.160 GeV are summarized in Table 6.3. 
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Figure 6.11: Evolution of the final state scaled variance u + of positively charged 
hadrons with the Monte Carlo parameter A = V\/V g for transverse momentum 
bins Apx t i (left) and for rapidity bins Ay^ (right). The solid lines show an an- 
alytic extrapolation from GCE results (A = 0) to the CE limit (A — > 1). Each 
marker, except the last, represents the result of 20 Monte Carlo runs of 2 ■ 10 5 
events. 8 different equally spaced values of A have been investigated. The last 
marker denotes the result of the extrapolation. 

Firstly the primordial scenario is considered. Charge conservation does not 
introduce any correlations in momentum space. Primordial results are thus in- 
dependent of where in momentum space the acceptance window is located. Only 
the fraction of the system observed matters. One therefore finds the scaled vari- 



9s 



Chapter 6: Multiplicity Fluctuations and Correlations 



0.25 



< 0.15 



0.1 
0.05 


■0.05 




- » Apr.3 o Apn 
o Apr, 5 



0.25 



+ 







0.2 



0.4 



0.6 



0.8 



0.15 
0.1 
0.05 


■0.05 




r T At/3 % 







0.2 



0.4 



0.6 



0.8 



A 



A 



Figure 6.12: Evolution of the final state correlation coefficient p + _ between posi- 
tively and negatively charged hadrons with the Monte Carlo parameter A = V\/V g 
for transverse momentum bins Apx,i {left) and for rapidity bins Ayi (right). The 
rest as in Fig. (6. 11). 
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Figure 6.13: Evolution of the final state scaled variance u± between all charged 
hadrons with the Monte Carlo parameter A = V\/V g for transverse momentum 
bins ApT,i (left) and for rapidity bins Ay^ (right). The rest as in Fig. (6. 11). 



ance and the correlation coefficient the same in all Apr,« and Ayi bins. All 
values converge to the one suggested by the CE acceptance scaling approxima- 
tion with acceptance probability q = 0.2 (i.e. 5 bins of equal average particle 
number content). Fully phase space integrated values u and p can be found 
in [101, 103, 141, 142] and later on in Chapters 8 and 9. 
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The primordial scaled variance u + of positively charged hadrons drops lin- 
early to its CE limit as the size of the charge bath is reduced. Here one also 
finds w™ ce < u/j 6 , however almost equal in Ap^i, in limited acceptance. Com- 
bined energy and charge conservation leads to more constrained multiplicity fluc- 
tuations. The primordial correlation coefficient p+_ between positively and neg- 
atively charged particles, on the other hand, rises as the size of the bath de- 
creases, A — y 1. The correlation coefficient is indeed positive, as one would 
expect from charge conservation. In full phase space p+_ and p+i e are rather 
close to unity, as discussed earlier. However the correlation coefficient in limited 
acceptance appears weaker in the MCE than in the CE, p™ c ^ < p c f_. Energy 
conservation de-correlates the joint multiplicity distribution. In a primordial 
neutral CE system one finds the width of the distribution of all charged hadrons, 
Eq.(6.1), u± ~ 1 , i.e. essentially a Poissonian, in any momentum bin [93, 101]. 
Fluctuations of the number of all charged particles appear to be unconstrained 
by charge conservation and, thus, to yield the same result as in the GCE [94]. 
But also, like for u, the observed value is independent of the position and shape 
of the acceptance window. This is in contrast to primordial MCE calculations 
where the u and p show a distinctive dependence on the bins Ap T ^ and A?/;. See 
Figs. (6.4-6.5). 

No additional correlations are introduced into the momentum space depen- 
dence of the final state multiplicity distribution, apart from the well-known corre- 
lations due to resonance decay. The CE final state scaled variance u + in transverse 
momentum bins in Fig. (6. 14) (left) shows qualitatively a similar behavior as in 
the GCE, Fig.(6.1) (left), and in the MCE, Fig.(6.8) (left). The scaled variance 
decreases as charge conservation is turned on, however not as strongly as in the 
MCE. Neutral hadrons lead to enhancement mostly at low p?, hence u+ is larger 
in Apr,i than in Apr,5- As resonance decay works more evenly across rapidity, 
one finds the Ai/i dependence mostly flat in Fig. (6. 14) (right). 

The CE final state correlation coefficient p_| in transverse momentum bins 

in Fig. (6. 15) (left) qualitatively shows an equally similar behavior as in the GCE, 
Fig. (6.2) (left), and in the MCE, Fig.(6.9) (left). However with a notably stronger 
correlation at larger bins in Fig. (6. 15) (left) than in GCE and MCE. The ex- 
trapolations in Figs. (6.7) (left) to the final state MCE and in Fig. (6. 12) (left) to 
the final state CE correlation coefficient, show hence different behavior. Starting 
from the final state GCE values the correlation coefficient grows as the charge 
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Figure 6.14: CE scaled variance cu + of multiplicity fluctuations of positively 
charged hadrons, both primordial and final state, in transverse momentum 
bins ApT,i {left) and rapidity bins Ay^ (right). Horizontal error bars indicate the 
width and position of the momentum bins (And not an uncertainty!). Vertical 
error bars indicate the statistical uncertainty quoted in Table 6.3. The markers 
indicate the center of gravity of the corresponding bin. The solid and the dashed 
lines show final state and primordial acceptance scaling estimates respectively. 
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Figure 6.15: CE multiplicity correlation coefficient p + _ between positively and 
negatively charged hadrons, both primordial and final state, in transverse mo- 
mentum bins ApT.i (left) and rapidity bins Ayi (right). The rest as in Fig. (6. 14). 



bath is removed, but drops when both, charge and heat bath, are removed. Res- 
onance decay introduces a positive correlation. Charge conservation makes this 
trend stronger. In the MCE energy conservation works in the opposite direction. 
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Figure 6.16: CE scaled variance oj± of multiplicity fluctuations of all charged 
hadrons, both primordial and final state, in transverse momentum bins Apx,i (left) 
and rapidity bins Ayi (right). The rest as in Fig. (6. 14). 



Primordial 


Ap T ,l 


Ap T , 2 


Apr,3 


Ap T ,4 


Ap T ,5 


P+- 


0.901 ±0.005 
0.108 ±0.004 


0.904 ±0.006 
0.109 ±0.004 


0.901 ±0.006 
0.109 ±0.005 


0.901 ±0.006 
0.111 ±0.004 


0.898 ±0.006 
0.115 ±0.004 


Primordial 




Ay 2 


A^3 


Ay 4 


Ay 5 


ce 
UJ + 
n ce 
P+- 


0.903 ±0.005 
0.108 ±0.004 


0.900 ±0.006 
0.110 ±0.005 


0.900 ±0.005 
0.113 ±0.004 


0.902 ±0.005 
0.111 ±0.004 


0.903 ±0.006 
0.108 ±0.005 


Final state 




Ap T ,2 


ApT,3 


Ap TA 


Ap T ,b 


ce 
P%- 


0.989 ±0.006 
0.214 ±0.004 


0.973 ±0.006 
0.168 ±0.004 


0.962 ±0.005 
0.178 ±0.004 


0.948 ±0.006 
0.151 ±0.004 


0.922 ±0.005 
0.151 ±0.004 


Final state 




Ay 2 


Ay 3 


Ay 4 


Ay 5 


ce 
UJ + 

p C + e - 


0.956 ±0.006 
0.172 ±0.004 


0.963 ±0.006 
0.185 ±0.004 


0.960 ±0.005 
0.187 ±0.004 


0.960 ±0.005 
0.186 ±0.005 


0.955 ±0.005 
0.172 ±0.004 



Table 6.3: Summary of the final state and primordial scaled variance oo + of pos- 
itively charged hadrons and the correlation coefficient p + _ between positively 
and negatively charged hadrons in transverse momentum bins Apx^ and rapidity 
bins Ayi. Both the GCE result (A = 0) and the extrapolation to CE (A = 1) 
are shown. Analytical primordial CE values are o/j 6 = 0.9019 and p c ?_ = 0.1088, 
while final state acceptance scaling values are u+ = 0.9401 and p+_ = 0.1485. 
The uncertainty quoted corresponds to 20 Monte Carlo runs of 2 ■ 10 5 events 
(GCE) or is the result of the extrapolation (CE). 
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The At/i dependence is again mild in the CE, Fig. (6.9) (right). 

The final state CE values of u± in Fig. (6. 16) can be compared to their GCE 
counterparts in Fig. (6.3). As charge conservation does not introduce any corre- 
lations in momentum space, and as the effect of charge conservation is not seen 
in u±, also final state values are equal (within error bars) to the GCE results. 
GCE and CE values for u± are the same in full phase space for a neutral system 
(both primordial and final state). As charge conservation ensures (in the absence 
of energy conservation) a positive correlation coefficient p+_, one finds u + < u± 
in the CE in any momentum bin. The acceptance scaling estimate gives a good 
average in transverse momentum bins Ap^i, while under-predicting correlations 
and fluctuations in rapidity bins Ai/i, due to correlations introduced by resonance 
decay. 

6.4 Discussion 

GCE joint distributions of particle multiplicity were extrapolated to the MCE 
and CE limits. A few remarks attempt to summarize. The final state transverse 
momentum dependence Apx t i of multiplicity fluctuations and correlations are 
qualitatively similar in MCE, CE, and GCE. On the other hand, their (final 
state) rapidity dependence Ai/i exhibits only a very mild dependence in GCE 
and CE, which is in contrast to the MCE (with momentum conservation), where 
effects are strong. In the primordial case there is no momentum space dependence 
of multiplicity fluctuations and correlations in ensembles without energy conser- 
vation. Resonance decay trends are enhanced in the MCE. Charge conservation 
effects cease in a neutral system, and thus final state values of u± are the same in 
GCE and CE, even after resonance decay in limited acceptance, which provides 
a valuable cross check with analytical solutions. 

The extrapolation scheme works rather accurate and efficient. The statistical 
error on the 'data' points grows as A — > 1. The extrapolation helps greatly to 
keep the statistical uncertainty on the MCE (or CE) limit low, which can be 
seen from a comparison of the right most markers on the extrapolation lines. 
The last point and its error bar denote the result of a linear extrapolation of 
variances and covariances, while the second to last data point and its error bar 
are the result of 20 Monte Carlo runs with A = 0.875. The analytical MCE (or 
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CE) values are well within error bars of extrapolated Monte Carlo results, and 
agree surprisingly well, given the large number of 'conserved' quantities (5 and 3, 
respectively) and a relatively small sample size of 8 • 20 • 2 • 10 5 = 3.2 • 10 7 events. 
In a sample-reject type of approach this sample size would yield a substantially 
larger statistical error, as only events with exact values of extensive quantities are 
kept for the analysis. As the system size is increased, a sample-reject formalism, 
hence, becomes increasingly inefficient, while the extrapolation method still yields 
good results. Error bars are mildly smaller for an CE extrapolation than for an 
extrapolation to the MCE, due to a smaller number of conserved quantities. For 
further discussion see Appendices F and G. 

The qualitative picture presented in Fig. (6.8) could be compared to similar 
analysis of UrQMD transport simulation data [131], or recently published NA49 
data on multiplicity fluctuations in limited momentum bins [59]. Both (data and 
transport simulation) appear qualitatively similar to the MCE scenario, which 
could hint at a potentially strong role of conservation laws on fluctuation and 
correlation observables. The arguments used are general enough to believe they 
might hold as 'rules of thumb 1 in non-equilibrium situations too. If, in one event, 
one has a larger (smaller) number of particles in a momentum bin with larger 
and positive p z , than also a larger (smaller) number should go in the negative p z 
direction. This effect should be considerable as long as the momentum bins 
contain a not too small fraction of the particles of the whole event, say 20% as 
here. Or even smaller (11%) as in Chapter 7. Or yet smaller (< 10%) as in the 
UrQMD transport simulation [131]. 

It should be hard to assess the importance of conservations laws in simulations, 
let alone in data. This model is simple as one can follow things step by step 
and turn on and off certain effects. This is much harder in considerabely more 
complicated transport simulations, where individual sources of correlations can 
not that easily be pinned down. 

The systems to be studied are in general not neutral, so mean values of pos- 
itively and negatively charged hadrons are not the same, hence the spectra will 
also vary. A simultanious measurement of the momentum space dependence 
of u;_, u + , and p±, and their mean values (AL), and (N + ) should help to disen- 
tangle effects. For instance uj + > u± in a high momentum bin, while uj + < uj± in 
a low momentum bins is hard to do, except with energy conservation. 
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Chapter 7 

Micro Canonical Ensemble 



In this chapter the MCE is studied more closely, as it appears to be the most 
interesting amongst the standard ensembles. A simplified physical system is cho- 
sen to allow for smoother discussion. Owing to available analytical solutions [127 
Fermi-Dirac (FD) and Bose-Einstein (BE) effects are included into the analysis. 
The Monte Carlo scheme only facilitated Maxwell-Boltzmann (MB) statistics. 
For the examples a gas with three degenerate massive particles (with positive, 
negative and zero charge) with mass m = 0.140 GeV in three different statistics 
(MB, FD, BE) has been chosen. Thus examples, in particular the FD case, are a 
little academic in the sense that there is no fermion of this mass. In a hadron res- 
onance gas, discussed in the previous chapters, the lightest fermion is the nucleon 
for which quantum effects are probably negligible. 

In Section 7.1 multiplicity fluctuations and correlations in transverse momen- 
tum bins Apx^ and rapidity bins Ayi are discussed. Section 7.2 will turn to mul- 
tiplicity correlations between bins disconnected in momentum space (i/a and ys, 
or pt,a and Pt,b), and correlations between bins separated by some distance 4>g ap 
or Ugap in azimuth and rapidity. 

7.1 Fluctuations and Correlations within one 
Momentum Bin 

The properties of a static thermal system will be discussed first. Joint distri- 
butions of multiplicities N+ and iV_ are measured in limited bins of transverse 
momentum Ap T)i and of rapidity Ay^. Results will, in particular, be compared to 
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the acceptance scaling approximation employed in [94, 101, 103], which assumes 
random observation of particles with a certain probability q, regardless of particle 
momentum (see also Appendix E). Corresponding results for scaled variance in 
MB statistics can also be found in [106]. 

Static System 

The momentum spectra are assumed to be ideal GCE spectra due to the large 
volume approximation. In Fig. (7.1) transverse momentum and rapidity spectra 
are shown for MB statistics. BE and FD statistics yield similar spectra, unless 
chemical potentials are large. Momentum bins Apx,i and Ayi are indicated by 
drop-lines. 




Pr [GeV] y 



Figure 7.1: Differential particle spectra for a 'pion gas' at T = 0.160 GeV. Trans- 
verse momentum spectrum (left) and rapidity spectrum (right). Both curves are 
normalized to unity. The bins are constructed such that each bin contains 1/9 of 
the total yield. 

In Fig. (7.2) the scaled variance of positively charged particles u + is pre- 
sented within different transverse momentum bins Ap^^ (left) and rapidity 
bins Ayi (right). The scaled variance in limited bins of momentum space is more 
sensitive to the choice of particle statistics than the spectra would suggest. BE 
and FD effects are particularly strong in momentum space bins in which occupa- 
tion numbers are large. Hence, at the low momentum tail one finds suppression 
of fluctuations for FD and enhancement for BE, while at the high momentum 
tail, one finds ujbe — wmb — ^fd, Fig. (7.2) (left). 
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Figure 7.2: Transverse momentum (left) and rapidity dependence (right) of the 
MCE scaled variance of positively charged particles at T = 0.160 GeV, for a MB 
(blue), FD (green), BE (red) 'pion gas' at zero charge density. Momentum bins 
are constructed such that each bin contains the same fraction q of the average ir + 
yield. The horizontal bars indicate the width of the Apx t i or Ay^ bins, while the 
marker indicates the position of the center of gravity of the corresponding bin. 
Dashed lines indicate acceptance scaling results. 
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Figure 7.3: Same as Fig. (7.2), but for the MCE correlation coefficient between 
positively and negatively charged particles. Dashed lines indicate acceptance 
scaling results. 



The rapidity dependence, Fig. (7.2) (right), has a different behavior. The rea- 
son is, that in any Ay, bin there is some contribution from a low p? tail of 
the differential momentum spectrum dN/dy/dpT where quantum statistics ef- 
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fects are pronounced. This leads to a clear separation of the curves and one 
finds ujbe > wmb > ^fd- In contrast to this, the fully phase space integrated 
(all particles observed) scaled variance is rather insensitive to the choice of statis- 
tics [100] (unless chemical potentials are large). There are in fact three different 
'acceptance scaling' lines in Fig. (7.2), which extrapolate the fully phase space in- 
tegrated scaled variance to limited acceptance. The differences are however very 
small and all three lines lie practically on top of each other. 

Fig. (7.3) presents the correlation coefficient p + _ between positively and neg- 
atively charged particles in transverse momentum bins Apx,i {left) and rapidity 
bins Ayt (right). The fully phase space integrated correlation coefficient between 

positively and negatively charged particles would be p H = 1 in the CE and MCE. 

In the GCE it would be 0. In the MB CE it would not show any momentum 
space dependence and would always be p + _ > 0. In the MCE the situation is 
qualitatively different: in low momentum bins particles are positively correlated, 
while in high momentum bins they can even be anti-correlated. Horizontal lines 
again indicate acceptance scaling. Quantum effects for the correlation coefficient 
remain small as there is no explicit local (quantum) correlations between particles 
of different charge. 

It should be stressed that the ApT,i dependence in Figs. (7.2,7.3) is a direct 
consequence of energy conservation. The Ayi dependence of u + and p + _, how- 
ever, is due to joint energy and longitudinal momentum (P z ) conservation. Dis- 
regarding P z conservation leads to a substantially milder Ayi dependence, see 
Fig.(7.4). 

This behavior can be intuitively explained: in a low momentum bin it is com- 
paratively easy to balance charge, as each individual particle carries little energy 
and momentum. In contrast to this, in a high momentum bin with, say an excess 
of positively charged particles, it is unfavorable to balance charge, as one would 
also have to have more than on average negatively charged particles, and each 
particle carries large energy and momentum. This leads to suppressed fluctua- 
tions and correlations in high momentum bins when compared to low momentum 
bins. 

In a small mid-rapidity window, with \y\ < 0.3, the effects of globally applied 
motional conservation laws cease to be important (see Fig. (7.5)). Local correla- 
tions due to BE and FD statistics begin to dominate, and MCE deviations from 
the GCE results, Eq.(7.1), are relevant only for the highest momentum bins. In 
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Figure 7.4: Rapidity dependence of the scaled variance of positively charged 
particles (left) and of the correlation coefficient between positively and nega- 
tively charged particles (right). Calculations are done for the same system as in 
Figs. (7.2,7.3), however disregarding momentum conservation. 
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Figure 7.5: Transverse momentum dependence of MCE scaled variance of nega- 
tively charged particles (left) and the MCE correlation coefficient between posi- 
tively and negatively charged particles (right). Only particles in a mid-rapidity 
window \y\ < 0.3 are measured. Dashed lines denote the GCE result, Eq.(7.1). 
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low pt, where occupation numbers are large. MCE calculations in Fig. (7.5) are 
close to the GCE estimate Eq.(7.1). In MB statistics one finds only a weak Ap T 
dependence in a small mid-rapidity window. Please note that the acceptance 
scaling procedure predicts a Poisson distribution with u + — > 1 and p + _ — > for 
all three statistics in the limit of very small acceptance. 

Collectively Moving System 

As established a long time ago, in order to properly define the thermodynam- 
ics of a system with collective motion, the partition function needs to be Lorentz 
invariant [125, 126]. The expectation values of observables need hence to trans- 
form according to the Lorentz transformation properties of these observables. In 
particular, the inverse temperature /3 = T~ l is promoted to a four- vector /3 M 
(combining local temperature with collective velocity). The entropy, as well as 
particle multiplicities, remain Lorentz-scalars. 

These requirements are in general not satisfied unless momentum conservation 
is put on an equal status with energy conservation. If the system is described 
by a MCE, then momentum should be conserved as well as energy [125, 126]. If 
the system is exchanging energy with a bath, it needs to exchange momentum as 
well. 

For ensemble averages, neglecting these rules and treating momentum differ- 
ently from energy is safe as long as the system is close to the thermodynamic 
limit, since there ensembles become equivalent. The same is not true for fluc- 
tuation and correlation observables, which remain ensemble-specific [94]. For a 
system at rest, these requirements are not apparent since the net-momentum is 
zero. Statistical mechanics observables in a collectively moving system, however, 
lose their Lorentz invariance, if this is not maintained in the definition of the 
partition function. 

To illustrate this point, the properties of a system moving along the z-axis 
with a collective boost yo are considered. The total energy of the fireball is 
then E = M cosh(y ) ) while its total momentum is given by P z = M smh(y ). The 
mass of the fireball in its rest frame is M = P^u^. The system four-temperature 
is = ftu^. Local temperature and chemical potentials remain unchanged. The 
second rank tensor (or co- variance matrix) k 2 , Eq.(2.18), will be discussed in this 
section. 
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The second order cumulant k 2 , Eq.(2.18), is given by the second derivatives 
of the cumulant generating function with respect to the fugacities. Essentially 
this is the Hessian matrix of the function Eq.(2.8), encoding the structure of its 
minima. The diagonal elements Kj<2 £1X6 the variances of the GCE distributions 
of extensive quantities X. For example, k 2 a ' Na measures the GCE variance of 
the distribution of particle multiplicity of species A, while k®'® denotes the GCE 
electric charge fluctuations, etc. The off-diagonal k 2 ' elements give GCE co- 
variances of two extensive quantities X and Y. 

For a boost along the z-axis the general co-variance matrix for a relativistic 
gas with one conserved charge reads: 
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Off-diagonal elements correlating a globally conserved charge with one of the 

momenta, well as elements denoting correlations between different 

p p 

momenta, 1.6. t\j<2 ^ - 9X6 equal to zero due to antisymmetric momentum integrals. 
The values of elements correlating particle multiplicity and momenta, i.e. k 2 a,Px , 
depend strongly on the acceptance cuts applied. For fully phase space integrated 
multiplicity fluctuations and correlations these elements are equal to 0, again due 
to antisymmetric momentum integrals. 

It is instructive to review the transformation properties of k under the Lorentz 
group: k 2 ' contains the correlations between four-momenta P 11 and, in gen- 
eral, (scalar) conserved quantities and particle multiplicities QK Hence, the ele- 



ments K 



i.e. (AP^AP"), will transform as a tensor of rank 2 under Lorentz 
transformations; k p , i.e. (AP^AQi), will transform as a vector (the rapidity 
distribution will simply shift); and the remaining n® 3 '® will be scalars. 

For a static system the one finds for the co-variances of energy and mo- 
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Figure 7.6: Same as Fig. (7.2) (right) and Fig. (7.3) (right), but for a system moving 
with yo = 2. 




Figure 7.7: Same as Fig. (7.6), but without P z exact conservation. 



eigenvalues of the matrix Eq.(7.2) factorize, and momentum conservation can be 
shown to drop out of the calculation [106]. For a boost along the z-axis (and 



arbitrary particle acceptance) non- vanishing elements = «V^ Z 7^ appear, 

and ensure that the determinant of the matrix n 2 , Eq.(7.2), remains invariant 



e.p z 



Py,E 
l%2 



0. Transverse momenta P T 



against such a boost. Note that still 
and P y remain un-correlated with energy E. 

In Fig. (7.6) multiplicity fluctuations (left) and correlations (right) are shown 
for a system with boost x/q = 2. The rapidity spectrum of Fig.(7.1) (right) is 
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simply shifted to the right by two units. The construction of the acceptance 
bins is done as before. Multiplicity fluctuations and correlations within a ra- 
pidity bin Ayi transform as a vector (i.e., its z component shifts in rapidity) as 
inferred from their Lorentz-transformation properties, provided both energy and 
momentum along the boost direction are conserved. 

This last point deserves attention because usually, starting from [1], micro 
canonical calculations only conserve energy and not momentum. Imposing exact 
conservation for energy, and only average conservation of momentum will make 
the system non-Lorentz invariant, since in a different frame from the co- moving 
one energy and momentum will mix, resulting in micro state by micro state 
fluctuations in both momentum and energy. This result is not obvious, because 
energy-momentum is a vector of separately conserved currents. It is therefore 
natural to assume that these currents can be treated within different ensembles; 
they are, after all, conserved separately. It must be kept in mind, however, that 
it is not energy or momentum, but particles that are exchanged between the 
system and any canonical or grand canonical bath. The amount of energy and 
momentum carried by each particle are correlated by the dispersion relation [126]. 

In the situation examined here (unlike in a Cooper- Frye formalism [143], where 
the system is frozen out on a hyper surface, a space-time four-vector correlated 
with four-momentum) all time dependence within the system under considera- 
tion is absent due to the equilibrium assumption. Furthermore, the system is 
entirely thermal: the correlation between particle numbers when the system is 
sampled at different times is a 5— function, that stays a 5— function under all 
Lorentz-transformations. Hence, unlike what happens in a Cooper-Frye freeze- 
out, energy- momentum and space-time do not mix in the partition function. 
Together with the constraint from the particle dispersion relations, this means 
that different components of the four-momentum need to be treated by the same 
ensemble. 

This situation is explicitly shown in Fig. (7.7). Here, multiplicity fluctuations 
and correlations are calculated for the same system as in Fig. (7.6), but with exact 
conservation of only energy (and charge), but not momentum. In the co-moving 
frame of the system, the fluctuations and correlations are identical to Fig. (7.4). 
When the system is boosted the distribution changes (not only by a shift in 
rapidity, as required by Lorentz-invariance), and loses its symmetry around the 
system's average boost. 
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This last effect can be understood from the fact that momentum does not have 
to be conserved event-by-event, but energy does. It is easier, therefore, to create 
a particle with less longitudinal momentum (energy) than the average, than with 
larger longitudinal momentum (energy), and still conserve energy overall. This 
leads to suppressed multiplicity fluctuations and a negative correlation coefficient 
for rapidity bins in the forward direction in comparison to rapidity bins in the 
backward direction. In Fig. (7.6), where the system also needs to conserve mo- 
mentum exactly, this enhancement is balanced by the fact that it will be more 
difficult to conserve momentum globally when particles having less momentum 
than the boost are created. 

7.2 Correlations between disconnected Momen- 
tum Bins 

Next, correlations between particles in different momentum bins are stud- 
ied. 'Long range correlations' between bins well disconnected in momentum 
have been suggested to arise from dynamical processes. Examples include color 
glass condensate [50, 51], droplet formation driven hadronization [144], and 
phase transitions within a percolation-type mechanism [50, 145]. The ellip- 
tic flow measurements, widely believed to signify the production of a liquid at 
RHIC [10, 146, 147, 148], are also, ultimately, correlations between particles dis- 
connected in phase space (here, the azimuthal angle). 

As will be shown, conservation laws will also introduce such correlations be- 
tween any two (connected or not) distinct regions of momentum space. No dy- 
namical effects are taken into consideration (only an isotropic thermal system). 

Considering correlations between the multiplicities of particles A and B, 
within two bins, each centered around the rapidities y A and ys, with (constant) 
widths Ay A = Ay B = 0.2, in Fig. (7.8) (left) the correlation coefficient p + _ 
between positively and negatively charged particles are shown as a function 
of ua and x/b- Fig. (7.8) (right) presents the correlation coefficients between like- 
charge unlike-charge and all charged particles p±±, as a function of 
the separation of the two bin centers y gap = Da — Ub- m Fig- (7-9) the correla- 
tion coefficients p ++ (left) and p±± (right) are shown as a function of i/a and x/b- 
The main diagonal contains multiplicities in the same momentum bin. In the 
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case of p ++ and p±± this leads to a double counting of each particle and hence 
to p = 1. 

Energy conservation always leads to anti-correlation between different mo- 
mentum space bins. Charge conservation leads to a positive correlation of 
unlike charged particles and anti-correlation of like-sign particles. Longitudi- 
nal momentum conservation, however, is responsible for the structure visible in 
Fig. (7. 8) (left). Having a small (large) number of particles of one charge in a bin 
with positive average longitudinal momentum, leads to a larger (smaller) number 
of particles of opposite charge in a bin with different but also positive P z , (blue 
dips). This makes also a state with smaller (larger) particle number with oppo- 
site average longitudinal momentum — P z more likely (red hills). At large values 
of Da the correlation coefficient p ~ for any y B , because the yield (Na) in Ay A 
is asymptotically vanishing. 




Figure 7.8: (Left): The correlation coefficient p + _ between the multiplicity of 
positively charged particles in a bin located at x/a with negatively charged particles 
in a bin located at y B , both with a 0.2 width in rapidity. (Right): The correlation 
coefficients p ++ , and p±± between multiplicities in two bins separated by y gap 
of like, unlike, and all charged particles. Both plots show MCE MB results. 

The correlation coefficients p ++ , Fig. (7. 9) (left), and p±±, Fig. (7. 9) (right), 
show strong anti-correlations around center-rapidity, as to be expected from 
charge, energy and momentum conservation. In the case of p ++ the correla- 
tion stays negative. For p±± one observes a positive correlation between the 
numbers of all charged particles going in opposite directions, while seeing an 
anti-correlation for particle multiplicities measured in the same direction. This 
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Figure 7.9: The correlation coefficient p ++ between the multiplicity of positively 
charged particles in a bin located at ua with positively charged particles in a bin 
located at ys {left), and the correlation coefficient p±± between the multiplicity 
of all charged particles in a bin located at x/a with all charged particles in a bin 
located at ys {right). Both plots show MCE MB results. 

behavior is similar to that in Fig. (7.8) (left). For p + _ energy, momentum, and 
charge conservation favor pairs of oppositely charge particles in opposite direc- 
tions over particles in the same direction. For like sign particles the correlation 
stays negative, however milder for bins located in opposite direction. For not 
charge separated multiplicities the effect of charge conservation cancels out and 
energy and momentum conservation determine the correlation function. 

In Fig. (7.8) (right) the correlation coefficient is presented along the diagonal 
from top left to bottom right as a function of separation y gap = yA~ Ub- Unlike- 
sign particles are positively correlated. Like-sign and all charged particles are 
negatively correlated at small separation y gap . For large separation the correlation 
becomes asymptotically zero, because the yield is zero. However, please note that 
in particular p±± > p + „ at large y gap - Here P z conservation is indeed dominant. 

Fig. (7. 10) shows the correlation coefficient p + _ between positive and negative 
particles as a function of yA and ys, this time disregarding momentum conserva- 
tion, for a static source (left) and a source with a collective boost of y$ = 2 (right). 
Disregarding P z conservation destroys the particular structure in Fig. (7.8) (left) 
and leads to a single peak at the origin. The correlation is then insensitive to the 
momentum direction, and only sensitive to the energy content of a bin Ay. The 
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Figure 7.10: The correlation coefficient p + _ between the multiplicity of positively 
charged particles in a bin located at i/a with negatively charged particles in a bin 
located at ys, without conservation of longitudinal momentum P z , for a static 
source (left) and a source with a collective boost of yo = 2 (right). 

observables in Figs. (7.8,7.9) do not, Fig. (7. 10) (right), transform under boosts 
(yA,B Va,b — Uo)i unless momentum along the boost axis is exactly conserved. 

To complement the discussion of correlations across rapidity intervals, corre- 
lations between two distinct transverse momentum bins are studied. Here the 
discussion can be reduced to energy and charge conservation. Fig. (7. 11) shows 
the correlation coefficient p + _ between the multiplicities of positively charged 
particles in transverse momentum bins of equal size (Apx = 0.2 GeV) centered 
around Pt,a with the one of negatively charged particles in transverse momen- 
tum bins located at pt,b (left), and the correlation coefficient p ++ between the 
multiplicities of positively charged particles in transverse momentum bins located 
at Pt,a with positively charged particles in bins located at Pt,b (right). 

From charge conservation alone it follows that particles with electric charge 
of unlike sign are positively correlated, while the multiplicities of particles of the 
same charge are negatively correlated. The main diagonal in Fig. (7. 11) (right) is 
excluded, as there p ++ = 1 by construction. In Fig. (7.11) (left) there is an island 
of weak anti-correlation. The correlation function has its maximum close to the 
peak of the transverse momentum spectrum, i.e. where the average particle yield 
is largest. In these bins, and charge, and in particular, energy conservation have 
strong effects. 
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Figure 7.11: (Left:) The correlation coefficient p + _ between the multiplicity of 
positively charged particles in a bin located at pt,a with negatively charged parti- 
cles in a bin located at pt,b- (Right:) The correlation coefficient p ++ between the 
multiplicity of positively charged particles in a bin located at pt,a with positively 
charged particles in a bin located at Pt,b- 

Angular correlations, lastly, could arise due to conservation of transverse mo- 
menta P x and P y . In Fig. (7. 12) the correlation coefficient between particles in 
different A0 bins are shown. The flat 1 angular spectrum dN/dcj) has been divided 
into 10 equal size bins and the correlation coefficient is presented as a function of 
separation of the centers of the corresponding bins. 

Considering Fig. (7. 12) one firstly notes that when disregarding exact conser- 
vation of P x and P y the correlation coefficients are insensitive to the distance 4>g ap 
of any two bins. Only the correlations due to energy and charge conservation 
affect the result. Charge conservation leads to correlation of unlike-sign particles 
and to anti-correlation of like-sign particles. Energy conservation always anti- 
correlates multiplicities in two bins. For p±± the effect of charge conservation 
cancels for a neutral system, however, effects of energy-momentum conservation 
are strong 61*. SbS £L larger number of particles (hence a larger part of the total 
system) is observed. 

Conservation of transverse momenta P x and P y is now responsible for the (f) gap 
dependence of p. The line of arguments is similar to the ones before: Observing 
a larger (smaller) number of particles in some bin at 0o implies that, in order 

Only globally equilibrated systems are considered, and elliptic flow is disregarded here. 
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Figure 7.12: The correlation coefficients of particles in distinct A<ft bins as a 
function of separation <p gap in azimuth. {Left:) integrated over all phase space. 
[Right:) only particles with \y\ < 0.3 are observed. Both plots show MCE MB 
results. No elliptic flow is considered. 

to balance momenta P x = P y = 0, one should also observe a larger (smaller) 
number of particles in the opposite direction n — O - A larger (smaller) number 
of particles in a bin with <ft gap = n/2 would do little to help to balance momentum, 
but conflict with energy conservation. 

The observable effects of energy-momentum conservation are weaker, if the 
experimental acceptance only covers a mid-rapidity region of \y\ < 0.3, as in 
Fig.(7.12) (right), rather than the whole rapidity distribution, Fig. (7.12) (left). 

7.3 Discussion 

In this chapter multiplicity fluctuations and correlations in limited momen- 
tum bins for ideal relativistic gases have been discussed in the MCE in the ther- 
modynamic limit. For the discussion a gas with three degenerate massive parti- 
cles (positive, negative, neutral) in three different statistics (Maxwell-Boltzmann, 
Fermi-Dirac, Bose-Einstein) was chosen. 

For the width of multiplicity distributions in limited bins of momentum space 
a simple and intuitive picture emerges. In the Maxwell-Boltzmann approxima- 
tion one finds a wider distribution for momentum bins with low average momen- 
tum when compared to bins with higher average momentum but same average 



7.3 Discussion 



119 



particle number. This qualitative behavior is a direct consequence of energy and 
momentum conservation. The results in Fermi-Dirac and Bose-Einstein statistics, 
furthermore, show pronounced effects at the low momentum tail of the spectrum. 

The correlation coefficient additionally shows a similar qualitative behavior. 
In bins with low average momentum the correlation coefficient between positively 
and negatively charged particles is indeed positive, as one would expect from 
charge conservation. However, in bins with large average momentum the effects of 
joint energy and momentum conservation can lead to anti-correlated distributions 
of unlike-charged particles. 

The role of exactly imposed motional conservation laws is particularly impor- 
tant for systems with collective velocity. Fluctuations and correlations transform 
under boosts, provided momentum conservation along the boost direction is taken 
into account. This ensures, in particular, that they become boost-invariant if the 
underlying system is boost-invariant. 

Lastly, it was found that even in the thermodynamic limit long range corre- 
lations between disconnected regions in momentum space prevail. Multiplicities 
in different bins in rapidity, transverse momentum, or in azimuth, can have a 
non-zero correlation coefficient. 
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Chapter 8 

The Phase Diagram 

In this chapter the statistical properties of a hadron resonance gas are studied 
in their dependence on 'freeze-out' parameters, and choice of ensemble. Fully 
phase space integrated extensive quantities are investigated, omitting limited ac- 
ceptance effects, but still considering resonance decay. 

Preceding chapters were mostly concerned with neutral and static gases. At- 
tention was given, in particular, to the consequences of limited acceptance in 
momentum space. Here, fully phase space integrated fluctuations and correla- 
tions observables for a hadron resonance gas will be studied in their dependence 
on energy and net-baryon density. For simplicity, remaining chemical potentials 
are chosen to be zero. Thus, mesons and anti-mesons contribute equally to the 
partition function. But as, for instance, the proton carries electric charge, and 
the A carries strangeness, the system is only for \xb = 0, i.e. equal baryon and 
anti-baryon contribution, also electric charge and strangeness neutral. For pos- 
itive \ib, while fiQ = fis = 0, the system has positive electric charge density 
and negative strangeness density. For this discussion the phase space occupancy 
factors [116, 117, 118] are taken to be 7 S = 7 9 = 1. Nevertheless, the constant 
average energy per average particle multiplicity (E) / (N) ~ 1 GeV freeze-out 
line [20] is shown for orientation. 

Above this freeze-out line the following discussion should be taken with a 
caveat mentioned: First and foremost, a phase transition is expected to occur. 
Degrees of freedom should hence be deconfined quarks and gluons, not hadrons. 
Furthermore, even on the hadronic side, inclusion of Hagedorn states might be- 
come important close to the phase transition line [149, 150]. The exponentially 
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rising Hagedorn mass spectrum [83] then would lead to a diverging partition func- 
tion, above a certain limiting temperature, and the phase is not even defined. This 
discontinuity could in principle be fixed [151, 152], by allowing for a transition to 
a bag model phase. Quantum number configurations or decay channels of these 
Hagedorn states or bags 1 are however unknown, and would have to be assumed. 
The chemical freeze-out line could be seen as an estimate of where in the phase 
diagram these effects would become important. 

It is stressed, again, that no hadronic phase should exist above the freeze-out 
line. It should also be mentioned that, in principle, the intensive variables T and 
fiB ar e an unconventional choice to present MCE (or CE) correlation functions. 
Nevertheless, this choice of presentation was made, to allow for a comparison of 
statistical properties of different ensemble formulations at different energy and 
net-baryon densities. To highlight certain aspects, phase diagrams are extended 
to negative values of \1b- 

In the CE or MCE charge correlation and fluctuations would be vanishing in 
full acceptance, since charges (or additionally energy and momentum) are exactly 
conserved in these ensembles. As discussed in Chapter 5, charge fluctuations are 
determined by the thermodynamic bath available, and remain, full acceptance 
assumed, unaffected by resonance decay. Charge, energy and momentum conser- 
vation is respected by each decay process (weak decays omitted). In contrast to 
this, multiplicity correlation coefficients in the primordial GCE are equal to zero, 
while multiplicity fluctuations are Poissonian in the Boltzmann approximation. 
Correlations in the GCE appear due to resonance decay, and are strongly mod- 
ified, along with a suppression of multiplicity fluctuations, by conservation laws 
in the CE or MCE, see Chapters 6 and 7. 

Acceptance effects are neglected here. The correlation coefficients and vari- 
ances studied here, with full acceptance assumed, would also exhibit the quali- 
tative behavior presented in Chapters 4 for charge correlations, and Chapters 6 
and 7 for multiplicity fluctuations and correlations. 

The T-fis phase diagram will be explored for GCE charge correlations and 

fluctuations in Section 8.1, while in Sections 8.2 and 8.3 multiplicity correlations 

and fluctuations will be discussed in different ensembles. 

■"■The terms Hagedorn state, string, cluster, bag, or resonance are used equivalently in this 
chapter. The only property they are required to have is that they decay into hadrons. 
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8.1 Charge Correlations and Fluctuations 

In Fig. (8.1) grand canonical correlation coefficients between charges are shown 
in their dependence on the freeze-out parameters T and ps- The correlation 
coefficients pbs, Pbq-, and psq, as well as related variances and covariances are 
symmetric around ps — in the T-ps phase diagram. The sum of baryon 
and anti-baryon yield is equal in systems with baryon chemical potentials of 
opposite sign, since here additionally remaining chemical potentials are chosen 
to be p s = p Q = 0. Further, one finds (ABAS), (ABAQ), ((AQ) 2 }, etc. to 
be the same for ps and ~Pb, due to equal contributions of particles and their 
anti-particles. The meson contribution then only depends on the temperature. 
The discussion will proceed by considering the four corners of the phase diagrams 
depicted in Fig. (8.1), starting from the bottom left. 

At low T and ps the systems is mostly composed of pions. Accordingly 
one finds the correlation coefficients pss — Pbq — Psq — in the bottom left 
corners of phase diagrams shown in Fig. (8.1). The system is too cold to populate 
states with hadrons which could mediate the correlation. Pions carry only electric 
charge, and thus cannot correlate to the strangeness and baryon number content 
of the system. 

As T is increased, but still at small pis, the system starts to equally excite 
baryon and anti-baryon states. The correlation between baryon number and 
strangeness content of the system, pss, therefore grows negative due to emerg- 
ing A and A states. Compared to pbs, the baryon number electric charge cor- 
relation coefficient pbq grows more mildly. Protons and anti-protons are being 
equally produced, yet abundant pions de-correlate the quantum numbers B and Q 
more strongly than heavier kaons can de-correlate the quantum numbers B and S. 
Lastly, in the top left corner of Fig. (8.1) (right), the strangeness electric charge 
correlation coefficient psq has a peak. All mesons and baryons (except the E + ) 
contribute positively to the covariance (ASAQ). The combined contribution of 
particles and anti-particles at a given temperature is strongest in charge neu- 
tral matter. Only at rather large temperature heavy S + resonances (and their 
anti-particles) start to again decrease psq- 

The bottom right corners in Fig. (8.1) represent a hadronic phase of low tem- 
perature and large net-baryon density. The correlation coefficient of baryon num- 
ber and strangeness pss grows in this direction due to A production. The system 
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Figure 8.1: Phase diagram showing grand canonical charge correlation coeffi- 
cients of baryon number and strangeness pss (left), baryon number and electric 
charge Pbq (center), and strangeness and electric charge Psq (right). The solid 
line indicates the (E)/(N) = 1 GeV chemical freeze-out line [20]. 
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Figure 8.2: Phase diagram showing normalized grand canonical charge fluctua- 
tions of baryon number ((AB) 2 ) (left), strangeness ((AS) 2 ) (center), and electric 
charge ((AQ) 2 ) (right). The color scale is logarithmic. 

is on the other hand is too cold to produce kaons, which could de-correlate the 
quantum numbers B and S, in sufficient numbers. As protons (and neutrons or 
light As) are the dominant particles in this phase, the correlation between baryon 
number and electric charge pbq is particularly strong. And, as kaons are not suf- 
ficiently produced at low temperature (and the A is electrically neutral), p$Q, the 
strangeness electric charge correlation coefficient, remains small at small T. The 
increase of psq for large pb is due to the onset of £ production. 
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Finally, moving into the top right corners in Fig. (8.1), with the above men- 
tioned caveat in mind. The correlation pbs becomes strongly negative in hot net- 
baryon rich matter, as strangeness carrying baryons are produced abundantly. 
The correlation pbq stays modest as heavy negatively charged baryons emerge 
too, and like pions, de-correlate the systems electric net-charge from its net- 
baryon number content. Much for the same reasons psq stays modest, too. Here 
heavy S + resonances decrease the correlation. 

In Fig. (8.2) grand canonical charge fluctuations are shown in their depen- 
dence on the freeze-out parameters T and /ig on a logarithmic color scale. The 
variance is, unlike the correlation coefficient, an extensive quantity. The vari- 
ances are therefore shown normalized to an unit volume. In a neutral system, 
with symmetric particle and anti-particle contribution to the partition function, 
charge fluctuations are largest. With increasing temperature, not only average 
occupation numbers, but also fluctuations in occupation numbers on individual 
momentum levels grow. Accordingly, net-charge fluctuations become stronger. 
For increasing fj,g, i.e. average baryon occupation numbers are enhanced, while 
occupation numbers of anti-baryons are suppressed by Boltzmann factors e ±Ats/3 , 
charge fluctuations decrease. At low temperature, few particles are produced, and 
the net-charge content of the system cannot fluctuate much. As p, B is increased at 
low T, baryon states are populated due to charge density. The net-charge content 
of the system can then fluctuate, by virtue of iib induced particle density. 

The baryon number variance ((AB) 2 ) in Fig. (8.2) (left) increases rapidly with 
baryon density and temperature. At low T and \xb essentially no baryon states 
are occupied, and hence their number cannot fluctuate strongly. The strangeness 
variance ((AS*) 2 ) in Fig. (8.2) (center) and electric charge variance ((AQ) 2 ) in 
Fig. (8.2) (right) are similar. Strangeness and electric charge density are intro- 
duced indirectly via the baryons p and A. Although the strangeness and elec- 
tric charge variances are initially larger, they grow more slowly than the baryon 
number variance. The mass difference between the lightest baryons and mesons 
becomes less important at very large temperature, and ((AB) 2 ) eventually over- 
takes ((AS) 2 ) and ((AQ) 2 ). 

In direction from bottom left to top right in Figs. (8. 1,8. 2), i.e. increasing 
both T and a transition from a mesonic to baryonic matter [153] occurs. The 
variance ((AB) 2 ) and the correlation coefficient pss grow as the entropy of the 
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system is increasingly carried by baryons. If one would instead follow the (E) / (N) 
freeze-out line one would observe a monotonous decrease in pbq as colliding beam 
energy is increased, due to emerging pion production de-correlating the charges B 
and Q. The correlation between S and Q, psQ, on the other hand, increases, 
because kaons, and not A, now carry the bulk of the system strangeness at large 
center of mass energy. Along the freeze-out line pss has a strong minimum 
around T ~ 0.13 GeV and \ib ~ 0.5 GeV. The strangeness and electric charge 
variances ((AS*) 2 ) and ((AQ) 2 ) increase along this line, while {(AB) 2 ) stays more 
or less flat. 

8.2 Multiplicity Correlations 

In this section the correlation coefficients p pn + and p pn ~ between the event-by- 
event multiplicities of protons, N p , and either of the charged pions, N n + and N n - 
are studied for different ensembles. The discussion is then extended to multiplic- 
ity fluctuations of positively and negatively charged hadrons u + and The 
hadrons p and tt^ 1 are the most abundant charged particles, and their relative 
yield determines structure of the T-ps phase diagram to a large extent. In turn, 
GCE, CE, and MCE correlation functions will be compared. Resulting (primor- 
dial) correlations are not due to local interaction amongst constituents, but due 
to globally implemented conservation laws for energy and charge. 

Grand Canonical Ensemble 

The multiplicities of any two distinct groups of primordial particles are uncor- 
rected in the grand canonical ensemble, p = 0, due to the infinite thermodynamic 
bath assumption. Their joint distribution factorizes into a product of two Pois- 
sonians with scaled variance u = 1. In the final state grand canonical ensemble 
the proton and charged pion multiplicities are correlated by the decay of parent 
baryons N. The following decay channels, containing protons and charged pions, 
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are considered: 



N~ < — > p + n~ 

iV° < — > p + ?T 

N + < > p + 7T^ 

N ++ < > p + 7T H 



7T 



7T 



(8.1) 
(8.2) 
(8.3) 
(8.4) 



Baryonic resonances can decay into protons and charged pions taking channels 
of kind (8.1) to (8.4). Negatively charged baryons N~ , could only decay into 
negatively charged pions and a proton via (8.1). The A~ resonances, however, 
mostly decay via A - — > n + 7r~. Strangeness carrying N~ also cannot strongly 
populate the channel (8.1), as either daughter particle, baryon or meson, has to 
pick up the strange quark. The channel (8.2) for decay of neutral baryons N° 
is also not too efficient. Two thirds of neutral baryons decay via N° — > n + 7r° 
due to "Clebsch-Gordan" coefficients. Although most low lying A + decay via 
the channel (8.3), A + — > p + 7r°, and do hence not contribute much to either 
correlation coefficient, heavier A + resonances, often produce several pions, and 
positively correlate p with both, ir + and 7r~. Finally, channels (8.4) are abundant 
due to a strong contribution of the A ++ (1232) resonance and its degenerate states. 
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Figure 8.3: Final state GCE correlation coefficient between the multiplicities of 
protons and positively charged pions p pvr + (left), and between the multiplicities 
of protons and negatively charged pions p pvr - (right). 
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The reverse arrows indicate back-reactions. The strong interaction between 
the ground state hadrons proton, neutron, and pion is conceptually taken into 
account in the hadron resonance gas model by inclusion of resonance states which 
mediate the interaction [154, 155]. In thermal and chemical equilibrium reaction 
rates are the same in both directions. The system is said to have attained "de- 
tailed balance". Only the existence of resonances states is relevant, and not how 
many particles are involved in each microscopic reaction. The ability of a non- 
equilibrium state to evolve into a chemical equilibrium state, on the other hand, 
is strongly increased by multi-hadronic (or Hagedorn) channels [150]. 

Hence, in the final state GCE, Fig. (8.3), both, ir + and tt~ are positively 
correlated with p. Phase diagrams showing multiplicity correlations are generally 
not symmetric around \Ib — 0. The correlation is particularly strong in hot 
baryonic matter. For negative /Ib either correlation coefficient p p7r + ~ p p7r - ~ 0, 
as protons are then suppressed compared to anti-protons. Charged pions would 
now be strongly correlated with the more abundant anti-protons. 

The fact that n + are more strongly correlated via resonance decay with p, 
than the 7r~, i.e. p pw - < p p7r +, can be understood considering reaction channels 
(8.1) to (8.4). Essentially one may add additional pions (ir + + n~) to each 
channel. Channels (8.1) and (8.2) are not very effective at correlating p and 7r~, 
since they are not very abundant. Channel (8.3) contributes equally to the p-7r _ 
correlation and the p-7r + correlation, while channel (8.4) strongly correlates p 
and 7r + , yet hardly p and n~ . The influence of the A ++ (1232) can clearly be seen 
in a comparison of the bottom right corners of Fig. (8.3) (left) to (right). 

Canonical Ensemble 

The event-by-event multiplicities of p and ir + are anti-correlated in the pri- 
mordial CE, Fig. (8.4) (left), while the multiplicities of p and ir~ are positively 
correlated, Fig. (8.5) (left), due to electric charge conservation. Additionally, mul- 
tiplicity fluctuations of either species are suppressed in the CE (see also Chap- 
ter 9). 

The multiplicities of the hadrons p and ir + are most strongly anti-correlated 
at low T and large u-b- The multiplicity of baryons cannot fluctuate strongly 
in cold net-baryon rich matter, due to suppressed anti-baryon states, compare 
Fig. (8.2) (left). Charge conservation now requires that once one removes (adds) 
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He [GcV] h b [GeV] 

Figure 8.4: CE correlation coefficient between the multiplicities of protons and 
positively charged pions p pn +. Both, primordial (left), and final state (right). 




Figure 8.5: CE correlation coefficient between the multiplicities of protons and 
negatively charged pions p p7r -. Both, primordial (left), and final state (right). 



a proton from (to) the system, to add (remove) a neutron and a positively charged 
pion, p -H- n + 7r + . The system is too cold to allow for strong baryon multiplicity 
fluctuations 2 . In a neutral, but hot, hadron resonance gas the anti-correlation 
is weaker, since channels involving anti-baryons are open. In a hot and baryon 
rich phase, the appearance of heavy A resonances makes the p-7r channel less 

2 Encrgy fluctuations in the CE are also suppressed in comparison to the GCE [93]. 
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influential, and accordingly p pn + ~ there. 

For the final state CE correlation coefficient p p7T +, in Fig. (8.4) (right), the 
anti-correlation between p and n + extents to high \xb and T. Primordial parent 
particles, in particular positively charged A resonances, are anti-correlated, due 
to baryon number and electric charge conservation, with primordial p and 7r + . 
These A decay into both, p and ir + , and partially re-correlate their multiplicities, 
yet an residual anti-correlation remains. At low T and low u< B the correlation is 
then weaker, since the system is too cold to produce sufficient baryonic resonances 
which could correlate p and 7r + via their decay, i.e. the primordial and the final 
state scenario p pn + are similar in this region of the phase diagram. At high T, 
but hb — 0, the final state correlation coefficient is positive, p pn + > 0, unlike 
one would expect from electric charge conservation. In a neutral and hot hadron 
resonance gas, baryon and anti-baryon number fluctuations are sufficient, to de- 
correlate primordial p and tt + . The positive resonance decay contribution causes 
then an overall positive correlation. In final state GCE one always finds p p7T + > 0. 

In the primordial CE, one finds the correlation coefficient between the multi- 
plicities of p and 7T _ to be positive, p pw - > 0, Fig. (8.5) (left), due to electric charge 
conservation. Their correlation is particularly strong for low T and large ri B , since 
there baryon multiplicity fluctuations are strongly suppressed, making the chan- 
nel n -H- p + 7r~ dominant. The correlation is weaker elsewhere in the T-/j,b 
phase diagram. The final state p pw -, Fig. (8.5) (right), is similar to the primordial 
scenario for low temperature. For larger T and ris, resonance decay then again 
correlates the event-by-event multiplicities of p and n~ , when compared to the 
primordial scenario. 

Micro Canonical Ensemble 

The primordial MCE correlation coefficient p pn +, Fig. (8.6) (left), is generally 
weaker than its CE counterpart \p p n ^+\ < {pp^+l, with the exception of a deeper 
pocket in the bottom right corner of the phase diagram, denoting cold baryonic 
matter, where \p™°+\ > \p^+\. The anti-correlation p <H> n + ir + is made stronger 
at low temperature in the MCE by combined energy and charge conservation. 
The sum of the multiplicities of p and n cannot vary much due to strongly sup- 
pressed p and n yields. The anti-correlation between the multiplicities of p and ir + 
is strongly mediated by the n multiplicity. At larger temperature other channels, 
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like A + <H> n + 7r + , seem preferable for energy conservation. Charge conserva- 
tion does not distinguish between these two channels, leading to \p p e 1T +\ > \p™+\ 
elsewhere in the phase diagram, except in the pocket. 




MGeV] MGeV] 



Figure 8.6: MCE correlation coefficient between the multiplicities of protons and 
positively charged pions p p7r +. Both, primordial (left), and final state (right). 
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Figure 8.7: MCE correlation coefficient between the multiplicities of protons and 
negatively charged pions p pn -. Both, primordial (left), and final state (right). 

The final state MCE correlation coefficient p p7r +, Fig. (8.6) (right), is negative, 
p p7r + < 0, and notably stronger than in the final state CE. Energy conservation, 
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and hence anti-correlation between primordial parent particles and primordial p 
and 7r + , is responsible. 

The primordial MCE correlation coefficient p pn -, Fig. (8.7) (left), is again pos- 
itive, p pn - > 0, like in the CE. The channel n -h- p + 7r~ is not much affected by 
additional energy conservation, with the notable exception of a stronger positive 
correlation p^f > in cold baryonic matter. 

From a comparison to the final state MCE correlation coefficient p pn -, 
Fig. (8.7) (right), again a positive resonance decay contribution is visible. One 
notes that the CE enhancement in net-baryon free and hot matter has vanished 
in the MCE. Energy fluctuations in the CE at a given temperature are strongest 
if the system is charge neutral. In the MCE the energy content is fixed, and 
anti-correlation amongst primordial parent resonance and primordial p and 7r + 
balances the enhancement due to resonance decay. 




Figure 8.8: Primordial MCE correlation coefficient between the multiplicities of 
the doublely charged A ++ (1232) and positively charged pions Pa++(i232)tt+ (left), 
and between the multiplicities of the doublely charged A ++ (1600) and positively 
charged pions Pa++(1600)tt+ (right). 

To conclude the discussion of multiplicity correlations primordial MCE cor- 
relation coefficient between the multiplicities of the doublely charged A ++ (1232) 
and positively charged pions Pa++(i232)tt+, and between the multiplicities of the 
doublely charged A ++ (1600) and positively charged pions Pa++(1600)tt+> are shown 
in Fig. (8.8) (left) and (right) respectively. 
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The A resonances have larger rest-mass than the ground state nucleons n 
and p. Hence, they start to contribute to the system partition function only at 
larger T, and preferably at large ps- The correlation between the event-by-event 
multiplicities of A resonances and 7r mesons is the stronger the larger a fraction of 
the total energy density they make up. The maxima in Fig. (8.8) (left) and (right) 
are hence in different positions. The heavier particle shows a weaker correlation. 
The correlation coefficients Pa++(1232)tt+ an d PA++(i600)?r+ should be compared to 
the primordial MCE p pw + in Fig. (8.6) (left). The CE yields a similar structure. 
In the GCE no primordial correlations appear in the ideal gas approximation. 

Meson Baryon Correlations 

Lastly, the correlation coefficients pmb between the event-by-event multiplic- 
ities of mesons plus anti-mesons, Nm = Am + N^, and the event-by-event multi- 
plicities of baryon plus anti-baryons, N B = A B + iVg-, are considered in different 
ensembles. In Fig. (8.9) the temperature - baryon chemical potential phase dia- 
gram is shown for the final state correlation coefficient pmb m the GCE (left) 
and CE (right). In Fig. (8. 10) the MCE correlation coefficient Pmb is shown, both 
primordial (left) and final state (right). 

No primordial correlation amongst distinct groups of hadrons exists in the 
GCE. Resonance decay is then again the only source of correlation, Fig. (8.9) (left), 
between mesons and baryons. Yet, resonances essentially do not decay into baryon 
anti-baryon pairs, B + B. At low temperature, where few resonances are formed, 
the correlation coefficient is then pmb ~ 0. On the other hand, when the net- 
baryon density is large (positive or negative), either baryon or anti-baryons are 
strongly suppressed, and the multiplicities N M = N M + N^, are essentially cor- 
related with Ab (pb > 0) or Ag- (ps < 0). 

Also in the CE primordial correlations are absent, due to the particular choice 
made for strangeness and electric charge chemical potentials pq = ps = 0. Reso- 
nance decay enhancement, Fig. (8.9) (right), is again visible at large T and ps- At 
large T, but small pb, however the correlation coefficient is somewhat stronger 
(than in the GCE). The expectation value of the sum of baryons plus anti- 
baryons Ag = JVb + Ag- is largest (at a given temperature) for charge neutral 
matter. However, the strangeness and electric charge content of the baryonic sec- 
tor is still fluctuating, leading to a weak and positive correlation to the mesonic 
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sector. 




Figure 8.9: Final state correlation coefficient Pmb between the event-by-event 
multiplicities of mesons plus anti-mesons and baryon plus anti-baryons in the 
GCE (left) and CE (right). 




Figure 8.10: MCE correlation coefficient Pmb between the event-by-event mul- 
tiplicities of mesons plus anti-mesons and baryon plus anti-baryons in the MCE, 
both, primordial (left) and final state (right). 



In the primordial MCE, Fig. (8. 10) (left) one observes a strong correlation 
between the event-by-event multiplicities of mesons plus anti-mesons and baryon 
plus anti-baryons. Unlike in the GCE and CE, the primordial MCE correlation 
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coefficient is non-vanishing and negative, Pmb < 0. The difference between the 
multiplicities of baryons and of anti-baryons, i.e. the system baryon number 
is conserved, B = N B — N#. Yet, fluctuations of their sum, N B = N B + jV-g-, 
however release or bind a large amount of energy. But as energy is conserved 
in the MCE, the multiplicity Nm = Nm + needs to adjust. Energy is spent 
either on the production of a baryon anti-baryon pair, or on the production of a 
few (> 2) mesons and anti- mesons. This effect is particularly strong in hot and 
neutral matter. Resonance decay in the MCE, Fig. (8. 10) (right), leads again to 
a positive correlation large T and \ib amongst final state hadrons, yet, a strong 
anti-correlation at large T and \ib ~ remains. 

Following the (E)/(N) ~ 1 GeV freeze-out line of Fig. (8.10) (right), one would 
observe a weak positive correlation at small ^/snn, and a strong anti-correlation 
at large ^snn- In the SPS energy range a maximum might emerge. In particular, 
for the region of the phase diagram accessible to RHIC experiments, the MCE 
formulation implies that some events will be mostly composed of mesons and 
anti-mesons, while other events should be mostly composed of baryons and anti- 
baryons. The bulk of the events should be quite normal, but some events could be 
either predominantly mesonic (bosonic) or predominately baryonic (fermionic). 

8.3 Multiplicity Fluctuations 

In this section the phase diagrams for primordial and final state multiplicity 
fluctuations of positively and negatively charged hadrons u + and U- are com- 
pared in CE and MCE in Figs. (8. 11-8. 14). GCE calculations will, in addition, be 
presented in Section 9.1. 

Some general comments attempt to summarize. At positive /xg, also the 
electric charge density is positive. The electric net-charge content of the sys- 
tem Q = N + — A_ is conserved. As now the average multiplicity of positively 
charged hadrons is larger than the one of negatively charged hadrons, (N + ) > 
(A_) (but as the variances ought to be equal, ((AA + ) 2 ) = ((AA„) 2 )) one finds 
the multiplicity distribution of positively charged hadrons to be more narrow than 
the one of negatively charged hadrons, u + < w_. At negative baryon chemical 
potential /ig, i.e. for negative electric charge density, the situation is reversed, 
and u + > ui~ . Their fluctuations are equal for iib — 0. 
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Figure 8.11: CE fluctuations of the multiplicities of positively charged 
hadrons u + . Both, primordial (left), and final state (right). 




H [GeV] hb [GcV] 

Figure 8.12: MCE fluctuations of the multiplicities of positively charged 
hadrons uj + . Both, primordial (left), and final state (right). 



With increasing temperature the system becomes more and more relativistic, 
and particle number increasingly fluctuates in the MCE. In the CE this is already 
less apparent. A very mild temperature dependence of the scaled variance of mul- 
tiplicity fluctuations arises only due to heavy hadrons carrying multiple charges. 
In a CE with only one species of particles in Boltzmann approximation, the scaled 
variance in neutral matter is u + = uj~ = 0.5 [94] independent of temperature. In 
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the primordial Boltzmann GCE fluctuations are Poissonian independent of either 
thermal parameter. Note the different color scales in Figs. (8. 11-8. 14). 




Figure 8.13: CE fluctuations of the multiplicities of negatively charged 
hadrons u~ . Both, primordial (left), and final state (right). 




Figure 8.14: MCE fluctuations of the multiplicities of negatively charged 
hadrons u>~. Both, primordial (left), and final state (right). 

Comparing primordial CE Fig.(8.11) (left) primordial MCE Fig. (8. 12) (left), 
one finds at low temperature u;™ ce < u/L e . In a cold gas particle multiplicity 
of positively charged hadrons cannot fluctuate much due to combined energy 
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and charge conservation. At large T and fig, however very mildly, u;™ ce > a/j_ e . 
Baryon multiplicity fluctuations induce meson fluctuations in hot baryonic mat- 
ter, see Fig. (8. 10). The net-baryon number is conserved, B = Ng — Ng-. Yet, 
fluctuations in the sum of baryon number and anti-baryon number, Ng + N-g, 
release or bind a large fraction of the available energy of the system, with the 
meson multiplicity adjusting. In the CE, with only charge conservation, meson 
multiplicity fluctuations are less affected. 

A weak, yet generally useful, rule of thumb is: Resonance decay in the CE 
Fig. (8. 11) (right) increases u, while in the MCE Fig.(8.12) (right) decreases u, 
when compared to their respective primordial scenarios (left). With the notable 
exception of very large T. There it seems more economical to produce particles 
via decay of resonances, rather than thermally with large kinetic energy due to 
temperature. The correlation between primordial parents and stable daughters is 
then weaker. And u+ nal > u p " m in the MCE is possible. In final state one always 
finds u™ ce < It is mentioned in this context, that even in the ultra-relativistic 
limit u mce < u ce [95]. 

Having stated this rule of thumb, its caveat shall be discussed. Its validity 
depends strongly on possible cluster formation and their decay modes 3 . Yet 
energy can only be spent once. Either on primordial particle production, or the 
formation of a cluster, bag, string, or Hagedorn state. If particle production 
via cluster is more economical, i.e. more energy spent "on mass" and less "on 
motion" , then final state fluctuations can be larger than primordial fluctuations 
in the MCE. 

Now, turning to the fluctuations of negatively charged hadrons in the posi- 
tive fig half of the T-fi B phase diagram, the diagrams depicted in Figs. (8. 11, 8. 12) 
could be connected to the reflected Figs. (8. 13,8. 14). The scaled variances of pos- 
itively and negatively charged hadrons are equal at fig = 0. For a change in sign 
of fiB, the roles of u + and U- are interchanged, i.e. the phase diagrams are not 
symmetric around fig — 0. 

The scaled variance cu_ is maximal where the ratio (N_) / (N + ) is minimal, 
i.e. when the electric charge density is strongly positive. For the primordial 
case, comparing CE, Fig.(8.13) (left) to MCE, Fig.(8.14) (left), one finds gen- 
erally uj™ C£ < wf_ e . The final state enhancement at low fig and large T in 

3 In this context the terms cluster, string, Hagedorn state, or bag are used equivalcntly. 
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the CE, Fig. (8. 13) (right) is due to resonance decay. For the final state MCE, 
Fig. (8. 14) (right) this enhancement is weaker, due to correlation of primordial 
parent resonances with stable particles via energy conservation. Similar to u + , at 
very large T and \ib one finds o;™ ce > cuf 62 in final state. Again, baryon multiplicity 
fluctuations induce meson multiplicity fluctuations. 

8.4 Discussion 

The T~hb phase diagram of the hadron resonance gas model has been explored 
in this chapter for fluctuation and correlation observables in their dependence on 
the chemical freeze-out parameters temperature and baryon chemical potential, 
omitting limited acceptance effects, but still considering resonance decay. 

Grand canonical baryon number, strangeness and electric charge fluctuations 
and correlations are strongly sensitive to the degrees of freedom, i.e. different 
particle species, available to the system. The contributions of different hadron 
species is rather different in different corners of the T-\i% phase diagram. Mesons 
dominate at low temperatures. With increasing temperature and net-baryon 
density however baryonic degrees of freedom take a larger fraction of the system 
entropy. Two charges are strongly (anti-)correlated if particles which carry both 
quantum numbers are abundant. The same argument was essentially used to 
explain the momentum space dependence of charge fluctuations and correlations 
in Chapter 4. 

Multiplicity correlations in the canonical and micro canonical ensembles follow 
suit. So is for instance the correlation between protons and charged pions strongly 
mediated by the neutron in cold baryonic matter. These correlations appear 
solely due to implementation of global conservation laws for charges and energy, 
and not due to a local interaction amongst constituent particles. With increasing 
temperature more channels, i.e. resonance formation and decay, become available 
to the statistical system, and subtle differences between canonical and micro 
canonical ensembles emerge. 

Multiplicity fluctuations highlight then the connection between energy spent 
on resonance formation and energy spent on primordial particle production. The 
event-by-event multiplicities of stable particles and unstable parent resonances are 
anti-correlated in the micro canonical ensemble, while being almost unaffected in 
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the canonical ensemble, and independent of each other in the grand canonical 
ensemble. 

Above the freeze-out line, no hadronic phase should exist. Even in the presence 
of critical phenomena close to the phase transition, baseline contributions of the 
kind discussed here, and in previous chapters, should remain. 
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Chapter 9 

The Chemical Freeze-out Line 



In this chapter the results of the hadron resonance gas for the scaled variances 
and correlation coefficients are presented for freeze-out parameters following the 
constant average energy per average particle multiplicity (E) / (N) ~ 1 GeV chem- 
ical freeze-out line [20]. Different ensemble predictions will be confronted with 
recent NA49 measurements of charged hadron multiplicity fluctuations. 

Mean hadron multiplicities in central heavy ion collisions at high energies 
are usually fitted within the GCE formulation of the hadron resonance gas 
model. The fit parameters then are: volume V, temperature T, baryon chem- 
ical potential and saturation parameters 7*, [116] and 7 g [117, 118]. The 
former allows for non-equilibrium strange hadron yields, while the latter de- 
scribes non-equilibrium light quark numbers. For reviews of the model see 
Refs. [3, 4, 5, 6, 7, 8, 9, 10]. 

Recent discussion of system size and energy dependence of freeze-out param- 
eters and comparison of different (chemical) freeze-out criteria can be found in 
Refs. [18, 19, 20, 21, 22]. The set of model parameters, V,T,/j,b, and j s , chosen 
in the following should correspond to chemical freeze-out conditions in heavy ion 
collisions 1 . The evolution of model parameters with collision energy is taken from 
previous analysis of mean hadron multiplicity data. Hadron resonance gas model 
fits should then be able to provide a link between experimental control parame- 
ters, center of mass energy, and size of colliding ions, and the region probed in 
the phase diagram of strongly interacting matter. 

1 Thc Bose-Einstein condensation model is not considered, as the analytic expansion used for 
the calculation here is inappropriate. However, one notes that energy conservation essentially 
prohibits diverging susceptibilities, if degrees of freedom have finite mass. 
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The dependence of the baryon chem- 
ical potential ji B on the collision energy 
is parametrized as [18]: \ib (^/sjvjv) = 
1.308 GeV ■ (1 + 0.273 y/s^)' 1 , where 
the center of mass energy per nucleon 
pair, y/SNN, is taken in units of GeV. The 
system is assumed to be net-strangeness 
free, S = 0, and to have the elec- 
tric charge to baryon density ratio of 
the initial colliding nuclei, Q/B = 0.4. 
These two conditions define the system's 
strangeness, /is, an d electric charge, /xq, 
chemical potentials. The average energy 
per average particle yield [20] (E) / (N) = 
1 GeV was chosen to fix the tempera- 
ture T. Finally, the strangeness satura- 
tion factor is parametrized [21], 7 S = 1 — 
0.396 exp(— 1.23 T//i B ). The asymptotic 
(or large volume) solutions for the scaled 
variance and correlation coefficients 
sumed. This determines all parameters of 
the model. 

Here, some further details should be mentioned. Quantum statistics effects 
are considered, while the non-zero (Breit-Wigner) width of resonances is disre- 
garded. The standard THERMUS particle table includes all known particles and 
resonances up to a mass of about 2.5 GeV and their respective decay channels. 
Heavy resonances do not always have well established decay channels. Branch- 
ing ratios given in THERMUS are re-scaled to unity, where it was necessary 
to ensure global charge conservation. To make correspondence with NA49 data, 
both strong and electromagnetic decays should be taken into account, while weak 
decays should be omitted. 

In Section 9.1 multiplicity fluctuations and correlations are considered for 
different ensembles along the chemical freeze-out line. In Section 9.2 a comparison 
to recent NA49 data on charged hadron multiplicity fluctuations is attempted. 




0.2 0.3 0.4 0.5 

H B (GeV) 

Figure 9.1: Chemical freeze-out 
points in the T-fiB plane for central 
Pb+Pb collisions. The solid line in- 
dicates the (E)/(N) = 1 GeV freeze- 
out line. Square markers denote SPS 
beam energies from 20A GeV (right) 
to 158 A GeV (left). Full and open 
circles are the best fit parameters 
from Refs. [21] and [22], respectively. 
The Figure is taken from [103]. 
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9.1 Fluctuations and Correlations 

In Fig. (9.2) the scaled variances of negatively charged hadrons, u~ (left), and 
the scaled variances of positively charged hadrons, u + (right), are shown, both 
primordial and final state, along the chemical freeze-out line for most central 
Pb+Pb (Au+Au) collisions. In Fig. (9.3) (left) the scaled variances of all charged 
particles, u ch , is shown 2 . The three standard ensembles are considered in turn, 
and some qualitative features are discussed. 
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Figure 9.2: The scaled variances of negatively charged particles, u~ (left), and 
the scaled variances of positively charged particles, u + (right), both primordial 
and final state, along the chemical freeze-out line for central Pb+Pb (Au+Au) 
collisions. Different lines present GCE, CE, and MCE results. Symbols on the 
lines for final state fluctuations correspond to specific collision energies, ranging 
from SIS and AGS (triangles), SPS (squares), RHIC (circles), to LHC (diamonds). 
Arrows indicate the effect of resonance decay. Figures are taken from [103]. 



The effects of Bose-Einstein and Fermi-Dirac statistics are seen best in the 
primordial GCE. At low temperatures most positively charged hadrons are pro- 
tons, and Fermi-Dirac statistics dominate, leading to sub-Poissonian multiplicity 
fluctuations, 0J+ ce , u c ^ e < 1. On the other hand, in the limit of high tempera- 
ture most charged hadrons are pions and the effects of Bose-Einstein statistics 
dominate. The multiplicity distribution is then wider than a Poissonian refer- 

2 Figures are taken form [103]. In previous chapters the scaled variances of positively, nega- 
tively, and all charged hadrons were denoted as uj + , cj_ , and oj± . In this chapter these quantities 
are refered to as u + , u~ , and uj ch . 
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Figure 9.3: The scaled variances of all charged particles, u ch , both primor- 
dial and final state (left), and the final state correlation coefficient px-n (right) 
along the chemical freeze-out line. HSD transport model results for pKn 
in p+p and Au+Au collisions are shown for comparison. Figures are taken 
from [103] (left) and [111] (right). 




Figure 9.4: The final state correlation coefficient px P (left), and the final state 
correlation coefficient p p7T (right) along the chemical freeze-out line. HSD trans- 
port model results for Au+Au collisions are shown for comparison. Figures are 
taken from [142]. 



ence, oof ce , io C g Ce > 1. Along the chemical freeze-out line, u gce is always slightly 
larger than one, as n~ mesons are the most abundant negatively charged hadrons 
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at, both, low and high temperatures, Fig. (9.2) (left). 

The bump in cu+ e for final state particles seen in Fig. (9.2) (right) at small colli- 
sion energies is due to correlated production of protons and ir + mesons from A ++ 
decay. This single resonance contribution dominates for u^ ce at small collision en- 
ergies (small temperatures), but becomes relatively unimportant at high collision 
energies. Along the freeze-out line the final state GCE correlation coefficient p p7r +, 
Fig. (8.3), has a maximum in the same region. For the final state id~ ce no bump is 
seen due to a weaker resonance decay contribution at low ^/snn- For comparison, 
the GCE resonance decay enhancement of multiplicity fluctuations is strongest 
for all charged hadrons, Fig. (9.3) (left). 

A minimum in u~ e for final particles is seen in Fig. (9.2) (left). This is 
due to two effects. As the number of negatively charged particles is relatively 
small, (iV_) <C (N + ), at low collision energies, both CE suppression and the ef- 
fect of resonance decay are small. With increasing y/s^^, the CE effect alone 
leads to a decrease of u~ e , but resonance decay leads to an increase of u~ e . The 
combination of these two effects, CE suppression and resonance decay enhance- 
ment, leads to a minimum of u~ e . The primordial u~ e decreases monotonously 
with v /J/vtv, compare also Fig. (8. 13). At large ^/snn the system is essentially 
charge neutral, CE suppression of w+ and enhancement of u~ e therefore cease, 
and one finds u~ e = u+ e accordingly. Both, primordial and final state, u+ e in 
Fig. (9.2) (right) hence increase with yj snn, compare Fig. (8. 11). For a charge 
neutral system the CE effects cancel exactly for the multiplicity fluctuations of all 
charged hadrons, both in final state and primordial, = w^ e , Fig. (9.3) (left). 
Compare also Fig. (6. 13). 

As to be expected, oo mce < u ce , along the chemical freeze-out line. Energy 
conservation further suppresses particle multiplicity fluctuations. A particular 
feature of the MCE is the additional suppression of multiplicity fluctuations after 
resonance decay, compare Figs. (8. 12, 8. 14). Energy conservation does not allow for 
strong multiplicity fluctuations at low temperature. With increasing s^n one 
finds oo^ ce and cu+ ce to grow. On the other hand, w~ ce drops. Negatively charged 
particles are outnumbered by positively charged particles at small collision ener- 
gies. Therefore, as the variances of the distributions of positively and negatively 
charged hadrons are the same, one finds u;~ ce t° grow towards smaller ^snn, es- 
sentially due to a dropping mean value (N~). One also notes, that the resonance 
decay enhancement of u~ e around y^sJvjv = 10 GeV is not seen for w~ ce , due 
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to a strong correlation of primordial parent resonances with primordial hadrons 
resulting from global energy conservation. This effect is even stronger for u^ ce . 

The correlation coefficient p^ n between the event-by-event multiplicities 
of charged kaons, N K + + N K ~, and charged pions, N n + + N n -, is shown in 
Fig. (9.3) (right) as a function of center of mass energy ^/snn for final state 
hadrons only. Fig. (9.4) additionally shows the correlation coefficient pK P (left) 
between the event-by-event multiplicities of charged kaons and protons plus anti- 
protons, iVp + Np (left), and the correlation coefficient p pn (right) between the 
event-by-event multiplicities of protons plus anti-protons and both charged pi- 
ons 3 . The three correlation functions are discussed in turn. 

In the final state GCE the correlation coefficient px-n remains weak. A very 
mild increase at large ^Jsnn is due to rare resonances decay feeding this chan- 
nel. In the final state CE p^ is somewhat weaker than in the GCE. Charge 
conservation effects remain small for the particle selections of both charged 
kaons, N K + + N K -, and both charged pions, N n + + N w -. At large ^/snn, i- e - 
for a neutral system, p 9 ^ ~ p c ^ n , CE suppression ceases entirely. As to be ex- 
pected, the event-by-event multiplicities of charged kaons and charged pions, are 
anti-correlated by energy conservation in the MCE. Due to a missing strong decay 
contribution, the correlation coefficient p 1 ^ stays negative at all y/s^N- 

The final state GCE correlation coefficient px P is also rather weak due to a 
weak resonance decay contribution. In the CE the event-by-event multiplicities 
of protons plus anti-protons and charged kaons are anti-correlated at low ^/Sjvjv- 
Their anti-correlation, p^p < 0, is mediated by correlation with primordial 
strangeness carrying baryons, in particular the A, at large baryon chemical po- 
tential fiB- With increasing A /J/v r /v the temperature rises and ri B drops. Charge 
conservation effects cease again for a neutral system, p 9 ^ ~ p^ p , for the particle 
selections of both charged kaons, N K + + N K ~ , and both charged pions, N n ++N n - . 
At small a/saT/v the results of CE and MCE are similar. At larger ^/saT/v, however, 
the correlation p^ p e remains negative, again due to energy conservation. 

Lastly, the final state correlation coefficient p pn is considered along the chem- 
ical freeze-out line in different ensembles. In Section 8.2 the correlation coef- 

3 Figures are taken form [141, 142]. In previous chapters the correlation coefficients denoted 
only the correlation between proton and either charged pion. In this chapter always both, 
proton and anti-proton, or both charged kaons, or both charged pions are correlated. 
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ficient p pn + and p p7r - were discussed. Here, both proton and anti-proton and 
both charged pions are counted, while in Section 8.2 only the proton was cor- 
related with either of the charged pions. In the final state GCE the correlation 
amongst not charge separated multiplicities is stronger than for charge separated 
multiplicities, > p^+,p^- 

The anti-proton multiplicity is strongly suppressed at low ^/snn, and hence 
their contribution to p p7r is weak. On the other hand, compare Fig. (8.3), one 
finds a strong resonance decay contribution at low ^snn- Approaching a neu- 
tral system, one finds all three final state GCE correlation coefficients to drop. 
Charge conservation, i.e. CE, effects again cease in a neutral system for the 
fluctuations of multiplicities iVp + Np and N w + + N n -. At low ^/snn a weak, 
but positive, final state correlation remains, compare Figs. (8.4,8.5). In the MCE, 
one finds a baryon density induced positive correlation. Anti-protons are essen- 
tially absent. The correlation coefficient p p7r + is negative, while the correlation 
coefficient p pvr - is positive. The correlation coefficient p p7r is then weakly pos- 
itive due to abundant neutrons and A resonances, allwoing for n -H- p + 7r~ 
and A -h- p + it. At large -x/snn, hence for a neutral system, the anti-correlation 
between primordial parent particles and primordial protons, anti-protons and 
charge pions takes over. For the charge separated counterpart, one finds on the 
other hand p p7r - ~ p p7r + ~ 0, compare Figs. (8. 6, 8. 7), for large T and p# ~ 0. 

The results of the relativistic transport model HSD on correlations coefficients 
between particle multiplicities of charged kaons, charged pions, and protons plus 
anti-protons are summarized in Refs. [141, 142]. 

9.2 Comparison with NA49 Data 

The fluctuations in nucleus-nucleus collisions are studied on an event-by-event 
basis: a given quantity is measured for each collision and a distribution of this 
quantity is constructed for a selected sample of these collisions. In the following 
a few aspects are discussed. Details of the experimental setup and analysis can 
be found in [58, 59]. 
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Centrality Selection 

The fluctuations in the number of nucleon participants provide a dominant 
contribution to hadron multiplicity fluctuations. In the language of statistical 
models, fluctuations of the number of nucleon participants correspond to volume 
fluctuations caused by the variations in the collision geometry. Mean hadron 
multiplicities are proportional (in the large volume limit) to the volume, hence, 
volume fluctuations translate directly to the multiplicity fluctuations. Thus, a 
comparison between data and predictions of statistical models should be per- 
formed for results which correspond to collisions with a fixed number of nucleon 
participants. 

Due to experimental limitations it is only possible to measure the number 
of participants of the projectile nucleus, Np OJ , in fixed target experiments (e.g. 
NA49 at the CERN SPS). This is done by measuring the energy deposited in a 
downstream Veto calorimeter. A large fraction of this energy is due to projectile 
spectators Ng ro3 . Using baryon number conservation for the projectile nucleus 
(A = N p ro:> + A r | roj ) the number of projectile participants can be estimated. 
However, also a fraction of non-spectator particles, mostly protons and neutrons, 
contribute to the Veto energy [58, 59]. Furthermore, the total number of nucleons 
participating in the collision can fluctuate considerably even for collisions with 
a fixed number of projectile participants, especially for peripheral collisions (see 
Ref. [156, 157]). This is due to fluctuations of the number of target participants. 
The consequences of the asymmetry in the event selection depend on the dynamics 
of nucleus-nucleus collisions (see Ref. [158] for details). Still, for the most central 
Pb+Pb collisions selected by the number of projectile participants an increase of 
the scaled variance can be estimated to be smaller than a few % [158] due to the 
target participant fluctuations. In the following hadron resonance gas results will 
be compared with the data of the NA49 collaboration on the 1% most central 
Pb+Pb collisions at 20A-158A GeV [58, 59]. 

Modeling of Acceptance 

In the experimental study of nuclear collisions at high energies only a fraction 
of all produced particles is registered. Thus, the multiplicity distribution of the 
measured particles is expected to be different from the distribution of all produced 
particles. In general, in statistical models, the correlations in momentum space 
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are caused by resonance decays, quantum statistics and the energy-momentum 
conservation law. In this section these correlations are neglected and the proce- 
dure of Appendix E is applied for different ensembles. This may be reasonably 
valid for u + and u~ , as most decay channels only contain one positively (or neg- 
atively) charged particle, but is certainly much worse for u) ch , for instance due to 
decays of neutral resonances into two charged particles. In order to limit corre- 
lations caused by resonance decays, the analysis is restricted to negatively and 
positively charged hadrons. 

The NA49 acceptance used for the fluctuation measurements is located in 
the forward hemisphere (1 < y(n) < ybeam, where y(ir) is the hadron ra- 
pidity calculated assuming pion mass, and shifted to the collision center of 
mass system [58, 59]). The acceptance probabilities for positively and nega- 
tively charged hadrons are approximately equal, q + ~ q~, and the numerical 
values at different SPS energies are: q ± = 0.038, 0.063, 0.085, 0.131, 0.163 
at v/i^F = 6.27, 7.62, 8.77, 12.3, 17.3 GeV, respectively. 

Comparison with NA49 Data 

Fig. (9.5) presents the scaled variances u~ and u + , obtained from applying the 
acceptance procedure of Appendix E to hadron resonance gas calculations in the 
GCE, CE, and MCE shown in Fig. (9.2). 
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Figure 9.5: The scaled variances for negatively charged (left) and positively 
charged hadrons (right) along the chemical freeze-out line for central Pb+Pb 
collisions at SPS energies. Full squares show the data of NA49 [58, 59]. Total 
(statistical+systematic) errors are indicated. Lines show the GCE, CE, and MCE 
acceptance scaling estimates. Figures are taken from [103]. 
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From Fig. (9.5) it follows that the NA49 data for u~ and u + extracted from 
1% of the most central Pb+Pb collisions at all SPS energies are best described 
by the results of the hadron resonance gas model calculated within the MCE. 
The data reveal even stronger suppression of the particle number fluctuations. 
Given that the NA49 acceptance is located in the forward rapidity hemisphere, 
one would expect further suppression of multiplicity fluctuations compared to the 
acceptance scaling estimate employed here. See also Chapters 6 and 7. 

As discussed, the multiplicity distribution in statistical models in the full 
phase space and in the large volume limit approaches a normal distribution. 
Particle detection is modeled by the simple procedure, Appendix E, which is 
applicable to any form of full acceptance distribution P^N). In the following 
multiplicity distributions in limited acceptance are discussed, and the statistical 
model results in different ensembles are compared with data on negatively and 
positively charged hadrons. The multiplicity distribution is, after the acceptance 
scaling procedure was applied, not Gaussian anymore. It is enough to mention 
that a Gaussian is symmetric around its mean value, while the new distribution, 
in the limit of very small acceptance, converges to a Poissonian. 

In the following for each beam energy the volume is adjusted to fit the mean 
values for negatively (V~) and positively (V + ) charged yields, separately. Note 
that values for the volume are about 10 — 20% larger than the ones in [21, 22], 
which were obtained using a much less stringent centrality selection (here only 
the 1% most central data is analyzed). One finds that the volume parameters V~ 
and V + deviate from each other by less than 10%. Deviations of a similar mag- 
nitude are observed between the data on hadron yield systematics and hadron 
resonance gas model fits. 

In order to allow for a detailed comparison of the distributions the ratios of 
data and model distributions to a Poissonian reference are presented in Fig. (9.6). 
The results for negatively and positively charged hadrons at 20 A GeV, 30A GeV, 
40^4 GeV, 80v4 GeV, and 158A GeV are shown separately. The convex shape of 
the data reflects the fact that the measured distribution is significantly narrower 
than a Poissonian. This suppression of fluctuations is observed for both charges, 
at all five SPS energies, and it is consistent with the results for the scaled variance 
shown and discussed previously. The GCE hadron resonance gas results are 
broader than the corresponding Poisson distribution. The ratio has a concave 
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Figure 9.6: The ratio of the multiplicity distributions to Poisson ones for neg- 
atively charged (top) and positively charged hadrons (bottom) produced in cen- 
tral (1%) Pb+Pb collisions at 20A GeV, 30A GeV, AOA GeV, 80A GeV, and 
158A GeV (left) to (right) in the NA49 acceptance. The experimental data (solid 
points) of NA49 [58, 59] are compared with the acceptance scaling estimates of 
the hadron resonance gas model obtained within different statistical ensembles, 
the GCE (dotted lines), the CE (dashed-dotted lines), and the MCE (solid lines). 
Figures are taken from [103]. 

shape. Introduction of quantum number conservation laws leads to the convex 
shape and significantly improves agreement with the data. Further improvement 
is obtained by additional introduction of energy conservation. 



9.3 Discussion 

In this chapter multiplicity fluctuations and correlations have been analyzed 
for thermal parameter sets following the chemical freeze-out line for most central 
heavy ion collisions. The evolution of model parameters was taken from previ- 
ous hadron resonance gas model comparison to experimental on average hadron 
production yields. These fits provide a connection between experimental con- 
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trol parameters center of mass energy and ion size, and the region of the phase 
diagram probed by the experiment. 

Taking the model in its present form, a naive comparison the NA49 data on 
charged hadron multiplicity fluctuations has been attempted. High resolution 
of the NA49 experimental data could in principle allow to distinguish between 
multiplicity fluctuations expected in hadron resonance gas model for different 
statistical ensembles. The measured spectra clearly favor predictions of the micro 
canonical ensemble. Worse description is obtained for the canonical ensemble and 
disagreement is seen considering the grand canonical one. 

All calculations were performed in the thermodynamical limit which is a 
proper approximation for the considered reactions. Thus, these results could be 
treated as an observation of the recently predicted [93, 94, 95, 99, 100, 115, 159, 
160, 161] suppression of multiplicity fluctuations due to conservation laws in rel- 
ativistic gases in the large volume limit. Similar observations have been made in 
detector physics [107], and in low temperature Bose-Einstein condensation [109]. 

The experimental resolution of the NA49/61 detector for the measurement 
of enhanced fluctuations due to the onset of deconfinement can be increased 
by increasing acceptance. This will give a chance to observe, for example, the 
dynamical fluctuations discussed in Ref. [41, 42]. The observation of the MCE 
suppression effects of the multiplicity fluctuations by NA49 was possible only 
because a selection of a sample of collisions without projectile spectators. This 
selection seems to be possible only in the fixed target experiments. 

In collider experiments nuclear fragments which follow the beam direction 
cannot be measured, and another measure has to be used to determine which 
centrality class a given event (or collision) would fall into. Often, a charged 
hadron multiplicity reference distribution is employed. This measure (like any 
other) provides, however, only indirect access to the initial collision geometry. 

Likewise, the energy deposited in the interaction region, and hence the amount 
of energy available for particle production cannot be considered as fixed across a 
whole sample of heavy ion collisions. Yet, for each individual collision the energy 
content available for particle production and (collective) motion is fixed. Thus, 
even if a centrality class of collision events cannot be seen as a representative a 
micro canonical sample with the same values of global observables, the final state 
of each collision will in general fulfill energy-momentum and charge conservation. 
One could then hope that the transverse momentum and rapidity dependence 
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will not be washed out by a less stringent centrality selection, but merely shifted 
up or down. 

More differential data on multiplicity fluctuations and correlations are re- 
quired for further tests of the validity of the statistical models and observation of 
possible signals of the phase transitions. Simultaneous differential measurement 
of mean values (N + ), and (iV_), variances ((AiV + ) 2 ), and ((AiV_) 2 ), and of the 
covariance (AN + AN_) should lead to an improved understanding. 

The momentum space dependence for scaled variances demanded by the micro 
canonical ensemble is qualitatively also seen in (non-equilibrium) transport sim- 
ulations [131] and recent data of the NA49 collaboration [59]. The experimental 
and theoretical study of multiplicity fluctuations should be complemented with 
a study of multiplicity correlation coefficients. This study could in principle be 
extended by a transverse momentum measurement. This would the be either a 
three by three covariance matrix, or if Pt is charge separated, a four by four co- 
variance matrix. The ^Jsnn and ion size program [71] should be complemented 
by an acceptance effect study. Critical and quantum phenomena might more 
strongly feature at low transverse momentum. 
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Chapter 10 
Summary 

In this thesis the statistical properties of ideal relativistic hadronic equilibrium 
ensembles have been studied. Apart from the three standard canonical ensembles 
being discussed, a class of ensembles with finite thermodynamic bath has been 
introduced. The influence of global conservation laws, resonance decay, quantum 
statistics, finite acceptance in momentum space, and thermal freeze-out parame- 
ters on fluctuation and correlation observables of a sample of hadron resonance 
gas events have been considered. It was argued that the emerging picture would 
at least qualitatively apply to heavy ion physics. 

In Chapter 2 grand canonical joint distributions of extensive quantities were 
obtained by Fourier integration of the grand canonical partition function. An 
analytical expansion method for calculation of distributions at finite volume for 
the canonical as well as the micro canonical ensembles of the ideal relativistic 
hadron resonance gas was presented. The introduction of temperature into the 
micro canonical partition function and chemical potentials into the canonical 
partition function have lead to the identification of the grand canonical partition 
function with the characteristic function of associated joint probability distri- 
butions. Micro canonical and canonical multiplicity distributions could then be 
defined through conditional distributions, i.e. the probability of finding a certain 
multiplicity, while other parameters (global charge or energy) were taken to be 
fixed. 

In Chapter 3 joint distributions of extensive quantities were then considered 
for statistical systems with finite rather than infinite thermodynamic bath. To 
introduce the scheme, a micro canonical system was conceptually divided into two 
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subsystems. These subsystems were assumed to be in equilibrium with each other, 
and subject to the constraints of joint energy-momentum and charge conserva- 
tion. Particles are only measured in one subsystem, while the second subsystem 
provides a thermodynamic bath. By keeping the size of the first subsystem fixed, 
while varying the size of the second, one can thus study the dependence of statis- 
tical properties of an ensemble on the fraction of the system observed, i.e. assess 
their sensitivity to globally applied conservation laws. The ensembles generated 
are thermodynamically equivalent in the sense that mean values in the observed 
subsystem remain unchanged when the size of the bath is varied, provided the 
combined system is sufficiently large. As the three standard canonical ensembles 
remain particular idealizations of physical systems, these intermediate ensembles 
might be of phenomeno logical interest, too. These two chapters form the mathe- 
matical basis for analytical calculations and Monte Carlo simulations performed 
in this thesis. 

The analysis of hadron resonance gas events started with the study of the 
grand canonical ensemble in Chapter 4. The grand canonical ensemble is consid- 
ered to be the most accessible amongst the standard canonical ensembles. Due 
to the assumption of an infinite thermodynamic bath, occupation numbers of in- 
dividual momentum levels are un-correlated with each other. Likewise, particle 
multiplicities of any two distinct groups of particles appear un-correlated. To- 
gether with occupation number fluctuations unconstrained, by global constraints, 
all extensive quantities, except for the volume, fluctuate on an event-by-event (or 
micro state-by-micro state) basis. Therefore, one finds the energy content and 
particle multiplicity to be strongly correlated, while the average energy per par- 
ticle is un-correlated with multiplicity. The net-charge content of baryon number 
and strangeness in the system are correlated because some hadron species carry 
both charges. Different hadron species have different quantum numbers and fol- 
low different momentum spectra. The correlation between baryon number and 
strangeness then depends, like the one of energy and momentum, or energy and 
particle multiplicity, on which part of the momentum spectrum is accessible to 
measurement. 

In Chapter 5 joint distributions of extensive quantities were extrapolated to 
their micro canonical limit. For this iteratively samples of events were generated 
and analyzed for decreasing size of the thermodynamic bath. By construction, 
the distribution of extensive quantities considered for re-weighting converges to a 
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5- function, while the positions of the mean values stayed constant. Whereas mean 
transverse momentum per particle and particle multiplicity were un-correlated in 
the grand canonical ensemble, this is not true anymore for systems with finite 
heat bath. Through successive focusing on events in the vicinity of a chosen 
equilibrium value the effects of global conservation laws became apparent. This 
proceeded in a systematic manner, allowing for extrapolation of observables from 
their grand canonical to their micro canonical limit. A caveat is that the statistical 
uncertainty, associated with finite samples, grows as the sample-reject limit is 
approached. 

In Chapter 6 multiplicity fluctuations and correlations were studied for a neu- 
tral hadron resonance gas in limited acceptance in momentum space. The effects 
of resonance decay were considered, and the extrapolation scheme was applied to 
obtain the canonical and micro canonical ensemble limits. Comparison to analyt- 
ical asymptotic solutions, available for primordial distributions in limited accep- 
tance, suggests good agreement. The larger the number of conserved quantities, 
the larger the statistical uncertainty associated with finite samples of events. Yet, 
micro canonical effects are accurately reproduced by the Monte Carlo approach. 
As conservation laws are turned on, multiplicity fluctuations and correlations are 
modified. Momentum space effects on primordial multiplicity fluctuations and 
correlations arise due to conservation laws. For an ideal primordial grand canon- 
ical ensemble in the Boltzmann approximation (the original sample), multiplicity 
distributions are just uncorrelated Poissonians, regardless of the acceptance cuts 
applied, as particles are assumed to be produced independently. The requirement 
of energy-momentum and charge conservation leads to suppressed fluctuations 
and enhanced correlations between the multiplicities of two distinct groups of 
particles at the 'high momentum' end of the momentum spectrum compared to 
the 'low momentum' end of the momentum spectrum, provided some fraction of 
an isolated system is observed. Resonance decay does not change these trends. 

Chapter 7 was dedicated to the micro canonical ensemble. A simplified phys- 
ical system was chosen to allow for smoother discussion. Owing to available 
analytical solutions, Fermi-Dirac and Bose-Einstein effects have been included 
into the analysis. Bose-Einstein enhancement and Fermi-Dirac suppression of 
multiplicity fluctuations are strong in momentum space segments where occupa- 
tion numbers are large. This effect has been found to be much stronger than 
previous calculations of fully phase space integrated multiplicity fluctuations of 
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Chapter 10: Summary 



Fermi-Dirac and Bose-Einstein systems suggested. For systems in collective mo- 
tion is was found that the role of exactly imposed kinematic conservation laws 
is particularly important. Fluctuation and correlation observables transform un- 
der boosts, provided momentum conservation along the boost direction is taken 
into account. Lastly, it has been found that even in the thermodynamic limit 
long range correlations between disconnected regions in momentum space prevail. 
Multiplicities in different bins of rapidity, transverse momentum, or of azimuth, 
can have a non- vanishing correlation coefficient. 

In Chapter 8 the temperature - baryon chemical potential phase diagram of 
the hadron resonance gas model has been explored. Grand canonical charge cor- 
relations and fluctuations are different in different corners of the phase diagram. 
Like in the acceptance analysis, the correlation between two charges is strong 
when particles carrying both charges are abundant. At low temperature and 
baryon chemical potential mesons dominate, at high temperature and baryon 
chemical potential baryons take over. Fluctuations and correlations of charges 
systematically evolve accordingly. Multiplicity fluctuations and correlations in 
canonical and micro canonical ensembles follow suit. The influence of a particu- 
lar conservation law on multiplicity fluctuations is strong if the analyzed species 
is abundant. A comparison between ensembles with and without energy, charge, 
or momentum conservation reveals subtle differences, such as for instance reso- 
nance decay in canonical and micro canonical ensembles. Resulting (primordial) 
multiplicity correlations are not due to local interactions amongst constituents, 
but due to globally implemented conservation laws for energy and charge. 

In Chapter 9 multiplicity fluctuations and correlations were analyzed for ther- 
mal parameter sets following the chemical freeze-out line. Model parameters were 
taken from previous hadron resonance gas model comparisons to experimental on 
average hadron production yields. These fits provide a connection between ex- 
perimental control parameters center of mass energy and ion size and the region 
of the phase diagram probed by an experiment. A first comparison to available 
experimental data suggests good agreement with hadrons resonance gas calcu- 
lations. In particular the micro canonical formulation of the model seems to 
accurately reproduce qualitative features. 

The calculations presented here provide qualitative effects affecting multi- 
plicity fluctuations and correlations. These effects arise solely from statistical 
mechanics and conservation laws. It will be interesting to see whether these qual- 
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itative effects are visible in further experimental measurements of the momentum 
dependence of fluctuations and correlations. If so, these effects might well be 
of similar magnitude to the signals for new physics. Disentangling them from 
dynamical correlations will then be an important, and likely non-trivial task. 
Equilibrium correlations are residual correlations which remain after the system 
has 'forgotten' its history. If the system breaks up before equilibrium is attained, 
then correlations due to the initial configuration will remain. Yet, even if the 
system stays far from any equilibrium point throughout its lifetime, correlations 
due to global conservation laws remain. 
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Partition Function 



Appendix A 



Partition Function 



This section serves to provide a connection between Eqs. (2.3-2.5) and Eq.(2.13); 
namely to prove the following relation: 

Z^&iV^^m) = e~ p ^ Z (V,P»,Q j ) , (A.l) 

where Z(V, P M , Q 3 ) is the standard MCE partition function for a system of vol- 
ume V, collective four- momentum P M and a set of conserved Abelian charges Q J , 
as worked out in [113, 114]. The MCE partition function Z(V, P^jQ 3 ) counts 
the number of micro states consistent with this set of fixed extensive quantities. 
Likewise, one could interpret the partition function Z pil,Q] (V, (3, u^, fij) as the 
number of micro states with the same set of extensive quantities for a GCE with 
local inverse temperature /3, four-velocity w M , and chemical potentials fij. 
The starting point for this calculation is Eq.(2.13): 



7T OO 



(Zir) J {2,7V ) 



x exp 



(A.2) 



Taking a closer look at the exponential of Eq.(A.2), one can spell out the single 
partition functions, Eq.(2.9), of the cumulant generating function, Eq.(2.8), after 
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having used the substitutions Eqs.(2.11) and (2.12): 
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(A.3) 



Expanding the logarithm yields: 
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Replacing now the momentum integration in Eq.(A.4) by the usual summation 
over individual momentum levels t-^tt f d 3 p — ^ Y^. gives: 
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Finally, expanding the exponential results in: 
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Only sets of numbers {ck n . } which meet the requirements: 
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have a non-vanishing contribution to the integrals. Therefore one can pull these 
factors in front of the integral: 
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Partition Function 



Reverting the above expansions one returns to the definition of Z(V, P M , Q J ) from 
Refs. [113, 114] times the Boltzmann factors: 



(27T) 



x exp 



(2tt) j 



d 3 p\n 1 ±e ip * a " e 
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-iPi*a L 



±1 



(A.9) 



which proves Eq.(A.l). Therefore, writing for the GCE distribution of extensive 
quantities: 

g-P'X./J e Q*HjP z(V P» O j ) Z p>i ' QJ (VBu u ) 

P 9ce {P^Q 3 ) = * Z[V ^' Q} = 7{ ^Y * ' (A ' 10) 

provides the promised connection between Eqs. (2.3-2.5) and Eq.(2.13). 
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Appendix B 



Second Derivative Test 



Considering the integrand of Eq.(2.13), after change of notation Eq.(2.16), 
one may use short hand notation: 



1(9) 



e~ 1 ^ 1 exp 



V W(6) 



(B.l) 



Taking the first derivative, i.e. the gradient, yields: 
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At the origin, using the definition of the cumulant tensor, Eq.(2.18), one finds: 
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V k[ - Q l 



(B.3) 



The GCE partition function Zgce = exp [V^ (0)] is the th cumulant. Therefore, 
the real part of Eq.(B.3) is identical to zero. In the large volume limit the 
imaginary part is also zero, for Q l = Vk\. Hence, Eq.(B.l) has a stationary 
point at Q\ = 0;. A second derivative test [119, 124] should be done in order to 
distinguish relative minimum, relative maximum or saddle point. Now, taking 
the second derivative, i.e. calculating the Hessian at the origin yields: 
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Second Derivative Test 



Here the imaginary part of Eq.(B.4) is equal to zero, while in the large volume 
limit (choosing intensive variables such that Q l = Vk\) the real part converges 
to: 



dl{6) 



89, d8 r . 



Zgce x Vk™> 1 (B.5) 



h = o ; 

— On 



In order to have a maximum for the integrand 1(6) at the origin, the Hes- 
sian, Eq.(B.5), has to be negative definite. Therefore, k™' 1 has to be positive 
definite. For a symmetric matrix it is enough to demand that all eigenvalues are 
positive. This is the case for all situations considered in this thesis. The deriva- 
tives are proportional to the GCE partition function Zgce- in conclusion, the 
integrand of Eq.(2.13) is strongly peaked at Q\ = 0/ in the thermodynamic limit. 
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Appendix C 



The Cumulant Tensor 



In this section the ft-tensor Eq.(2.18) is examined. Only a static source in 
full acceptance is considered here. Resonance decay effects are included. For 
calculation of asymptotic multiplicity fluctuations and correlations, integrated 
over full momentum space, it is not necessary to take momentum conservation 
into account. The choice of ensemble naturally defines the cumulants needed 
for calculations. Limited acceptance and moving sources, requiring momentum 
conservation, are considered in Refs.[106, 127]. 

The analogs of the matrix «2 and the vector Ki, however in different notation, 
were used in both previously published methods for calculation of scaled variance 
under the thermodynamic limit; the micro-correlator approach, see [93, 95, 99, 
100, 101, 102], and saddle point expansion method [115]. 

Cumulant structure 

Cumulants of order one give GCE expectation values, hence average baryon, 
strangeness, electric charge, and particle density. The second cumulant contains 
information about GCE fluctuations of some quantity (diagonal elements), as 
well as correlations between different quantities (off-diagonal elements). The first 
two cumulants are: 





and 




(C.2) 
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The Cumulant Tensor 



The index I denotes one of L conserved quantities which form the vector Q l . For 
the asymptotic solution to the distributions P(Q l ) the first and second cumu- 
lants are sufficient. The set of conserved quantities could be B, S, Q, E, with 
additionally particle multiplicities Na and N B included, hence: 
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For clarity some elements are explicitly given. A general shorthand notation 
for the derivatives is used: 
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The mean values of the electric charge density and energy density nf are: 
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For instance = in an electric charge neutral system, due to equal contribu- 
tions of particles and anti-particles with — > —q t . The final state mean value k^ a 
of particle density of species A is given by: 
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More about resonance decay and the definition of decay channels T™ A later in this 
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section. In case two quantities are un-correlated, as for example primordial ir 
multiplicity and globally conserved strangeness (n + does not carry strangeness) 
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0. Energy fluctuations k 
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then the corresponding elements k 
correlations between baryonic charge and energy and correlations between 

strangeness and baryonic charge ftf' 5 , are given by: 
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Here for instance k b ' e = for a neutral system, yet in general fvf'' B ^ 0. The cor- 
relations between particle number Na and baryonic charge k b,Na , or energy k e ' Na 
are: 
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(C.13) 



If, for instance Na = N p +N p , then k B ' Na = in a neutral system, as contribution 
of particles and anti-particles cancel each other out, and baryon number conser- 
vation does not affect multiplicity fluctuations of Na = N p + N p . Fluctuations in 
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The Cumulant Tensor 



the density of particle species A, k 
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are: 



K 



N A ,N A 



8 



89 



N A 



8i=0, 



E 



En 



n A 2 
lit 



II A 



(i,i;0) 



4 2,m ■ (C.14) 



Lastly, correlation between the multiplicities Na and N B are only due to resonance 
decay: 
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Resonance Decay 



In this section resonance decay is included analytically directly into the system 
partition function. This has proved far more efficient than the definition of a 
generating function [101, 103], which requires a rather cumbersome calculation 
of all possible primordial correlators. Particle decay is itself a random process. 
Nevertheless, one can assign a particular volume in phase space, given by the value 
of its single particle partition function ipi, Eq.(2.9), to one type of resonance i. 
Resonance decay will now populate this volume in phase space according to some 
weight factor, the branching ratio, for each of the possible decay modes. This 
weight can be assigned to the particle type(s) one is set to investigate. Based 
on the assumption that detected particles are drawn in the form of a random 
sample from all final state particles, e.g. disregarding correlation in momentum 
space, this procedure leads to the acceptance scaling approximation employed in 
Refs.[94, 101, 103]. (see also Appendix E) Conservation laws can be imposed on 
the primordial state (rather than the final state), since decay channels, which are 
experimentally measured, do not only obey charge conservation, but all relevant 
conservation laws (omitting weak decays). 
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For calculation of final state distributions, one needs to determine final branch- 
ing ratios of a resonance type into only stable particles. As an example the 
decay channel A — > B + X is considered with branching ratio 1 Ta^b+x = a,. 
Resonance B could itself be unstable and subsequently decay via the chan- 
nel B — > Y + Z with branching ratio T b->y+z = b. So the final branching ratio is 
defined T a-+x+y+z = ^a^b+x ^b^y+z = oh. Decay tables in [101, 103, 141, 142] 
have also been generated according to this prescription. 

For resonances it seems economical to define further absolute branching ra- 
tios r™ j4 '" s as the sum over all final decay channels of resonance % with given 

2 

numbers ua and % of selected daughters. Hence T { ' is the sum over all final 
decay channels with two daughter particles A and none of species B which are of 
interest, i.e. two positively charged particles in case one wants to calculate u + . As 
a consequence of this definition, branching ratios r" a,Tls will depend on which uj 
or p one is set to calculate. 

For the final state one has to take all absolute decay channels of resonance 
type i into a number ua and of selected stable particles into account. For the 
sake of a common treatment for all particles and resonances are assigned a 'decay' 
channel to stable particles as well, either T]'° = 1 if selected, or Ylm B r°' ns = 1 if 
not selected. The single particle partition function then reads: 

fpi (0e, 0j , On A , 0n b ; Ti A ' nB ) 



where the form of the vector q\ depends on ones choice of ensemble. For in- 
stance, one could have qj = (qi, bi, Sj) for a hadron resonance gas in the CE with 
three conserved charges, or q\ = (qi, bi, Sj, £j) in the MCE without momentum 
conservation. The sum over all decay channels of the selected types needs to be 
one: 

j2 r; A ' nB = i . (c.i7) 

n A ,n B 

This is somewhat of a practical challenge, since decay chains of heavier resonances 

^^The the letter T is used instead of the more conventional Br or br to denote branching 
ratios. 
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The Cumulant Tensor 



are not always well established [162] and respective thermal models codes [88, 89, 
90, 130] struggle to implement this. There are several ways to deal with this, the 
two extreme ones are 1) rescale all known channels according to Eq.(C.17), to 
unity, or 2) assign the missing fraction to the 'channel' r i '°, e.g. to the channel 
without stable particles of interest. 



169 



Appendix D 



Distributions 



From the assumption that the distribution P gce (Q l , N A , Nb) of a GCE with L 
charges Q l and particle numbers of species A and B converges to a multivariate 
normal distribution, it also follows that the marginal distribution P gce {Q l ), a s 
well as the conditional distribution P gce (N A , Nb\Q 1 ), are normal distributions. 
Hence, P m ce{N A , Nb) should have a good approximation in a bivariate normal 
distribution Pbnd(N a , N b ) in the large volume limit (where particle multiplicity 
can be appropriately treated as continuous): 

Pbnd(N a ,N b ) = - Tr - =L= = == (D.l) 



x exp 
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where AN X = N x - (N x ), with X = A, B and: 

Va\ = (N 2 A ) - (N A ) 2 , (D.2) 
Vol = (N B ) - (N B f , (D.3) 
Va AB = (N A N B ) - (N A )(N B ) . (D.4) 

Here Va A and Va B are the variances of the marginal distributions of particle 
multiplicities N A and N B . The term Va AB is called the co- variance. 

The MCE joint multiplicity distribution P mce (N A ,N B ) is conveniently ex- 
pressed by the ratio of two GCE joint distributions, i.e. is given by the conditional 
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probability distribution: 



Pmce(N A , N B ) 



= P gce (N A ,N B \Q l ) 
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In the thermodynamic limit the distributions P gce (NA, N B , Q!) and P gC e(Q!) can 
be approximated by multivariate normal distributions, Eq.(2.29). The extended 
vector Q 7 summarizes particle numbers Na and N B and extensive quantities Q l 
to a vector of length J = 2 + L, hence Q> = (N A , N B , Q ). The vector denoting 
the deviation of the mean values AQ J = Q 3 — Vk{ would then read: 



AQ j = (AN a, AN B , AQ l ) 
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Evaluating the multivariate normal distribution, Eq.(2.29), around its peak 
for (Q l ) yields: 

AQ j = (AN a, AN B , l ) . 
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The vector Eq.(2.26) then becomes: 
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where Ajj are the elements of the matrix Eq.(2.23). Therefore: 

j j j 

t'ti = v 1 [(an a ) 2 J2 ^, 1 + 2(AN a )(AN b )J2 ^,2 + (AN b ) 2 J2 a* 

3 I 3 I 3 I 

(D.10) 

with J = 2 + L = 9, where L = 3 + 4 for a MCE hadron resonance gas with 
momentum conservation. The micro canonical joint multiplicity distribution of 
particle species A and B can thus be written as: 
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where is the 7-dimensional inverse sigma tensor of the distribution P gce {Q 1 )- 
Comparing this to a bivariate normal distribution, Eq.(D.l), one finds: 



£ A « = = A • (EU2) 
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After short calculation one finds for the co-variances: 

< = (D.15) 

= Jg^TcS. (D.17) 
and additionally for the correlation coefficient, Eq.(3.49): 

P = = , , (D.18) 

where the terms A,B,C are given by Eqs.(D.12 - D.14). For the normalization 
in Eq.(D.ll) (from a comparison with Eq.(D.l)) one finds: 
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An analytic formula for the scaled variance, without the need to invert matrices, 
can be found in [105]. 
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Appendix E 



Acceptance Scaling 



To illustrate the 'acceptance scaling' procedure employed in [94, 101, 103] 
uncorrelated acceptance of particles of species A and B is assumed. Particles are 
measured, or observed, with probability q regardless of their momentum. The 
distribution of measured particles ua, when a total number Na is produced, is 
then given by a binomial distribution: 



The same acceptance distribution is used for particles of species B. Independent 
of the original multiplicity distribution P(Na, N b ), the moments of the measured 
particle multiplicity distribution are defined by: 



(n A -n b B ) = n An b BPacc(nA\NA) Pacc(n B \N B ) P(N A ,N B ). (E.2) 
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For the first moment (ua) one finds: 
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The second moment (n 2 A ) and the correlator (tia ■ tib) are given by: 




q 2 (N 2 A ) + q(l-q)(N A ) , 
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For the scaled variance of observed particles one now finds: 



(n 2 A ) ~ (n A ) 2 A 

w> - 1 - q + g ^' 



(E.6) 



where cuf n is the scaled variance of the distribution if all particles of species A are 
observed. Lastly, the correlation coefficient p q is: 
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with (An A AriB) = (n A • n#) — (n A )(riB), and ((An A ) 2 ) = (n A ) — (n A ) 2 . Substi- 
tuting the above relations, one finds after a short calculation: 
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In case ojf n 



47r — cj^, Eq.(E.8) simplifies to: 
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Both lines are independent of the mean values (N A ) and (Nb)- 
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Appendix F 
Convergence Study 

Not only for the sake of completeness, the convergence of various quantities 
with the sample size, i.e. the number of events, N events , in our Monte Carlo scheme 
is discussed in this section. Here only final state (stable against electromagnetic 
and weak decays) particles are analyzed. A closer look is taken at the data sub- 
set of 20 ■ 2 • 10 5 events, with A = V\/V g = 0.875 for the size of the bath, which 
already has been discussed in Chapter 6. 

There is a degree of freedom at so how to estimate the statistical uncertainty 
on the moments of a distribution of observables of a finite sample. The approach 
taken here is straight forward, but could, however, certainly be improved. 

In Fig.(F.l) the evolution of the mean values (N + ) (left) and the vari- 
ances ((AiV + ) 2 ) (right) of the distributions of positively charged hadrons for the 
five transverse momentum bins Apx,i, defined in Table 4.1, with the sample size 
is shown. Mean values of particle multiplicities in respective bins are in rather 
good approximation equal to each other, but are, however, not identical due to 
finite resolution on the underlying momentum spectrum, even for A = 0.875 (bins 
were constructed using GCE events from an independent run). Variances con- 
verge steadily and are different in different bins. The event output was iteratively 
stored in histograms, which were then evaluated after steps of 2 • 10 4 events. 

In Fig.(F.2) the evolution of the scaled variance u + of positively charged final 
state particles (left) and the correlation coefficient p + „ between positively and 
negatively charged particles (right) with the sample size is shown. The results for 
the respective transverse momentum bins can be compared to the second to last 
markers Figs. (6.6,6.7) (left), which denote the corresponding results of grouping 
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Figure F.l: Step histogram showing the convergence of the mean values (N + ) 
(left) and variances ((AiV + ) 2 ) (right) for positively charged final state hadrons in 
transverse momentum bins ApT,i for a hadron resonance gas with A = V\/V g = 
0.875. 
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Figure F.2: Step histogram showing the convergence of the scaled variance u + 
(left) of positively charged hadrons and the correlation coefficient p + _ be- 
tween positively and negatively charged hadrons (right) in transverse momentum 
bins Ap Tji for a final state hadron resonance gas with A = Vi/V g = 0.875. 



the same data into 20 Monte Carlo sets of 2 • 10 5 events each. 

The distribution of scaled variances of positively charged particles u + (left) 
and correlation coefficients between positively and negatively charged parti- 
cles p_| (right), resulting from grouping again the same data set into 200 samples 

of 2- 10 4 events each are shown in Fig.(F.3). The transverse momentum bin Apr,5 
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Figure F.3: Histogram showing the results for the scaled variance u + (left) of 

positively charged hadrons and the correlation coefficient p H between positively 

and negatively charged hadrons (right) in the transverse momentum bin Aj»t,5 



for a final state hadron resonance gas with A 
runs of 2 • 10 4 events each are analyzed. 



V 1 /V g = 0.875. 200 Monte Carlo 
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Figure F.4: Evolution of the primordial scaled variance u + of positively charged 
hadrons (left) and the primordial correlation coefficient p + _ between positively 
and negatively charged hadrons (right) with the Monte Carlo parameter A = 
Vi/Vg in different rapidity bins A?/j. The solid lines show an analytic extrapolation 
from GCE results (A = 0) to the MCE limit (A — > 1). The 5 leftmost markers 
and their error bars represent the results of 20 Monte Carlo runs of 2 • 10 5 events. 
Three additional values of A have been investigated with 20 Monte Carlo runs 
of 1 • 10 7 events. The 5 rightmost markers denote the results of the extrapolation. 
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for a final state hadron resonance gas with A = V\/V g = 0.875 was chosen. 

Monte Carlo results for A = 0.875 of the analysis shown in Fig.(F.3), are 
for the scaled variance oj + (Apr 5 ) = 0.8069 ± 0.0514, and the correlation coef- 
ficient p + _(Apr 5 ) = —0.0026 ± 0.0421. They are nicely scattered around the 
mean values, denoted by the bottom lines in Fig.(F.2), cj + (Apr j5 ) = 0.8082, 
and p + _(Ap T5 ) = —0.0028, respectively. 

They are also compatible with the analysis shown in Figs. (6. 6, 6. 7), of Chap- 
ter 6, u + {Ap T)5 ) = 0.8081 ± 0.0149, and p + _(ApT i5 ) = -0.0022 ± 0.0125, at 
the same value of A. The comparatively large statistical error on the analysis 
in Fig.(F.3) is due to the splitting up into many small sub-samples. The mean 
values of different analysis agree rather well. 

Lastly, in Fig.(F.4) the results of additional Monte Carlo runs are shown for 
values of A closer to unity. Here 20 additional runs of 1 • 10 7 primordial events 
for A = 0.925, 0.950, and 0.975 were performed. As discussed above, error bars 
diverge, but convergence seems to be rather good. The additional data has not 
been used for the extrapolation, so it can serve as an un-biased cross-check. 
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Appendix G 

The Canonical Boltzmann Gas 



An analytical and instructive example is the canonical classical relativistic 
particle anti-particle gas discussed in [93, 94, 159, 160]. This example is used 
to show that, although the procedure is formally independent of ones choice of 
Lagrange multipliers, it is most efficient for those defined by Maxwell's relations. 
Starting off with Eqs.(3.1), this section will discuss, in turn, the first and second 
moments of the multiplicity distribution of particles, and the first four moments 
of the Monte Carlo weight factor. 

The canonical partition function Z^^Vi, /3, Q\) of a system with volume V±, 
temperature T = charge Qi, particle number Ni, and anti-particle num- 
ber Mi = iVx — Qx, is given by: 

Z (V 8 0) - {Vl ^ )Nl ^f^l (G u 

Zm(Vi,P,Qi) - -^r~ {Nl _ Ql)] ■ ( ) 

The single particle partition function in Boltzmann approximation is given by 
Eq.(3.31), ip = ^2 m2 K 2 (m/3). The canonical partition function with arbi- 
trary particle number, but still fixed charge Q\, is obtained by: 

oo 

ZQfuPM = ZnAVu^Qi) = IqAZV^) . (G.2) 

A r i=Qi 

Here Iq 1 is a modified Bessel function. Temperature is the same in both subsys- 
tems; the bath and the observable part. The partition function of the bath is 
therefore: 

Z(V 2 ,(3,Q 2 ) = I Q2 (2V 2 ^) . (G.3) 
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Imposing the constraints V 2 = V g — Vi, and Q 2 — Q g — Qi, similar to Eq.(3.4), one 
finds [119] for the canonical partition function, Eq.(3.5), of the combined system: 

oo 

Z(V g ,P,Q g ) = T Qt( 2V ^) ^0.-0i (2 - Vi) ^) = I Qg (2V 9 i>) , 

Qi=— 00 

(G.4) 

as required. The weight factor is then: 

L Q g \ AV gV) 

Analogous to Eq.(3.7) one obtains for the joint particle multiplicity and charge 
distribution: 

P(Qi,7Vi) = WfaQxVvQglP) Z Nx {Vx,PM . (G.6) 

Monte Carlo Weight 

Next, introducing Eq.(3.12), the joint GCE distribution of charges and particle 
multiplicity: 

PgceiQuNj = - Z Nl (Vi,/3,Qi) . (G.7) 

Z{Vi, 

The Monte Carlo weight, Eq.(3.15), is then given by: 

VV^O/^IA/i) = W(Vi,Qi;Vg,Q g \{3) Z(Vi,f3,n) e~ Ql ^ . (G.8) 

In accordance with Eq.(3.11), the distribution Eq.(G.6) is then equivalently writ- 
ten as: 

PiQuNj = W 0li0 «(Vi;V>,/3) PgceiQuNj . (G.9) 
The GCE partition function is: 

Z{Vx,P,fi) = eQl ^ z (Vuf3,Qi) = exp^cosh^/i)] . (G.10) 

Qi=— 00 
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Moments of Distributions 

The moments of the multiplicity distributions Eq.(G.6) or Eq.(G.9) are given 



by: 



oo oo 



( N i) = E E N ? p (Ni,Qi 

N 1= Qi=-oo 



(G.n; 



Additionally, the moments of the weight, Eq.(G.5), are defined through: 



oo oo 



(W n ) = E [w(V 1 ,Q 1 ;V g ,Q g \P)] n Zn^^Q,) , (G.12) 



JV 1= 0Qi=-oo 

and of the Monte Carlo weight, Eq.(G.£ 



oo oo 



(W n ) = J2 E [^ Ql ' QB (Vi;V g \(3^)\ PgceiQ^Nx). (G.13) 

AT 1= 0Qi=-oo 

Next, attending to the first two moments of the multiplicity distribution, substi- 
tuting Eq.(G.6), or Eq.(G.9) into Eq.(G.ll) yields: 



and 



(G.14) 



(G.15) 



Canonical suppression of yields and fluctuations acts on the global volume V g . In 
the GCE the first two moments are (Ni) = V^e^, and (Nf) = (iVi) 2 + (Ni), 
respectively. The CE limit is obtained by V g — > V\, and Q g = (Qi). Substituting 
Eq.(G.14) and Eq.(G.15) into Eq.(3.50), and using Eq.(3.9), A = VjV g , yields: 



uj = A uj ce + (1 — A) oj 



gee i 



(G.16) 



where the CE scaled variance u ce of the combined system is given by [93, 94]: 



i - ■ (y g i>) 



/ qg -i (2V g rJ>) _ I Qa - 2 (2V g ^) 
I Qg (2V g ^ I Qa -, (2V g ^) 



(G.n) 



and u gce = 1 is the GCE scaled variance, as the particle number distribution is a 
Poissonian. 
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The Monte Carlo scheme is now applied to an observable subsystem of vol- 
ume V\ = 50 fm 3 embedded into a system of volume V g = 75 fm 3 , charge Q g = 10, 
and temperature T = = 0.160 GeV. Particles and anti-particles have 

mass m = 0.140 GeV and degeneracy factor g — 1. The average charge content in 
the observable subsystem is then (Qi) — 6.667. The mean particle multiplicity, 
Eq.(G.14), is (Ni) ~ 7.335, and the scaled variance of particle number fluctua- 
tions, Eq.(G.16), is oj ~ 0.3896. The GCE in V\ is sampled for various values 
of \iq. The Monte Carlo weight, Eq.(G.8), is then employed to transform these 
samples to have the statistical properties required by Eq.(G.6) or Eq.(G.9). For 
each value of \iq fifty samples of 2000 events each have been generated to allow 
for calculation of a statistical uncertainty estimate. 




0.05 0.1 0.15 0.2 0.25 0.3 0.05 0.1 0.15 0.2 0.25 0.3 



Figure G.l: The first four moments of the Monte Carlo weight, Eq.(G.8) (left) 
and the first two moments of multiplicity distributions (right) , as described in 
the text. 

In Fig.(G.l) (right) mean value (Ni) and variance ((AiVi) 2 } of the particle 
multiplicity distribution of the original GCE samples for different values of chem- 
ical potential \xq are shown in open symbols. The closed symbols denote mean 
value and variance of these samples after the transformation Eq.(G.8) was ap- 
plied. Independent of the original sample the result stays (within error bars) the 
same. However, the statistical error is lowest for a chemical potential close to: 



(G.18) 
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i.e. when the initial sample is already similar (at least in terms of mean values) to 
the desired sample. This is reflected in the moments of the Monte Carlo Weight 
factor, Fig.(G.l) (left). Higher moments have a strong minimum around fiQ = 
0.1896 GeV, i.e. the weights are most homogeneously distributed amongst events, 
and most efficient used is made of them. 
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